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Abstract. We consider triangular holes on the hexagonal lattice and we study their inter- 
action when the rest of the lattice is covered by dimers. More precisely, we analyze the joint 
correlation of these triangular holes in a "sea" of dimers. We determine the asymptotics of 
the joint correlation (for large separations between the holes) in the case when one of the holes 
has side 1, all remaining holes have side 2, and the holes are distributed symmetrically with 
respect to a symmetry axis. Our result has a striking physical interpretation. If we regard 
the holes as electrical charges, with charge equal to the difference between the number of 
down-pointing and up-pointing unit triangles in a hole, the logarithm of the joint correlation 
behaves exactly like the electrostatic potential energy of this two-dimensional electrostatic 
system: it is obtained by a Superposition Principle from the interaction of all pairs, and the 
pair interactions are according to Coulomb's law. 



1. Introduction 

Monomer-monomer and especially dimer-dimer correlations 1 on a plane bipartite lattice 
(especially the square and hexagonal lattice) have been studied quite extensively (see for 
instance [12], [17], [Kl] and [K2]). Color the vertices of the lattice white and black so 
that each edge has one white and one black endpoint. From the point of view of this paper, 
there is a fundamental difference between studying dimer-dimer and monomer-monomer 
correlations: the former have the same number of vertices of each color, while the latter 
have an excess of either a white or a black vertex. 

In this paper we consider correlations of triangular holes on the hexagonal lattice. 

This type of hole has the convenient feature that the difference between the number of 
its white and black constituent vertices is equal to the length of its side. We will be lead 
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1 The monomers, respectively the dimers, are interacting via a sea of dimers that cover all lattice sites 
not occupied by them. 



by our results to interpret the white vertices as elementary negative charges ("electrons"), 
and the black vertices as elementary positive charges ( "positrons" ) , so that the triangular 
plurimers become charges of magnitude given by their side-length, and sign given by their 
orientation (up-pointing or down-pointing). The main result of this paper, Theorem 2.1 
(see also its much simpler restatement (2.6)), implies that in the fairly general situation in 
which it applies (namely, when the holes are symmetrically distributed about an axis, and 
all holes have side 2, except for one of side 1, on the symmetry axis) the logarithm of the 
joint correlation of triangular holes behaves exactly like the two-dimensional electrostatic 
potential energy of the corresponding system of charges: it is obtained by a Superposition 
Principle from the interaction of all pairs, and the pair interactions are according to 
Coulomb's law. (It is now clear why dimer-dimer and monomer-monomer correlations are 
fundamentally different: dimers are neutral!) 

To present our results in the background of previous related results in the literature, 
we point out the following facts. 

First, we mention that there is an alternative approach for expressing joint correlations 
of holes on the hexagonal (or square) lattice due to Kcnyon ([21, Theorem 2. 3]; see also 
[20]). When it applies, it provides an expression for the joint correlation as a k x k 
determinant (where 2k is the total number of vertices in the holes). However, while not 
requiring symmetry, the set-up of Kenyon's approach limits its applicability to the case 
when all holes have even side (thus not accommodating our hole of side 1 on the symmetry 
axis), and, more restrictively, to the case when the total "charge" of the holes is zero 2 . 
The main advantage of our approach is that it sets no restriction on the the total charge. 

Furthermore, our result (14.9) gives the joint correlation of a general distribution of 
collinear monomers on the square lattice, a situation in which none of the holes satisfies 
the even-sidedness required by Kenyon's approach. 

Second, there are other discrete models in the physics literature (see e.g. the survey 
[27] by Nienhuis) believed to behave like a Coulomb gas. However, there are several im- 
portant differences between them and our model: (1) we do not require that the total 
charge is 0, a fact built into the definition of the Coulomb gas model in the survey by 
Nienhuis [(2.7), 27]; (2) by studying correlation of holes, our discrete model seems quite 
different from the others in the literature; (3) for those models surveyed in [27] for which 
the (believed) Coulomb behavior is only asymptotic (as it is the case for our model), the 
arguments for their equivalence with the Coulomb gas model are only heuristic, while our 
results are proved rigorously; (4) in our model all states have the same energy, so the 
emergence of the Coulomb interaction is entirely due to the number of different geomet- 
rical configurations (unit rhombus tilings) compatible with the holes — in the language of 
physicists, our model is stabilized by entropy; by contrast, in all models surveyed in [27] 

2 Indeed, Kcnyon [20] defines the correlation of holes as the limit of M{Hm tn )/M{H m ^ n ) when m, n — + 
oo, where H m ^ n is a toroidal hexagonal grid graph, and H m ,n is its subgraph obtained by removing the 
holes (M(G) denotes the number of perfect matchings of the graph G). In order for M(H myn ) to be 
non-zero one must have the same number of vertices in the two bipartition classes of the vertices of H my „. 
Thus, since M(H m ,n) also has to be non-zero in order for Kenyon's correlation to be non-zero, the number 
of vertices in the two bipartition classes that fall in the holes must also be the same; i.e., the total charge 
of the holes must be zero. 
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different configurations have different energies, specified by a Hamiltonian — they are sta- 
bilized by energy. Furthermore, in our model the electrical charge has a purely geometric 
origin: the charges are holes in the lattice, and their magnitude is the difference between 
the number of right-pointing and left-pointing unit triangles in the hole. 

Third, in a recent paper [22] Krauth and Moessner study numerically (using Monte 
Carlo algorithms) a very special case of the problem considered in this paper, namely just 
the two-point correlations of monomers on the square lattice. Their data leads them to 
conjecture that the two-point correlations on monomers behave like a Coulomb potential 
(the case of monomers on different sublattices appeared already in [12]; the new part is 
the simulation data for monomers on the same sublattice). Krauth and Moessner also 
state that they could not find these correlations worked out in the literature. Since in 
the current paper we address the case of (2m + 2n + l)-point correlations, showing that 
they satisfy the much more general Superposition Principle, their remark suggests that 
our results are also new. 

And fourth, shortly after the current paper was posted on the preprint archive (web 
address arxiv.org/abs/math-ph/0303067, March 2003), physicists D. A. Husc, W. Krauth, 
R. Moessner and S. L. Sondhi posted an article (arxiv.org/abs/cond-mat/0305318, May 
2003) presenting numerical simulations that suggest a positive answer to Question 15.1 
(which concerns a three dimensional analog of the model presented in this paper) for the 
case of two monomers. 

2. Definitions, statement of results and physical interpretation 

We will find it more convenient to present our results on dimer coverings of portions of 
the hexagonal lattice in the equivalent, dual language of lozenge tilings of lattice regions 
on the triangular lattice. In this language a monomer is just a unit triangle of the lattice. 
A lattice triangle of side two will be called a quadromer. The results proved in this paper 
involve only monomers and quadromers. More generally, triangular plurimers are lattice 
triangles of arbitrary side. A dimer becomes a lozenge — a unit rhombus covering precisely 
two unit triangles. We will often refer to lozenges as dimers. 

Let m and n, and N be nonnegative integers. Consider also the nonnegative integers 
Ri, V{, i = 1, . . . , m and R^, v[, i = 1, . . . , n. Define the region 

/ Ri R m R[ R' n \ 
V Vl v m v 1 v n ) 

as follows. 

Consider a lattice hexagon H whose sides alternate between 2N + An + 1 and 2N + 4m, 
starting with the base. Denote by I the vertical symmetry axis of H . Let u be the up- 
pointing unit triangle on I whose base is (N + 2m)v / 3 units above the base of H (that is, 
the base of u is along the horizontal diagonal of H). 

Let D(R, v) be the down-pointing quadromer (i.e., down-pointing lattice triangle of side 
2) whose base is centered R units to the left of £, and (2v + l)\/3/2 units below the base of 
u (two instances of this appear in Figure 2.1). Let U(R, v) be the up-pointing quadromer 
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(2.1) 



whose base is centered R units to the left of I and (2v + l)V3/2 units above the base of u 
(see Figure 2.1 for three instances of this). 

Finally, let D(R,v) and U(R,v) be the mirror images in £ of the above quadromers. 

We define H N ■ ■ ■ m ; / ■ ■ ■ to be the region obtained from the hexagon H 

\ Vl V m V 1 v n J _ 

by removing u (a monomer), D(Ri,Vi), D(Ri,Vi), i = 1, . . . ,m and U(R'j,v'j), U(R'j,v'j), 
j = 1, . . . , n (a total of 2m + 2ra + 1 holes). 

We assume throughout the paper that Ri > 1, i = 1, . . . , m, and ii^ > 1, j = 1, . . . , n, 
so that I separates the original plurimers from their mirror images. Figure 2.1 shows the 

regions ^ ^ 2 4 j. 

Define the correlation at the center (or simply correlation) of these 2m+2n + l plurimers 

by 
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(2.2) 



where M(D) is the number of dimer coverings (i.c, lozenge tilings — tilings by unit rhom- 
buses each covering precisely two unit triangles) of the lattice region D. The existence of 
the limit on the right hand side of (2.2) follows from Proposition 3.2. 

In the above definition, the region at the denominator is obtained from the one in the 
numerator by packing the quadromers as tightly as possible around u — for the example 
of Figure 2.1 this is illustrated in Figure 2.2. As indicated by the latter figure, because of 
forced lozenges, the situation is equivalent to removing just two large plurimers from H, 
one down-pointing of side 2m and one up-pointing of side 2n + 1. 

Assume that the midpoints of the bases of the quadromers have coordinates 



Ri = AiR 

Vi = q z Ri + q = qiAiR + c l 
Rj = BjR 

v'j = q'jR'j + c'j = q'jBjR + c' p (2.3) 



where: Ai > 0, qi > 0, i = 1, . . . , m are fixed rational numbers chosen so that the pairs 
(Ai, q\), . . . , (A m , q m ) arc distinct; Bj > 0, q'j > 0, j = 1, . . . , n are fixed rational numbers 
so that the pairs (B\, q[), . . . , (B n , q' n ) are distinct; Cj > 0, i = 1, . . . , m and c'j > 0, 
j = 1, . . . , n are fixed integers; and R is an integer parameter. 
The main result of this paper is the following. 
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Theorem 2.1. If the coordinates of the quadromers are as in (2.3) and R — > oo, t/ieir 
joint correlation with the fixed monomer u is given asymptotically by 

Rl R-m R'l R' n 

vi v m ' v[ v' n 

c f\(2R^T\(2R' ^=^ ? + ^ ? 

2. j ^ J- J- / — 1 ? 1 

x II [(^J -^i) 2 + 3(uj - v,) 2 } 2 [(Rj + Rif + 3(v 3 - v,f 

l<i<j<m 

YlKi^Wj ^) 2 + 3(4 - ^) 2 ] 2 [Ogj + ^) 2 + 3 (4 - 4) 2 ] 2 

x uti n?=i[(^ - ^) 2 + 3(4 + ^) 2 ] 2 [(^- + i?o 2 + 3K- + «. t ) 2 p 

+ ^ R A(m-nf-4 m -l^ ^ 



\2i2 



where 



Cfc ; 



n(i)^n-(372)r) • (2 - 5) 



Remarkably, each factor in (2.4), except for the multiplicative constant c 2m ,2 n , is exactly 
equal to the Euclidean distance between the midpoints of the bases 3 of some pair of our 
2m + In + 1 plurimers (one monomer and 2m + 2n quadromers). Indeed: 

(i) 2Ri is the distance between D(Ri,Vi) and D(Ri,Vi); 

(ii) 2R' J is the distance between U(R'j,v'j) and U(R'j,v'j); 

(Hi) [(R'j) 2 + 3(4) 2 ] 1 ^ 2 i s both the distance between U(Rj,v'j) and monomer u, and 
the distance between U(R'j,v'j) and monomer u; 

(iv) [R 2 + ivf] 1 / 2 is both the distance between D(Ri,Vi) and monomer u, and the 
distance between D(Ri, v{) and monomer u; 

(v) [(Rj — Ri) 2 + 3(vj — vi) 2 ] 1 / 2 is both the distance between D(Ri, vi) and D(Rj,vj), 
and the distance between D(Ri,Vi) and D(Rj,Vj); 

(vi) [(Rj + Ri) 2 + 3(vj — Vi) 2 } 1 / 2 is both the distance between D(R i7 Vi) and D(Rj,vj), 
and the distance between D(Ri,Vi) and D(Rj,Vj); 

(vii) [(R'j - R' t ) 2 + 3(4 - 4) 2 ] 1/2 is both the distance between U(R[, v[) and U(R' j ,v' j ), 
and the distance between U(R'i,v'i) and U(R'j,v'j); 

(viii) [(R' 3 +R'i) 2 +Z(v' J -v' i ) 2 ] 1 / 2 is both the distance between U(R' i ,v' i ) and U(R' j ,v' j ), 
and the distance between U(R'i,Vi) and U(Rj,v'j); 

(ix) [(R'j - Ri) 2 + 3(4 + ^) 2 ] 1/2 is both the distance between D(R l ,v i ) and U(R' J ,v' j ), 
and the distance between D(Ri,Vi) and U(R'j,v'j); 



3 For the monomer u — involved in the factors of the fraction on the second line of (2.4) — instead of the 
midpoint of its base we need to consider the closest lattice point above it or below it, according as the 
factor in question is at the numerator or denominator, respectively. 
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0) [(R'j + Rj) 2 + 3(v'j + Vi) 2 ] 1 ' 2 is both the distance between D(R ll v i ) and U {R'^v'^), 
and the distance between D(Ri, Vi) and U(Rj, v'j). 

This remarkable phenomenon goes even further: if one defines the charge ch(Q) of a hole 
to be the number of its up-pointing unit lattice triangles minus the number of its down- 
pointing ones (this clearly gives 1 for the monomer u, and —2 and 2 for the quadromers 
of type D and U, respectively), the exponent with which each such factor occurs in (2.4) 
is precisely half the product of the charges of the corresponding two plurimers. Therefore, 
the statement of Theorem 2.1 can be rewritten as follows. 

Theorem 2.1 (equivalent restatement). Denote the (2m + 2ra+l) triangular pluri- 
mers removed in region (2.1) by Q\, . . . , Qi m +in+\> a nd define the charge ch(Q) of Q to be 
the number of up-pointing unit lattice triangles of Q minus the number of its down-pointing 
unit triangles. Let d(Qi,Qj) denote the Euclidean distance between the midpoints of the 
bases of plurimers Qi and Qj . 

Then as the coordinates of the midpoints of the bases of the 2m+2n quadromers approach 
infinity as specified by (2.4), the asymptotics of their joint correlation with the monomer 
u (which by definition is kept fixed at the origin) is given by 

U ) (Ql,...,Q 2 m+2n+l) = C2m,2n J] d(Qi, Q j ) ch < '> dl «'>/ 2 

l<i<j<2m+2n+l 

+0{R 4{m - n)2 - Am - 1 ), (2.6) 

where the constant Ck,i is given by (2.5). 

Our result has the following striking physical interpretation. 

Take the logarithm of both sides of (2.6), and assume 4(m — n) 2 — Am ^ 0. Since then 
the main term on the right hand side of (2.4) approaches either or infinity as R — > oo, 
the resulting contribution from the multiplicative constant is negligible. We obtain 

log W (Q 1) ...,Q 2ro+2n+1 )~ Ch(0i) 2 Ch(Qj) logd(Q„Q 3 ), (2.7) 

l<i<j<2m+2n+l 

where d(Qi, Qj) is the Euclidean distance between the plurimers Qi and Qj. 

But the logarithm of the distance is, up to a negative multiplicative constant, just 
the potential in two dimensional electrostatics! So if we view the plurimers as electrical 
charges of (signed) magnitude given by the operator ch (in this case, just their side- 
length if they point upward, or minus their side-length if they point downward), the 
logarithm of the joint correlation (multiplied by a negative constant) acts precisely like the 
electrostatic potential energy of this two dimensional electrostatic system: it is obtained 
by a Superposition Principle from the interaction of all pairs, and the pair interactions are 
according to Coulomb's law. 

Both for our proof of Theorem 2.1 and from the point of view of physical interpretation 
we find it useful to define another kind of plurimer correlation, in which the boundary of 
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Figure 2.3. W 2 J; * g 4)- 



the defining regions makes its presence felt. To this end, define another family of regions 
as follows. 

(R R R f R f \ 
1 • • • m ; } ■ ■ ■ ? ) along a zig-zag line that starts at the 
ui w m ' v[ v' n J 

midpoint of its top side, follows t as closely as possible on its right until it reaches the 

central monomer u, crosses £ along the western edge of u and then follows I as closely as 

possible on its left, ending at the midpoint of the base (this cut is pictured in Figure 2.3). 

We define 

Rl Rm R'\ R' n \ 

vi v m ' v[ v' n J 

to be the region obtained this way to the west of our cut, in which in addition the N + 2n 
lozenge positions above u and immediately to the left of the cut are distinguished and given 
weight 1/2 (i.e., each lozenge tiling T of the regions Wjv is weighted by l/2 k , where k is 
the number of distinguished lozenge positions occupied by a lozenge in T). An example 
is illustrated in Figure 2.3 (the distinguished lozenge positions are marked by shaded 
ellipses) . 

When N — > 00 and the coordinates Ri, Vi, R[, v[ of the quadromers are fixed, they 
maintain a fixed relative position with respect to the right boundary of the regions Wjv. 
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We define the boundary-influenced correlation u>b by 



Rl Rm R'l R'n 

V-L V m ' v[ V' n 



M W N 



Rl i?2 Rm . R'l R'2 R'n 

V-L V 2 V m ' V[ v' 2 V n 



Alfu-vf ' :i ... 2 " i '•' ') .2.- 

' 1 '00 

(the fact that this limit exists follows by Lemma 5.1). 

An important part of the proof of Theorem 2.1 — and a result interesting on its own — is 
the determination of the asymptotics of LOb- 

Theorem 2.2. If the coordinates of the quadromers D(Ri, Vi), i = 1, . . . ,m and U(Rj,Vj), 
j = 1, . . . , n are given by (2.3), the large R asymptotics of their boundary-influenced cor- 
relation is 



where 



Rl Rm . R'l R'n \ 

vi v. m ' v'i v' n J ~ 

™ " / n"=i[(^-) 2 + 3(^) 2 ] 1/2 

hm ' 2n l\ {2Rj) YlZARU^Y/i 

x II [(Rj - Rif + 3fe - Vi) 2 } [(Rj + Ri) 2 + 3( Vj - v,) 2 } 

l<i<ij<rn 

n i<t<3 <J(^- - Rtf + 3Qj - vj) 2 } [{R' 3 + Ri) 2 + Sty - v[) 2 ] 

x uz 1 n;=i m + 3 K- + «o 2 ] + ^) 2 + 3 K + ^) 2 ] 

+ 0( - jR 2(m-™) 2 -2m-l^ ^g) 



2 fc 3 (fc+o/4-(fe-o 2 /4 fe -i ( 2 )j ^(j + 2)fc 



/i (l)i(3/2), fJi (3/2), ' 

This result also has an interesting physical interpretation. Relabel the m+n quadromers 
by Qi, . . . , Qm+n- Assuming 2(m — n) 2 — 2m 7^ 0, one deduces from (2.9) — just as (2.7) 
was deduced from (2.4) — that 

log £ ,; 6 (Q 1 ,...,Q ro+n )~ £ ch(Q)ch(Q0 logd(Q QOi (2n) 

Q,Q'eQ,Q^Q' 

where ch denotes the charge, d the Euclidean distance, and Q contains, in addition to 
our m + n plurimers, their mirror images in the vertical line touching all zig-zags on the 
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right boundary of Wjv, and also an up-pointing monomer just outside the lower of the two 
aligned consecutive segments of this right boundary. 

Relation (2.11) also has a strong electrostatic reminiscence. It shows that the quadro- 
mers near a zig-zag boundary can be interpreted as electrical charges near a straight line 
whose effect is to bring in the mirror images of our charges. A well-known situation in 
which this happens is for electrical charges near a conductor — but in that case the image 
charges need to be taken of opposite sign. In fact, there is another physical situation in 
which the images are taken with the same sign: that when we have fixed charges inside a 
dielectric and near a straight boundary separating this dielectric from another one, of a 
dielectric constant negligible in comparison with the first (see e.g. [11, §12-2]). Our relation 
(2.11) models this situation — with the additional specification that the slight irregularity 
of the right boundary near its middle behaves like an extra charge of magnitude +1. 

We conclude this section by presenting a view of our results as a possible random tiling 
model for two dimensional electrostatics. 

Suppose one pictures the two dimensional "universe" being spatially quantized, con- 
sisting of a very fine lattice of triangular "quanta of space", each of them responsible 
for creating a unit charge (of sign given by its orientation) when displaced by a "body" 
(which, since we are considering space to be quantized, can be placed only in such a way 
that it consists of whole unit triangles). Suppose also that once a number of bodies are 
placed in this space — thus acquiring the charge determined by the above interpretation — 
the unoccupied quanta of space (i.e., unit triangles) have the tendency to pair up with 
their neighbors to form a sea of dimers, according to the uniform distribution on all such 
possible pairings. (Assuming such a tendency for unoccupied unit triangles to pair up 
with a neighbor seems natural from the point of view of quantum mechanics, as the quan- 
tum fluctuations of the vacuum ceaselessly cause particle-antimatter-companion pairs to 
erupt into existence only to be annihilated after an instant. One particular rhombus tiling 
corresponds to one particular way for these virtual particles to annihilate. To consider av- 
eraging over all such possible tilings is analogous to Feynman's "sum-over-paths" approach 
to quantum mechanics; see the remark at the end of Section 15.) Then (2.7) shows — under 
the assumption that the charges are symmetrically distributed and are all of magnitude 
±2, except for a single +1 on the symmetry axis — that these charged bodies would interact 
precisely like charged bodies in two dimensional electrostatics. We would get a proof of 
the Superposition Principle (for such charge distributions), and very strong evidence for 
Coulomb's law! 4 

To make this parallelism more explicit, let S be the family of all (2m + 2n + l)-element 
sets of triangular plurimers (one monomer, and the rest quadromers) satisfying the hy- 
pothesis of Theorem 2.1, and being contained in some fixed disk D centered at the center 
of the regions (2.1) (the radius of D is some large absolute constant, the "diameter" of 
our system of plurimers). 



4 The reason we don't obtain a proof for Coulomb's law from Theorem 2.1 is that we are assuming 
there a symmetric charge distribution, and therefore the case of two charged bodies in general position 
is not covered. This case is addressed in [6] for opposite sign charges of magnitude 2. The general 
even-magnitude case will be presented in a sequel of the present paper. 
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Let A = {Ax, . . . , A 2m +2n+i} and B = {Bi, . . . , B 2m +2n+i} be two members of S. Re- 
denote the charges of their elements by ch(Aj) = ch(Bj) = Qi, i = 1, . . . , 2m + 2n + 1, and 
denote their corresponding regions (2.1) by Hn(A) and H^{B), respectively. Denote also 
their corresponding correlations (2.2) by lu(A) and u>{B), respectively. 

Define a probability distribution on S by requiring the ratio of the probabilities Pjs, and 
Pb to be 

EA = y m ( h n(A)) _ to(A) 
Pb 



:= lim 

iV->oc 



M(H N (B)) u(B)' 



(2.12) 



the second equality being valid by (2.2). 
By (2.6) we obtain from (2.12) that 



Pa 
Pb 

exp I 



53 ^(-lnd^,^)) - 53 q l q 1 (-\nd(B l ,B j )) 

l<i<j<2m+2n+l l<i<j<2m+2n+l 



(2.13) 



for large mutual distances between the elements of A and B. 

Consider, on the other hand, a two dimensional physical system of 2m + 2n + 1 charges 
Qi, ■ ■ ■ , Q2m+2n+i of magnitudes q\q e , . . . , q2m+2n+iq e (expressed as integer multiples of 
the elementary charge q e ). Then assuming that there are only electrical forces between 
them (governed by the two dimensional electrostatic potential, — l/(27reo) ln(rf), where d 
is the distance and eo is the permittivity of empty space [11, §4-2, §5-5]), we obtain by 
the Fundamental Theorem of Statistical Physics (see e.g. [10, §40-3]) that the probability 
P e (d) of finding the charges at mutual distances d = (dij)i<i < j<2 m +2n+i, relative to the 
probability P e (d') of finding them at mutual distances d' = (rf^)i<i<j<2m+2n+i, is 

Pe(d') ~ 



exp 



2ire kT 



53 9iQj{- 

l<i<j<2m+2n+l 



Indij)- 53 H1j(- lTLd 'ij) 

l<i<j<2m+2n+l 



(2.14) 



where k is Boltzmann's constant (see e.g. [10, §39-4]) and T is absolute temperature. 

Relations (2.13) and (2.14) show that, at least in the case when the two dimensional 
physical system of electrical charges satisfies the magnitude and geometrical distribution 
requirements in the hypothesis of Theorem 2.1, their physical electrostatic interaction at 
temperature T = ql/(Treok) is correctly modeled by our model 5 . We present in Section 15 
a possible three dimensional analog of this parallelism that allows any temperature T. 



5 In fact, if one accepts Conjecture 14.1, our set-up allows one to introduce an extra parameter in 
(2.13), a positive integer x, which makes the parallelism between (2.13) and (2.14) go through for any 
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Moreover, what plays the role of electrostatic potential in our model is just the entropy — 
the logarithm of the number of dimer coverings — which had been considered before, in a 
different context [2] [23] (it involved the number of dimer coverings of the inside structure 
of a molecule, not of the outside, as in our situation), as a possible measure for chemical 
energy (one of the original motivations was simply its satisfying the addition principle, for 
isolated systems). 

Our interpretation has the advantage that Coulomb's law emerges as an asymptotic 
law valid when the distances between the charges are large in comparison with the size of 
the charges. This would explain why Coulomb's law seems to break down at very small 
distances (cf. [11, §5-8], at distances lower than about 10~ 14 ). 

Furthermore, in this context (2.11) can be interpreted as describing the interaction of 
electrical charges with an "edge" of this two dimensional "universe." The same sign for 
the image charges ensures a repelling effect, and this helps keeping the charges inside the 
"universe" ! 

An appealing feature of such a model is that it is discrete, and therefore models the 
physical electrostatic interaction by considering a discrete ambient space. 

The possibility of a discrete "machinery" underlying the electrostatics of the real phys- 
ical world is mentioned by Feynman in [11, §12-7]. When indicating how several equations 
in physics, like for instance that for neutron diffusion, are true only as approximations 
when the distance over which one looks is large (for neutron diffusion, large in comparison 
to the mean free path), Feynman goes on to ask: 

"Is the same statement perhaps also true for the electrostatic equations? Are they also correct only as 
a smoothed-out imitation of a really much more complicated microscopic world? Could it be that the real 
world consists of little X-ons which can be seen only at very tiny distances? And that in our measurements 
we are always observing on such a large scale that we can't see these little X-ons, and that is why we get 
the differential equations? 

Our currently most complete theory of electrodynamics does indeed have its difficulties at very short 
distances. So it is possible, in principle, that these equations arc smoothed-out versions of something. 
They appear to be correct at distances down to about 10~ 14 cm, but then they begin to look wrong. 
It is possible that there is some as yet undiscovered underlying 'machinery,' and that the details of an 
underlying complexity arc hidden in the smooth-looking equations — as is so in the 'smooth' diffusion of 
neutrons. But no one has yet formulated a successful theory that works that way." 

The random tiling model presented in this paper (and pursued further in a subsequent 

temperature. To this end, refine our triangular lattice T so that each unit triangle is subdivided into 
equilateral triangles of side 1/x; denote the new lattice by T x . The hypotheses of Conjecture 14.1 are 
clearly satisfied when the plurimers in A and B are regarded as plurimcrs on T x . Both their charges and 
their mutual distances get multiplied by x. By (2.13) one readily sees that the effects of the change in 
distance cancel out, while the effect of changing the charges is that the fraction —1/2 in front of the square 
brackets in (2.13) changes to —x 2 /2. 

A second "calibration" parameter a, a given positive integer, can be introduced as follows. Arrange, 
by a-fold lattice refinement, for all q^s in (2.13) to be divisible by a. Factor out the common factor of a? 
in front of the square parenthesis in (2.13) and re-denote each leftover qi/a by The formula obtained 
this way from (2.13) still parallels (2.14), but now the elementary physical charge q e corresponds to a 
plurimcr of charge a. In particular, a should have a fixed value. 

The overall effect of rr-fold lattice refinement and "calibration" by a is that the fraction —1/2 in front 
of the square brackets in (2.13) changes to —(ax) 2 /2. 
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paper) seems to be a possible such "machinery" that would produce, in a two dimensional 
world, precisely the effects of electrostatics. As presented in Section 14, we conjecture 
that the validity of our model does not depend on the particular choice of the hexagonal 
lattice, but it holds in fact for any plane bipartite lattice under a suitable embedding. 
This would imply that the electrostatic effects we obtain depend only on the space — in 
accordance with Feynman's suggestion [11, §12-7] that it might be the space itself, the 
common "framework into which physics is put," that is responsible for the emergence of 
such a simple equation governing electrostatics (as well as other physical phenomena). 

The model presented in this paper has natural analogs in higher dimensions — for in- 
stance, by considering the cubic lattices. We believe that the three dimensional analog 
presents the effects of electrostatics in the real three dimensional world. Some details 
about these considerations are presented in Section 15. 



3. Reduction to boundary-influenced correlations 



We will find it convenient to express the number of dimcr coverings of the region 
Hn ( ^ 1 • • • ^ m ; ^} ••• ^7 ] in terms of the number of tilings of two "halves" of it. 

One of them is the region Wn ( ^ ■ ■ ■ ^ m ; ^} ••• J defined in Section 2 — the por- 

V V l V m v l v n J 

tion to the west of the cut illustrated in Figure 2.3. The other is the part of Hn to the east 
of that cut — more precisely, what is left from the region east of the cut after removing the 
forced dimers (see Figure 3.1 for an example; the forced dimers are shaded), and weighting 
by 1/2 the N + 2m — 1 dimer positions below the monomer u that are closest to the cut; 

(Ft Ft R' Ft* \ 
1 • • • m ; / ■ • • ™ I ■ Figure 3.1 illustrates an example (for clarity, 

the two regions are pictured after separating them horizontally by one unit; the dimcr 
positions weighted by 1/2 are indicated by shaded ellipses). 

The regions Wn and En were defined in such a way that they are precisely the ones 
that result from applying the Factorization Theorem [3, Theorem 1.2] to the region Hn- 
Therefore, the quoted Factorization Theorem yields the following result. 



Proposition 3.1. 
M ( H N 



Rl Rm . R'l R' n | | 

vi v m ' v[ v' 



2 N+m+n M f w (Rl ... Rm R'l ... KX\ M ( £ ( Rl ... Rm R[ . . . R'n 

V \ v l v m V [ V' n J J \ \V 1 V m V[ V' n 



(3-1) 



In particular, (3.1) holds for the reference position of the plurimers used in definitions 
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(2.2) and (2.8). We obtain 



M H 



N 



2 N+m+n M[W N 



2m -1 1 

' ' 

1 2m 
"' 



2n-l 




1 







2n- 1 




MLB 



2m 



1 1 

' 



2n- 1 




(3.2) 



Dividing (3.1) and (3.2) side by side and letting N — > oo, one sees by (2.2) and (2.8) that 
lo and LOb can be related by considering another boundary-influenced correlation, defined 
by means of the regions Em'- 



Ri R m R[ R' n 



M E 



:= lim 

N— too 



Ri 

Vl 



Rm , R\ 



R 'n 



M E N 



2m 



2n-l 




(3.3) 



(the fact that this limit exists follows by Lemma 13.1). 

By (3.1)-(3.3), the definitions (2.2) and (2.8), and the existence of the limits in the latter 
two definitions (guaranteed by Lemma 5.1 and Lemma 13.1), one obtains the following 
result. 

u 



Proposition 3.2. The limit (2.2) defining the correlation at the center exists and its 
value is the product of the two boundary-influenced correlations: 



Ri 

Vl 



R-n 
Vrr 



R'l R'n 



Ri 

Vl 



Rm , Rl 

v m ' v[ 



R'n 



Ri 

Vl 



Rm . R'l 

V m ' Vi 



R'n 

,1 



(3.4) 



Therefore, in order to obtain the asymptotic behavior of the correlation at the center 
u>, it is enough to study the boundary-influenced correlations u>b and u>b- 

4. A simple product formula for correlations along the boundary 

Our calculations are built upon an explicit product formula that we present in this 
section for certain correlations along the boundary of the regions Wn- To state this we 
need to introduce a new family of regions, closely related to the Wjv's, and in terms of 
which the boundary-influenced correlation uob turns out to be expressible. 

Let W be the region determined by the common outside boundary of the regions 

W N ( Rl ■ ■ ■ Rm : R '} ■ ■ ■ R P ) , for fixed N, m and n. Then W is the half-hexagonal 
V vi v m v[ v' n ) 

lattice region with four straight sides — the southern side of length N + 2n, southwestern 

of length 2N + 4m, northwestern of length 2N + 4n + 1, and northern of length N + 2m — 

followed by N+2n descending zig-zags to the lattice point O, one extra unit step southwest 

of O, and TV + 2m more descending zig-zags to close up the boundary (the boundary of 

the region W corresponding to m = 4, n = 6, and N ~ 2 can be seen in Figure 4.1). In 

addition, the N + 2n dimer positions weighted 1/2 in the regions Wn are also weighted 

so in W. 

The eastern side of W can be viewed as consisting of bumps — pairs of adjacent lattice 
segments forming an angle that opens to the west: N + 2m bumps below O, and N + 2n 
above O. Label the former by 0, 1, ... , iV+2m— 1 and the latter by 0, 1, ... , N+2n—l, both 
labclings starting with the bumps closest to O and then moving successively outwards. 

Since everywhere in the description of W the parameters m and n appear with even 
multiplicative coefficients, we re-denote, for notational simplicity, 2m by m and 2n by 
n. Therefore we consider the four straight sides of W to have lengths N + n, 2N + 2m, 
2N + 2n + 1 and N + m, while the number of bumps below and above O is N + m and 
N + n, respectively. The results of this section do not assume that m and n are even (even 
though we will only use them for even m and n). 

We allow any bump above O to be "removed" by placing an up-pointing quadromer 
(lattice triangle of side two) across it and discarding the three monomers of W it covers. 
Similarly, a bump below O can be removed by placing a down-pointing quadromer across 
it and discarding the three monomers it covers. 

We are now ready to introduce the family of regions mentioned in the first paragraph 
of this section: define Wn\^i, ■ ■ ■ , k m ] h, ■ ■ ■ Jn] to be the region obtained from W by 
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Figure 4.1. W 2 [1,2, 3, 5; 0,1, 3, 4, 5, 6] 



removing the bumps below O with labels < k\ < k 2 < ■ ■ ■ < k m < N + m — 1, and 
the bumps above O with labels 0<Zi<Z2<---<^n<-^ + ^ — 1- Figure 4.1 shows 
W 2 [l,2,3,5;0,l,3,4,5,6]. 

The product formula referred to in the title of this section is stated in the following 
result. 

PROPOSITION 4.1. For m, n > and fixed integers < k\ < k 2 < ■ ■ ■ < k m and < h < 
l 2 < ■ ■ ■ < l n we have 



M(W N [h,..., 
N™oo M (W N [0, . . . , m - 1; 0, . . . , n - 1]) ~~ 

m (3/2), n (3/2), n n ^ 

TT l^Jfcj TT j-j-jjj l<a<j<m 1<»<J<" / , -, n 

i=li=l 

n n 

— v- TT ( 3 / 2 )fc. TT ( 3 / 2 )'. i<»<3<" i<»<J<" a9 N 

i=l i=l 

16 



where (a)k ■= a(a + 1) • • • (a + k — 1) (by convention, (a) := 1) 



W-ll ( i )i(3 /2) i 11 (3/2), • 



Proof. An explicit product formula for M (WV[fci, . . . , k m ; h, . . . , l n ]) follows from the 
results of [5]. Indeed, our regions ITjv[fci, . . . , k m ; h, . . . , l n ] are just a different notation 
for the regions -Rk,i(aO defined in [5, Section 2]: In the notation of [5], x is the length of 
the base — N + n in the current notation — while k and 1 are the (strictly increasing) lists 
of labels — incremented by 1, due to a shift of one unit in labeling bumps in the present 
paper compared to [5] — of the bumps that remain in the region. Therefore, 

W N [kx, . . .,k m ;h, ...,l n ] = 

R[i,...,N+ m ]\[k 1 +i,...,k m +i],[i,...,N+n]\[i 1 +i,...,i n +i](N + n). (4.4) 

Relations [5, (6.1)], [5, (6.2)] and [5, (6.4)] state that for any pair of lists p = [pi, . . . ,p s ], 
1 < Pi < P2 < ■ ■ ■ < Ps and q = [q lt . . . , q t ], 1 < q x < q 2 < . . . < q t , and for any 
nonnegative integer x>q t — p s — t + s (see [5, p. 10]), one has 

M (i? P , q (x)) - c p , q F p , q (x), (4.5) 

where 

- .,; ■ If 1 TT 1 ^<i<i<s^i - P^i<i<i<t^ '"ft) U (\\ 

p ' q 11 (2P, - 1)1 ll ( 2ft )i nun5=i(Pi+*) 

and the polynomials F Piq (x) satisfy 



fp^qW (g) 



= (x-pi +p s )(x+pi +p s - s + t+ 1), for 1 < i < s (4.7) 
= (x + qi+p s + l)(x- qi+p s - s + t), for 1 < i < t (4.8) 



(here p'^ is the list obtained from p by increasing its i-th element by 1 — in particular, it 
is defined only if l i+ i — h> 2). 

To deduce the limit (4.1) from these formulas, it will be convenient to consider first the 
limit 

M(W N [k u . . . ,ki- U ki + l,k i+1 , . . . ,k 
n^oo M(W N [ki,...,ki-i,ki,ki +1 ,---,k m ;h,...,l n ]) 
for ki + 1 < fcj+i. 

Use (4.4) to view the regions involved in this fraction as -R p , q (x)'s. The lists p and q 
corresponding to the regions at the numerator and denominator in (4.9) are illustrated in 
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fc 2 +l 


i f +2 




k m +l 






1 • 


• k x *2- 


• -k 2 k 2 +2 ■ ■ 


• kj+l kt+3 
kf+2 


' ' -k M k M *2 • 


■ ■ k k - 


a- ■ -N+m 






k?\ 


k 2 +l 












Z 1+ l 




(a). 


The p's. 




Z„+l 




1- ■ 


■h 'i + 2 ■ 


■ ■ h '2+2 


(b). 


The q's. 




/„ l„ + 2 



Figure 4.2. Incremental change in the first index list. 

Figure 4.2 (the shaded squares indicate the entries removed at the indices on the right 
hand side of (4.4)). Since the q- lists are the same, the limit (4.9) can be found applying 
formulas (4.6) and (4.7) to the lists illustrated in Figure 4.2. 

Since the only difference between the p-lists in Figure 4.2(a) is the one indicated by 
the dot in that figure, we obtain from (4.6) that the contribution to the fraction in (4.9) 
coming from the ratio of the c Pjq 's is 



fcj [N + m - {h + 



(2h + 1)! (kj - fci) ■ ■ ■ (kj - kj-i) (kj+i - kj) ■ ■ ■ (k m - h) 

1 (h + iy. [N + m-(k l + 2)}\ 

VK+W (fci - fei + 1) ■ ■ • (h - fc 4 _x + 1) {k l+1 - h - 1) • • • (k m -h-i) 

{h + 3) ■ ■ ■ {kj + 2 + N + to) 
(fci + h + 3) ■ ■ ■ (kj + L + 3) 
X (h + 2) ■ ■ ■ (kj + 1 + N + to) 
(ki + h + 2) • • • {ki + l n + 2) 

{{ki + 1) - fci) • • • {{ki + 1) - fci_i)(fci+i - {ki + !))••• (fc m - (fc, + 1)) 



{ki fci) • • • {ki fci — 1 ) (fci+ 1 fcj) • • • {k m fci) 
(fci + li + 2) ■ ■ ■ {kj + In + 2) 

(fci + fi + 3) ■ ■ ■ (fci + Z„ + 3) ' 



(4.10) 



On the other hand, the contribution to the fraction in (4.9) coming from the polynomials 
Fp,q{%) is readily seen from (4.7) to be 

; (4.11) 



{2N + m + n- {k t + 1))(2N + to + n + (fc, + 1) + 1) 
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By (4.10) and (4.11) one readily finds the limit (4.9) to be 



M {Wj^ku + 1, k i+1 , . . . , fc m ; l 1; . . . , l n }) = 

n-,oo M (W N [ki, ki-i,ki, fcj+i, . . . , k m ;h, . . . , Z n ]) 

fc, + 3/2 ((fc, + 1) - fci) ■ ■ ■ ((fcj + 1) - fcj-iXfcj+i - (fcj + !))■■■ (fc m - (fcj + 1)) 
k{ H~ 2 (fe^ — &i) • • • (/^ ) (A^_|_i — foj) • • • (/c m fcj) 

(fci + /l + 2) • • • (fej + l n + 2) 

((fci + 1) + ii + 2) ■ ■ ■ ((A* + 1) + i„ + 2) ' 1 ' J 



Regard (4.12) as giving the effect of decrementing the entry fcj + 1 in the argument to 
fcj. A further decrementation of this entry to ki — 1 will produce, by (4.12), the factor 



+ 1/2 



(fci - fcl) • • • (fci - fcj_i)(fc; + i 



(fcm fci)) 



fc, + l ((fc; -1)- fci)- ••((&; -1)- 
(kj + h + !)■■■ + /» + 1) 
X (ki + h+2)--- (fci + Z„ + 2) ' 



fci-l)(fci+l - (fci - I))-" (fcm - (fci - 1)) 



(4.13) 



One notices that telescoping simplifications occur when multiplying together the second 
and third fractions in (4.12) for successive decrementations of each argument fci. In partic- 
ular, if the argument fci is decremented all the way to 0, we obtain by repeated application 
of (4.12) that 



lim 



M(W N [k u k 2 



■ j ^1 ; • • • ; t>n\ ) 



M(W N [0,k 2 ,...,k m ;h,...,ln]) 
(3/2) fcl (fc 2 - fci)(fc 3 - fci) • • • (fcm - fci) (h + 2) • • • (l n + 2) 



(2) fel (fc 2 -0)(fc 3 -0)- 



(fcm 
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0) (fci+Zi + 2)---(fci+Z n + 2)' 



(4.14) 



Similarly, we obtain 



M(W N [0,k 2 ,k 3 ,...,k 
n^oo M(W N [0,l,k 3 , . . . ,k m ;h, . . . ,l n ]) 

(5/2) fc2 _i (fc 2 - 0)(fc 3 - k 2 ) ■ ■ ■ (fc m - fc 2 ) (Zi + 3) ■ ■ ■ (/ n + 3) 

(3) fe2 _! (l-0)(fc 3 -l)"-(fcm-l) (k 2 + l 1 + 2)---(k 2 + l n + 2)' 
M (W N [0, 1, fcgj fc 4 , ■ ■ ■ , k m ;h, l n }) = 

jv^oo M (VKjv[0, 1, 2, fc 4 , . . . , fc m ; h, . . . , /„]) 

(7/2) fc3 _ 2 (fc 3 - 0)(fc 3 - l)(fc 4 - fc 3 ) ■ ■ ■ (k m - k 3 ) (h + 4) ■ ■ ■ (Z» + 4) 

(4) fe3 _ 2 (2-0)(2-l)(fc 4 -2)---(fc m -2) (fc 3 + ; 1 + 2)---(fc 3 + ;„ + 2)' 



M(Wjy[0,l,...,m-2,A: m ;ii,...,f n ]) = 
n"^ M(Wjv[0,l,...,m- l;Zi,...,Z n ]) 

((2m + l)/2) fcm _ m+ i (fc m - 0)(fc m -!)■■■ (fc m - (m - 2)) 

(m + l)fe m -m+i ((m - 1) - 0)((m - 1) - 1) • • • ((m - 1) - (m - 2)) 
(Zi + ro + 1) ■ ■ ■ (In + m + 1) 
X (fcm + ii + 2) ■ ■ ■ (fc m + Z„ + 2) ' 



(4.15) 



Now multiply together the equalities in (4.14) and (4.15). The first fractions on their 
right hand sides combine to give I12=i((3/2)fc i /(2)fc ( )/n2=i((3/2)i-i/(2) i _i). Due to 
simplifications, the second fractions on the right hand sides yield Y\i <i< j <m (kj—ki)/(j—i). 
The third fractions give IlILift + 2 W II™ i lXj=i( k i + h + 2 )- We obtain 

M(W N [k 1 ,k 2 ,...,k m ;h,...,l n }) _-fj n , ^ 
N^oM(W N [0A,..., m -l;h,...J n })-l\ (k + 2)m 

yi™ (3/2) fci kj - kj 

(2)ki lll<i<j<m j-i 

n"i nri n"-i(fci + ^- + 2)' (4 . 16) 

The effect of decrementing arguments in the list [Zi, . . . , l n ] can be analyzed in a similar 
way. Indeed, compare the situation when the lists of kj's at the indices on the left hand 
side of (4.1) are the same, and the lists of Ij's are [Zi, . . . , h-i, h + 1, h+i, ■ ■ ■ , l n ] an d 
[Zi, . . . , Zj_i, li, Zj + i, . . . , l n ], respectively. These lists are illustrated in Figure 4.3. We need 
to find 

lim M(W N [ki, ... ,fc m ;Zi, ...... . ,k-i,k + l,Z m ,Z n ]) ^ 

n^oo M(W N [ki, . . . ,k m ;h, . . . ,k-i,li,l i+ i, . . . ,l n ]) 
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■ k 2 k 2 +2 


' ' ' k m k m +2 ' ' ' N+m 










(a). The p's. 










/;+2 /,-+!+! /„+l 






1 • 


■ ■ h '1+2 ■ • 


• h l 2 *2- ■ 


i; + 3 ■ ■ ■ /, +1 +2 • • ■/„ l m+ 2 


■ ■ N+n 








h+2 




• 






; 2 +i 


h + l + l l n *\ 





(b). The q's. 

Figure 4.3. Incremental change in the second index list, 
for h + 1 < l t+1 . 

By (4.4), view the regions in (4.17) as i? p . q (x)'s and use formulas (4.5), (4.6) and (4.8). 
Repeating the reasoning that proved (4.12), we see that the formulas we obtain now are 
almost exactly (4.10) and (4.11), with the roles of the lists [k\, . . . ,k m ] and [Zi,...,Z„] 
interchanged. 

Indeed, the only difference from (4.10) of the contribution coming from the c Piq 's of 
(4.6) is that the first fraction on the left hand side of the analog of (4.10) is (l/(2(Zj + 
l))!)/(l/(2(Zi + 2))!), as opposed to the (l/(2l i + l)\)/(l/(2l i + 3)\) that results from (4.10) 
by interchanging the lists [k\, . . . , k m ] and . . . , l n }. The consequence of this difference 
is that the numerator of the first fraction after the equality sign in the present situation 
analog of (4.10) is (2Zj + 3)(2Zj + 4), and therefore the first fraction on the right hand side 
of the analog of (4.12) is (h + 3/2)/(Z; + 1). 

On the other hand, the contribution to the left hand side of (4.17) coming from the 
fpq's of (4.6) is readily seen to have, after letting N — > 00, exactly the same effect as 
(4.11). 

Therefore formulas (4.14) and (4.15) have perfect analogs when changing the arguments 
of the second list, with the only difference that the integers between the round parentheses 
at the denominators in the first fractions on their right hand sides are now decremented 
by one unit: they are (l) (l , (2)j 2 _i,..., (n); n _„+i. 

Just as we deduced (4.16) from (4.14) and (4.15), we obtain from the above analysis of 
the differences between decrementing the fci's and the Zj's that 



,. M {W N [k u ...,k m ; h,l 2 , l n ]) _ 
n^oo M (W N [h, . . . , k m ; 0, 1, . . . , n - 1]) " 11^ + Z)n 



n" (3/2);, „ Ij-U 
i=1 (l); t lll<J<i<n j-i 

mu nri n^fe+^ + 2)' (418) 



(1). 

Combining (4.16) with the specialization of (4.18) for ki = i — 1, i = 1, ...,m, and 

1;= 



using — 1) + 2) n = JTiLi n?=i(* + i)i we obtain the first equality in (4.1). The 
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second equality is just a repackaging of the first, using Yli<i<j< m U ~ *) = an d 

n£in?=i(<+i) = n"=o(j+2)n. □ 

5. A (2m + 2n)-fold sum for uii, 

In this section we present an expression for the boundary-influenced correlation ujf, as 
a (2m + 2n)-fold sum. We deduce this by expressing the generating function for dimer 

(R R R f R f \ 
1 • • • m ; } ■ ■ ■ " I in terms of generating functions for 
vi v m v 1 v n J 

the dimer coverings of the regions Wjv[fci, . . . , k m ; h, . . . , l n ] introduced in the previous 

section, and then using the explicit product formula (4.1). The remaining part of the 

paper will consist mainly of analyzing the asymptotics of this multiple sum. 

The result mentioned in the title of this section is the following. 

Lemma 5.1. For fixed Ri, . . . ,R m ,R[, . . . , R' n > 1 and Vi, . . . ,v m ,v[, . . . ,v' n > we have 



K 

V" (_l)E£i(«*+6i)+Er=i( c *+*) 



g (Ri + <n - 1)1 (Ri + bt- 1)1 (3/2)„ (+0i (3/2)„ (+6i 



Ri 




.R'i 


v m 




Ri 

E •• 

oi,6i=0 


R m 

E 


E 

c 1 ,d 1 = 



1 (2a,)! (Ri - a,)! (26,)! (Ri - h)\ (2) Vi+ai (2) Vi+bi 

A (R l + c l -l)\(R l + d l -l)\ (3/2)„, +Ci (3/2)„,+ d< 

X 11 (2c, + 1)! (Ri - ci)! (2d, + 1)! (ifc - rf,)! (1) <+Ci (l) <+di 

x j [ (vj -Vi + a,j - a,i)(vj -Vi + a,j - bi)(vj -Vi + bj - di)(vj -Vi + bj - h) 

l<i<j<m 

X Yl ( V 'j - V 'i + C 3 ~ C i)( V 'j ~ V 'i+ C 3 - d *)( V 3 - V 'i + d 3 - C i)( V 'j ~ V 'i + d 3 ~ *) 
l<i<j<n 

n?=M-bi) 2 n? =1 (ci-di) 2 



II™ 1 ]l?=l( W ii + a i + C j)( U ij + a i + d 3)( u i3 + b i + c j)( u i3 + b i + d j) 



l 3 

where 



(5.1) 



Uij =Vi + v'j+2 



for i = 1, . . . , m, j = 1, . . . , n, and Xm,n is given by (4.3). 

To prove this Lemma we will need the following special case of the Lindstrom-Gessel- 
Viennot theorem on non-intersecting lattice paths (see e.g. [13] or [29]). 

Consider lattice paths on the directed grid graph Z 2 , with edges oriented so that they 
point in the positive coordinate directions. We allow the edges of 1? to be weighted, and 
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define the weight of a lattice path to be the product of the weights on its steps. The weight 
of an TV-tuple of lattice paths is the product of the individual weights of its members. The 
weighted count of a set of iV-tuples of lattice paths is the sum of the weights of its elements. 

Let u = (ui, . . . , ujv) and v = (vi, . . . , Vn) be two fixed sets of starting and ending 
points on Z 2 , and let W(u, v) be the set of non-intersecting TV-tuples of lattice paths with 
these starting and ending points. For P £ Af(u, v), let erp be the permutation induced by 
P on the set consisting of the N indices of its starting and ending points. 

Theorem 5.2 (Lindstrom-Gessel-Viennot). 

J2 (-1) CTP wt(P) = det ((a 4 ,)i< 4 ,,<„) , 
PeAf(u,v) 

where a,ij is the weighted count of the lattice paths from Ui to vj . 

What makes possible the use of this result in our setting is a well-known procedure 
of encoding dimer coverings by families of non-intersecting "paths of dimers:" given a 
dimer covering T of a region R on the triangular lattice and a lattice line direction d, the 
dimers of T parallel to d (i.e., having two sides parallel to d) can naturally be grouped into 
non-intersecting paths joining the lattice segments on the boundary of R that are parallel 
to d, and conversely this family of paths determines the dimer covering (see Figure 5.1 for 
an illustration of this, and e.g. [5] for a more detailed account). 

Consider such a path of dimers P. Let T denote our triangular lattice. Clearly, P can 
be identified with a lattice path on the lattice C of rhombi formed by the midpoints of 
the segments of T that are parallel to the encoding direction d. In turn, by deforming the 
lattice of rhombi £ to a square lattice, the path of dimers P can be regarded as a lattice 
path on Z 2 . It is in this sense that we will view the paths of dimers as lattice paths on Z 2 
in the remainder of this section. 

Proof of Lemma 5.1. Choose the lattice line direction d in the above encoding pro- 
cedure to be the southwest-northeast direction of T. Encode the dimer coverings of 

W N ( Rl ■ ■ ■ ^ m ; R } ■ ■ ■ R P ) by (27V + 2m + An + l)-tuples of non-intersecting paths— 

2N + An + 1 starting at the northwestern side, and 2m starting at the m down-pointing 
removed quadromers — consisting of dimers parallel to d (see Figure 5.1; there and in the 
following figures the N + 2n dimer positions weighted by 1/2 are not distinguished, but 
are understood to carry that weight). 

Let P be such a (2N + 2m + An + l)-tuple. Consider the permutation op induced by 
P on the set of the 2N + 2m + An + 1 indices of its starting and ending points. We claim 
that the sign of ctp is independent of P. Indeed, denote by Di, i = 1, . . . ,m and Uj, 
j = 1, . . . , n, the removed down-pointing and up-pointing quadromers, respectively. While 
the way in which the starting and ending points of P — clearly independent of P — are 
matched up depends on P, it is always the case that the two paths starting at Di end at 
consecutive ending points, and the two ending on Ui start at consecutive starting points. 
It is easy to see that this implies that all op 's have the same sign. 
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5 2 4 2 3 

a tiling of W2 I ; „ ) • sion over the rows indexed by a and /?. 



Figure 5.1. Lattice path encoding of Figure 5.2. The effect of Laplace expan- 
'5 2.4 2 3" 
1' 1 2 4 / 

Consider the lattice £ defined just before the beginning of this proof. Weight by 1/2 its 

Ri R m R 1 R' n 



edges corresponding to dimer positions weighted by 1/2 in Wn 

\ v i v m v'i v'„ 
Weight all other edges of £ by 1. Then the weight of the dimer covering encoded by P is 
just wt(P), and we obtain by Theorem 5.2 and the constancy of the sign of crp that 

M (w N ( Rl ■■■ Rm ; R '} ■ ■■ R '?\\ = Idet A\ , (5.2) 
\ \vi v m v[ v' n JJ 1 

where A is the (27V + 2m + An + 1) x (2N + 2m + An + 1) matrix recording the weighted 
counts of the lattice paths on C with given starting and ending points (note that the right 
hand side of (5.2) is independent of the ordering of these starting and ending points). 

We deduce (5.1) by applying a sequence of Laplace expansions to the determinant in 
(5.2). Recall that for any m x m matrix M and any s-subset S of [m] := {1, . . . ,m}, 
Laplace expansion along the rows with indices in S states that 

det M = Y,(- 1 Y (K) dct M s det MgJ* (5.3) 

K 
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where K ranges over all s-subsets of [m], e(K) :— ^ keK {k — 1) and Mj is the submatrix 
of M with row-index set J and column-index set J. 

The rows and columns of the matrix A in (5.2) are indexed by the starting and ending 
points of the (2JV + 2m+4n + l)-tuples of non-intersecting lattice paths encoding the dimcr 

(R R R r R f \ 
1 • • • m ; / ■ ■ ■ " ) . The starting points are the 2N + An + 1 unit 
V\ v m v 1 v n J 

( R R R! R! \ 
segments along the northwestern boundary of Wn I 1 • • • m ; / ■ ■ ■ " J , together 



with the segments a* and /3, of D t that are parallel to d, i = 1, . . . , m (see Figure 5.1). 
The ending points are the 2iV + 2m + 2n + 1 segments parallel to d on the eastern boundary 

(R R R f R 1 \ 
1 • • • m ; / • • • " ] , together with two more such segments on each Uj , j = 
Vl v m v 1 v n J 

1, . . . ,n. 

Label the bottommost N + 2m ending segments on the eastern boundary, from top to 
bottom, by 0, 1, . . . , N + 2m — 1 (this labeling is illustrated in Figure 5.1). 

Apply Laplace expansion to the matrix A of (5.2) along the two rows indexed by a\ 
and 0i (see Figure 5.1). The first determinant in the summand in (5.3) is then just a 
two by two determinant. Its entries are weighted counts of lattice paths on C that start 
at «i or 0i and end at some labeled segment on the eastern boundary. There are only 
Ri + 1 labeled segments that can be reached this way: those with labels ranging from 
vi to vi + R\. We can restrict summation in (5.3) to the two-element subsets K of this 
set of segments, since all other terms have at least one zero column in the two by two 
determinant. Therefore we obtain from (5.2) that 



M W N 



Rl Rm . R[ R' n 
Vl 



E/_ i\ai+6i i , .{m+«i,tii+()i} i , 4 [2JV+2m+4n+l]\{vi+ai,i;i 
^ l > aCtyi {aiA} aCT/1 [2Ar+2m+4n+l]\{ Ql ,/3 1 } 



+bl} 



0<ai<bi<7?i 



(5.4) 



(here [2N + 2m + 4n + 1] denotes the set of starting, respectively ending points of the 
families of non-intersecting lattice paths encoding tilings of the region (5.2)). Choosing 
the origin of the lattice C to be at ai, and its positive axis directions to point east and 
southeast, one sees that 0i has coordinates (—1,1), and the segment labeled Vi +j on the 
eastern boundary has coordinates (i?i — 1 — j, 2j), j = 0, . . . , Ri. Since the lattice paths 
counted by the entries of A^^p^ Vl+hl ^ have all steps weighted by 1, the determinant of 
this matrix is 



dfitA{" 1+ "Y Bl+6l} =det 

{ai,Pi} 



r fiJi-l+aA (R 1 -l+bi\ i 
\ 2ai )\ 26i 
/ili-l+oi\ /JJi-l+6i 
V 2ai-l )\ 26i-l 



(bi-ai)(Ri+ai-l)\(Ri +b 1 -l)\ 

1 (2ai)\(Ri-ai)\(2bi)\(Ri-bi)\ ' [ ' } 
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On the other hand, the second determinant in the summand in (5.4) can be interpreted 
as being the weighted count of dimer coverings of the region 



W N 



(R 2 

V v 2 



R m R[ 



R'n 

vL 



Vi + Cbi 

vi + bi 



(5.6) 



obtained from W, 



N 



Rn 



R[ 



R'n 



by placing back quadromer D\ and remov- 



Ri 

v\ v m ■ v\ v„ 

ing the two monomers that contain the segments labeled v\ + a\ and v\ + bi on the eastern 
boundary (see Figure 5.2 for an illustration). Indeed, the Lindstrom-Gessel-Viennot ma- 
trix of this region is precisely ^[27v+2m+4n+i]\{ai"^i 1 } 1 ' 1+bl ^' anc ^ tne ar g umcnt that 



proved (5.2) we obtain that M I Wn 



detA 



[2AT+2m+4n+l]\{j;i+ai,i;i+6i} 



R'n 



Vl 



Oi 



is equal to 



R2 Rm . R[ 

V 2 V m ' V[ 

, , . .. n , , : ; . . r , up to a sign that is independent of ai and b\ (to check 
the latter statement, in the labeling of the starting and ending points of P that defines 
crp , label the bottommost TV + 2m — 2 northwest facing unit segments on the boundary of 
(5.6), say from bottom to top, by 1, . . . , N + 2m — 2, irrespective of the values of a\ and 
61 ; the argument in the second paragraph of the current proof shows that the sign of ap 
is independent of not just P, but also of a\ and 61). Therefore, using (5.5) we can rewrite 
(5.4) as 



M W N 



Ri 

Vl 



2R 1 





R[ R' n 


V m 


v'i v' n 


E 


(— l)° 1+bl 



61-1)! 



0<ai<6i<fli 

X M I W N 



R2 

V2 



(2a 1 )\(R 1 -a 1 y.(2b 1 y.(R 1 -b 1 )\ 



Rn 

V„. 



R'n 

n,' 



v\ - 

Vl 



01 
■61 



(5.7) 



This way the matching generating function of M I Wj 



'N 



.R'i 

■A 



R' 



Vl V m ' V\ V n 

pressed as a double sum involving the matching generating functions of regions having 
one less down-pointing triangular hole. By the same reasoning the second down-pointing 
quadromer D 2 can be removed from the regions (5.6) and the last factor in the summand 
of (5.7) is expressed as a double sum involving regions with two less down-pointing holes. 
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One obtains 



M W N 



4RiR 2 



Ri 

Vl 



Rm . R'l 



0<ai<6i<fii 0<a 2 <b 2 <R 2 



(-l) ai+bl+a2+b2 sgn(ui +oi,«i + 61, v 2 + a 2 ,v 2 + b 2 ) 



(6i-ai)(i?i 



ai 



!)!(#! + &!-!)! 



(2oi)! (Jii - 0l )! (26x)! (fl! - 6x)! 

(62 - a 2 )(i?2 + a 2 - 1)! (^2 + frz - 1)! 
(2a 2 )!(i? 2 -a 2 )!(26 2 )!(i? 2 - 6 2 )! 



x M Wn 



(R 3 

\ u 3 



vL 



Vl 
Vl 



■ ai 



v 2 + a 2 
v 2 + b 2 



(5.8) 



where the summation extends only over those indices for which (t>i + 01, t>i +&i, t> 2 +a 2 , v 2 + 
6 2 ) consists of distinct components, sgn denotes its sign when regarded as a permutation, 
and the region on the right hand side of (5.8) is obtained from (5.6) by placing back the 
quadromer D 2 and removing the two monomers containing the segments labeled v 2 + a 2 
and v 2 + b 2 on its eastern boundary. 

Indeed, perform Laplace expansion for the region (5.6) over the rows indexed by a 2 
and (3 2 . Because of the geometry of the triangular lattice T, the paths starting at these 
segments can end on the eastern boundary only at segments with labels in the range 
[v 2 ,v 2 + R 2 \; denote these labels by v 2 + a 2 and v 2 + b 2 . Moreover, because the segments 
labeled vi + a\ and vi + bi are not present in the region (5.6), the labels v 2 + a 2 and 
v 2 + b 2 must in fact also be different from both v 1 + ai and vi + b\. This explains why the 
summation range in (5.8) needs to be restricted to vi + a, v\ + bi, v 2 + a 2 , V2 + b 2 being 
distinct. 

The only other needed explanation for justifying (5.8) is the factor sgn(i;i + ai,vi + 
bi,v 2 + a 2l v 2 + b 2 ). The reason it appears is the following. In the summand of the Laplace 
expansion (5.3) one has the factor (— 1)^ K \ In our situation this becomes (— l) fel+fe2 , 
where the "column indices" fcj indicate that the two path-ending segments on the eastern 
boundary are the kith and fc 2 th from, say, the top, in the current region. This column 
index coincides — up to a translation, which, pertaining to both ki and k 2 , does not affect 

the sign (-l) e W— in the case of M [W N ( Rl ■ ■ ■ Rm ; R '} ■ ■ ■ R 'p ) ) , with the label of 

\ \Vl V m V 1 v n J J 

the path-ending segment. However, when we do Laplace expansion for the region (5.6), 
this column index is affected — unless v 2 + a 2 < v 2 + b 2 < v\ + ai < vi + bi — by the 
previous removal of the segments with labels vi + ai and vi + bi . 

More precisely, say v 2 + a 2 < v\ + a\ < v 2 + b 2 < vi + bi. Then while v 2 + a 2 
correctly gives the column index of the segment where the path starting at a 2 ends, 
v 2 + b 2 is, due to the absence of the segment labeled vi + 01, one unit more than the 
column index of the ending segment of the path starting at (3 2 . It is easy to see that in 
general the effect of this "interference" is precisely the multiplication of the summand by 
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sgn(v 2 + a,2, V2 + b 2 ,v 1 +a 1 ,v 1 + bi). Since clearly sgn(v 2 + a 2 , v 2 + b 2 ,v 1 + ai,vi + b x ) = 
sgn(«i + a\, Vi + 61, v 2 + a 2 , v 2 + b 2 ), (5.8) is completely justified. 
Applying the same reasoning m — 2 more times, we obtain 



M W N 



R-n 
V„ 



R[ 



R' n 



Y (-l) ET = l(ai+bl) sgn(v 1 +a 1 ,v 1 +b 1 ,...,v m + a m ,v m + b rn ) 

0<oi<6i<_Ri a<a m <b m <R m 



yr (bj - aj^Ri + a % - 1)! (R l + h~ 1)! 
fj (2a i )\(R i -a i )\(2b i )\(R i -b i )\ 



x M Wjv 



1 . R 'l R 'n 



vi+ai v m + a 
Vi + 61 v m + b„ 



(5.9) 



where the region on the right hand side of (5.9) is obtained from the region on the left by 
placing back all m down-pointing quadromers, and removing instead the 2m monomers 
containing the segments with labels Vi+ai and Vi+b i} i = 1, . . . , m, on its eastern boundary. 

We complete expressing the left hand side of (5.1) in terms of regions with no holes by 
repeating our hole-removing procedure for the regions 



W, 



N 



R'l R'n 



Vl 



ai 



+ a r , 



(5.10) 



on the right hand side of (5.9). 

To this end, we encode the tilings of the regions (5.10) by paths of dimers parallel to the 
direction d now chosen to be the southeast-northwest direction (an example is illustrated 
in Figure 5.3). As in the previous encoding, weight by 1/2 those segments of the resulting 
lattice C that correspond to dimer positions weighted 1/2 in the region (5.10), and weight 
all remaining segments of C by 1 . 

Each tiling of the region (5.10) gets encoded this way by a (2N + Am + 2n)-tuple P 
of non-intersecting lattice paths on £, starting at the unit segments on its southwestern 
boundary or at the unit segments 7^ and 8i of Ui (i = 1, . . . , n) that are parallel to d, and 
ending at the unit segments parallel to d on its eastern boundary (see Figure 5.3). 

Label the topmost N + 2n ending segments on the eastern boundary, from bottom to 
top, by 0, 1, . . . , N + 2n — 1 (this labeling is illustrated in Figure 5.3). 

As in the argument that proved (5.2), the sign of the permutation up is independent 
of P. Therefore, we obtain by Theorem 5.2 that 



M W N 



0. R'l K 



vi + ai 
vi + h 



+ a r , 
. + b ri 



= €1 det B, 



(5.11) 



where B is the (2N + 4m + 2n) x (2N + 4m + 2n) matrix recording the weighted counts 
of the lattice paths with specified starting and ending points, and the sign ei in front of 
the determinant is the same for all choices of and 6j, i = 1, . . . , m. 
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Figure 5.3. Lattice path encoding of 



a tiling of W 2 I jjj ; * ^ 4 



2 
4 3 



FIGURE 5.4. The effect of Laplace expan- 
sion over the rows indexed by 7 and S. 



Apply Laplace expansion in det B along the two rows indexed by 71 and 8\ . The first 
determinant in the summand of (5.3) is again two by two, and records weighted counts 
of lattice paths starting at 71 or 8\ and ending at some segment on the eastern boundary 
(see Figure 5.3). There are + 1 segments on the eastern boundary that can be reached 
this way; in the labeling of the paragraph before (5.11), they are the ones with labels 



R[. As with our previous Laplace expansions, we can restrict the 



summation range in (5.3) to obtain 



M W N 



'. R 'i 



R 'n 

vL 



Vi + (li 

vi + bi 



v m + b, n 



(- 1 \a+di J„ + R {«i+ci,fi+di} , „[2N+4m+2n]\{v' 1 +c 1 ,v' 1 +d 1 } 
tl 2^ V ' {7i,5i} QCL£i [2Ar+4 ; 



4m+2n]\{ 7l ,5i} 



(5.12) 



Centering £ at 71 and choosing the positive directions in the lattice C to point east 
and northeast, 8\ has coordinates (—1, 1), and the segment labeled v[ + j on the eastern 
boundary has coordinates (R[ — 1 — j, 2j + 1), j = 0, 1, ... , R[. The weighted counts 
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involved in the entries of £?r^^i'" 1+dl ^ (which involve this time some steps of weight 
1/2) are easily calculated and one obtains 



dctB { { ;j+^ +dl} =det 



l(R{-l+ci\ , fiii-l+ciN l/iiJ-l+dA , ("Jii-l+dA 

2 V 2ci / ' V 2ci+l / 2 V 2d! / ' V 2di + l ) 

1 /fli-l+cA , /Ri-l+ci\ l/i?i-l+di\ , /«i-l+di\ 

2^ 2ci-l )'\ 2ci / 2 V 2di-l / ' V 2di / 



= 2i?' 1 (i?' 1 -l/2)(i?' 1 + l/2) 



(2 Cl + 1)! (R[ - ci)! (2^ + 1)! (fli - d!)! ' 



(5.13) 



On the other hand, by applying Theorem 5.2 one more time one sees that 



, R [2JV+4m+2n]\{t;i+ci,t)i+di} _ 
aCT £> [2Ar+4m+2n]\{ 7 i,5i} 



M W N 



where 



W N 



R'o 



i?2 

^2 



v\ + a\ 
vi + h 



■ Ci 



+di 



Oi 

6i 



- b„ 



ci 
di 



(5.14) 



(5.15) 



is the region obtained from the region (5.10) by placing back quadromcr U± and removing 
the two monomers containing the segments labeled v[ + c\ and v[ + d\ on the eastern 
boundary, and the sign e[ is independent of c\ and d\. Therefore, by (5.12)-(5.14) we 
obtain that 

' • ... R '"\ \ v i+ a i ,_v m + a„ 
vi + bi v m + b 



M W, 



N 



0' v[ 



= 2i?' 1 (i?' 1 -l/2)K + l/2) 



X t\€ 1 



V , lYl+dl (di-cQ^+d-^i^+d!-!) 
^ " (2 Cl + 1)! (R[ - ci)! (2d! + 1)! (R[ - di)! 



0<ci<di<ili 

x M f Wat 



i? 2 

«2 



«1 



Oi 
■&1 



1>„ 



+ «r, 



■ Ci 

dj 



(5.16) 

By further Laplace expansions applied on the pairs of rows indexed by 7, and Si, i = 
2, . . . , n, one can successively remove all remaining holes U2, ■ ■ ■ , U n . The argument that 
proved (5.9) yields now 

R[ R' n \ Vi + Oi _ _ _ v m + a n 
v[ v' n J vi+ bi v m + b„ 



M W N 



x £ 

0<ci<di<_Ri 



v (-i)S?=i(ci+di) sgn (^ + Cl ,v' 1 + di 

0<c„<d„<H' 



n^n2^(^-i/2)(^+i/2) 

X + C n ,D^ + d n ) 



i=l 



n 



(d z -d)(^ + d- l)!(i^ + dj-l)! 
(2d + 1)! {R[ - d)! (2d 4 + 1)! {R[ - *)! 

v A/r f w ^ \ [«i + ai u m + fl m «i 
' V[| WjV (,0'0j [v,+b 1 -- 
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ci 



w m + b 

in ' wi + di 



< + d„ 



(5.17) 



Figure 5.5. The regions W N 



2 

3' 3 



and Wjv[0,2,3,4;0,l,2,3,5,6] 



have the same matching generating function. 



where the multiple sum extends only to the summation indices for which v[ + ci , v[ + 
di, . . . , v' n + Cn, v' n + d n are distinct, the signs a and e-, i = 1, . . . , n, are independent of 
oi, h, . . . , a m , b m , and 



W N 



vi + ai 
vi + h 



+ a r , 

, + br, 



v[ + a 
v[+d! 



v' n + c n 
v' n + d n 



(5.18) 



is the region obtained from the region (5.10) by placing back all its removed quadromers 
Ui, i = 1, . . . , n, and removing the 2n monomers containing the segments labeled 7, and 
Si, i = 1, . . . , n, on its eastern boundary. 

However, the region (5.18) differs from the region WAr[{«i+ai, v\+b\, . . . , v m +a m , v m + 
b m }<', {v'i + Ci, v[ + d\, . . . , v' n + c n , v' n + d„}<] (where A < denotes the list of increasingly 
sorted elements of the set A of distinct nonnegative integers) defined in Section 4 only in 
that the former contains 2m + 2n more unit rhombi near the eastern boundary (sec Figure 
5.5; the additional rhombi are shaded). Moreover, all these additional dimer positions have 
weight 1 and are forced to be part of any dimer covering of the former region. Therefore 
the two regions have the same matching generating function, and replacing (5.17) in the 

31 



right hand side of (5.9) one obtains 



M (w N ( R v * • • • ^ ; ^ • • • ^ ) ) = 2 m +" ^ f[ ^ ( R i V2) + 1/2) 

... v ^ ••• V (-i)Sr=i(«i+i'i)+Er=ife+*) 

0<ai<6i<i?i 0<o m <6 m <iJ m CKc^difC-R'j 0<c n <d n <R' n 

X f = \ (2a i )\(R i -a i )\(2b i )\(R i -b i )\ 

tt (^-c a )(^ + c a -l)!(^ + d t -l)! 
11 (2c 4 + 1)! (^ - a)\ (2d t + 1)! - di)\ 

x sgn(ui + ci, ^ + di, . . . , v' n + c„, v' n + d n ) sgn(wi + a u v 1 + &_, . . . , v m + a m , v m + b m ) 
x M(W N [{ Vl +a u ...,v m + & m }<;K + ci, + d n }<])| , (5.19) 

the summation extending only over those summation indices for which both lists of argu- 
ments in the region on the right hand side have distinct elements. 

After removing the forced dimers, the normalizing region at the denominator of (2.8), 

fl 3 2m -1.1 3 2n-l\ , . 

^(o o- o 'o o'" o J' (5 - 20) 

is seen to be precisely W_v[0, 1, . . . , 2m — 1; 0, 1, . . . , 2n — 1]. 
Moreover, it follows from (4.2) that 

,. M {W N [{vi + ai, . . . + b m }<; +C!,...,v' n + d„}<]) _ 
w™oo M(W_v[0,...,2m-l;0,...,2n-l]) ~ X2m ' 2n 

x sgn(^ + a,v[ +di, . . . ,i4 + c„,^ + rf„)sgn(i;i + oi,Ui +6i,.. . , w m + a m , w m + 6 m ) 

^ (3/2)„ i+ai (3/2)„ i+6< A (3/2) <+e , (3/2) <+di 
X fJ [2) Vi+ai (2) Vi+hi f = l (1) B , +Cl (!)„,+_, 



Q («j - u. + % - aO( w i ~Vi + dj - h)(vj - v t + bj - Oi)(vj - Vi + bj - h) 

II ( V 'l ~ V 'i + C 3 ~ Ci ^ V 'j ~ V '* + C 3 ~ di ^ V 'j ~ V 'i+ d 3 ~ C ^( V 3 - V 'i+ d 3 ~ d i) 



X 

l<i<j<m 



X 

\<i<j<n 



x _______[_____ / 521 x 

Ilili nj=i( u ii + a i + c j)( u ij + a « + dj)(u tj +bi + Cj)(uij +h + dj) ' 

with x given by (4.3) and u;j = + v'j + 2, _ = 1, . . . , m, j = 1, . . . , n. Indeed, (5.21) is 
a direct consequence of (4.2) when v\ + a\ < v\ + b\ < ■ ■ ■ < v m + a m < v m + b m and 
v[ + c\ < v[ + di < • • • < v' n + c„ < v' n + d n - Since each "elementary move" changing the 
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relative order of just two consecutive elements in either of these two lists (except when 
these two elements are Vi + at and Vi + bi, or v\ + c, and v^ + di, for some i, case in which 
their relative order is fixed by construction) has the effect of multiplying both the sign part 
and the product part of the right hand side of (5.21) by — 1, it follows that (5.21) is valid 
in general. 

Therefore, divide (5.19) by the matching generation function of the region (5.20) and 
let N — > oo. By (2.8), the left hand side becomes the left hand side of (5.1), while by 
(5.21) the summand on the right hand side becomes precisely the summand on the right 
hand side of (5.1). (In particular, it follows that the limit (2.6) defining u>b exists.) We 
obtain 



Rl Rm . R\ R' n 

v 1 v m ' v[ v' 



2 m+n x m ,„ n ^ n R i( R i - V2)(^ + 1/2) 



■■■ H ••• Y (-i)S™i( a ' + ^)+s"=i( ci + d ^ 

0<ai<6i<fli <3<a m <b m <R m 0<c 1 <d 1 <R' 1 0<c n <d n <R' n 



n 



(Ri + a, - 1)! {Ri + bi- 1)! (3/2)„ i+04 (3/2)„ i+6( 



L = l (2a,)! (Ri - ai)\ (26,)! (ifc - 6,)! {2) Vi+ai {2) Vi+bi 

yr (Ri + a-iy.jRi + di-iy. (3/2)„> +Ci (3/2) <+di 

X 11 (2c, + 1)! (Rt - a)\ (2d, + 1)! (Ri - di)\ (l) v > +Ci (l) vl+di 

x } [ (vj -Vi + ctj - cti)(vj -Vi + aj - bi)(vj -Vi + bj - Oi)(vj -Vi + bj - h) 

\<%<j<m 

X II ( W J ~ V 'i + C 3 ~ C »)K' ~ V 'i + C l ~ d i)( V 'j - V 'i + d 3 - c »)(«j - v 'i + d J - di) 

l<i<j<n 

Y[ZM-t>i?K=M-di? 



Tlili Hj=i( u ij + a i + c j)( u ij + a i + dj)(iiij + hi + Cj)(uij + bi + dj) ' ^ ^ 

where the summation range is restricted to those summation variables for which vi + 
ai,vi +bi,...,v m + a m , v m + b m , as well as v[ + c\, v[ + di, . . . , v' n + c„, v' n + d n , are 
distinct. 

The fortunate situation is that, on the one hand, when a, = bi, for some i — 1, . . . , to, 
or when Ci = di, for some i = 1, . . . ,n, the summand in (5.22) becomes zero. On the 
other, the summand in (5.22) is invariant under independently interchanging a, with bi, 
i = 1, . . . , m, and Cj with di, i = 1, . . . , to (because the differences bi — cii and di — Ci end up, 
by the combination of (5.5), (5.13) and (5.21), appearing at the second power). Therefore 
the summation range may be extended in (5.22), at the expense of a multiplicative factor 
of l/2 m+ ™, to the summation range in (5.1). This leads precisely to the multiple sum 
given in the statement of the Lemma. □ 

6. Separation of the (2m + 2n)-fold sum for uii, in terms of 4mn-fold integrals 



If it weren't for the denominator in the last fraction on the right hand side of (5.1), one 

33 



could expand the product at the numerator of the summand in terms of monomials in the 
summation variables, and the (2m + 2n)-fold sum (5.1) would separate: one could express 
it in terms of simple sums. 

We can get around this obstacle by expressing each factor at the denominator as an 
integral: 



Uij 


+ a t 


+ Cj 




1 




Uij 


+ di 


+ dj 




1 




Uij 


+ h 


+ Cj 




l 




Uij 


+ h 


+ dj 



Uij 1 7 

X ij ax ij 




y-r-'" 1 ** l d Vij 




Uij+bi+cj-l , 
*ij 






(6.1) 



for i = 1, . . . , m, j — 1, . . . , n. Then the multiple sum inside the absolute value signs in 
(5.1) can be expressed as a sum of 4mn-fold integrals of products of simple sums. More 
precisely, expand the product of the numerators on the last three lines of (5.1) as 

m n 

s :=n^-^) 2 n^-^) 2 

i=l i=l 

x J| ((vj - v{) + aj - a,i)((vj - v t ) + a 3 - bi)((vj - v t ) + bj - a l )((v j - v{) + bj - b { ) 

1<z<j <rn 

X 11 {{v'j - V'i) + Cj - Ci){{v'j - V'i) + Cj - di)((v'j - v'i) + dj - Ci )((v'j - v'i) + dj - di) 
l<i<j<n 

= Y, e(C)aT {C) b^ C) ■ ■■a^ c ht (C) ^ {C) d 5 ^ C) ■ ■■c^ c U s n ^ c \ (6.2) 
cec 

where C is the collection of all 2 2m + 2 "3 4 ( 2 ) +4 W signed monomials in the (vj — «i)'s, 
(v'j — w-)' s j ai' s , bi's, Cj's and dj's obtained by expanding the left hand side of (6.2), and 
for such a monomial C <G C, ai(C), (3i{C), 7j(C) and Sj(C) are the exponents of a i: b i7 Cj 
and dj, respectively, while e(C) is the part of C besides Jlti a* 1 ^ b^ ^ n£=i Cj 3 ^ d^ 3 ^ 
(so e(C) is a signed monomial in the (vj — «i)'s and (v'j — v-)'s). 
Define also 

1 ^ (-R) a (R) a (3/2) v+a f xy n 
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n>(n)fn „..„"> ._ ly(--R)c(fl)c(3/2),+c (X 

We have the following result. 



Proposition 6.1. The boundary-influenced correlation uo^ can be written as 



' Jb ( t ' ' ' I! 5 5 ' ' ' < ) = X2m ' 2 " A ^ ft ^ - V2)(^ + 1/2) 

n i— 1 z— 1 



^ e(C)M Ql(c)A(c)i ... ;am(c . ))/ 3 m(c) . 7l(c)i5l(c)i ... ;7n(c)i5n(c . ) 



(6.5) 



where \ *s <?wen (4.3), i/ie collection C and e(C), ai(C), (3i{C), 7j(C), <5j(C) are as in 
(6.2), and the "moment" M ai;/ a 1; ... iami| g m;7li 5 li ... ;7rli 5 n equals the 4mn-fold integral 

~l ~l m n 

M «i,ft,..,«»,fc;7i,«i,..,7»,s»= / •••/ n n(%-^z, 3 ic 3 )"' + " ;+i 

i=lj . =1 

n n 

x r( ai >(fli, «i; [] aiyyij) T^\R uVl ; J[ z ljWlj ) ■ ■ ■ 

3 = 1 J = l 

n n 



2 = 1 2 = 1 

m m 

■ ■ •T /(7 " ) (i?^,^; Y[x in z in )T'( Sn \R' n ,v' n ; Y\_yinW in ) dx lx ■ ■ ■ dw mn , 

i=i i=i (6.6) 

withT( n \R,v;x) andT' {n) {R,v-x) defined by (6.3) -(6.4). 
Proof. We note first that, since 

(R + a-l){ (-I)- (-R)a(R)a 



(2a)\(R-a)\ R 4° (1)„ (l/2)„ 

and 

(i2 + c-l)! (-l) c (-R)c(R)c 



(2c+l)\(R-c)\ R 4 C (1) C (3/2) C ' 
the sums in (6.3)-(6.4) can also be written as 

fW '^"-S'-" , gl^fr - ,6 - 8 > 
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Denote by S the multiple sum in (5.1). Express the factors at the denominator of 
the summand of S as integrals, using (6.1). Expand the factors of the summand of S 
contained in the left hand side of (6.2) as a sum of monomials in the summation vari- 
ables a\,b\, . . . , a mi b m and c\,d\, ... ,c n ,d n of S, as indicated by (6.2). Bring all 4mn 
integration signs in front of S. We obtain 

r i r i ( Rl R m R'i K 

[ai,6i=0 a m ,6 m =0ci,di=0 c„,d„=0 

X 11 (2a,)! {Ri - a t )\ (26,)! {Ri - h)\ {2) Vi+ai (2), i+bi 

A (R l + c l -l)\(R l + d l -l)\ (3/2)„> +Ci (3/2)„. +d< 

X 11 (2cj + 1)! (i? 4 - cj)! (2d, + 1)! {Rt - d t )\ {1) <+Ci {\) <+di 

m n 



Wi+t^+di+cy + l Uj+j^.+a.+dj + l Ui+ti^ +6i+d.j + 1 

I 1 "I x " Vij Z ij W ij 



i=l j=l 



X e(C)aT {C) b^ C) ■ ..a^bi^cT iC) 4 l{C) ■ ■ .c^#<°>} d*n ■ ■ ■ d«,„ 

CdC ) 



X 

cec 



Reversing the order of the two multiple summations in the integrand above yields 

ii=0 c„,d„=0 

(Ri + ai-l)\{Ri + bi-l)\ 



,1 r i ( Ri R m R'i K 

•/-/ E X - E E ••• E(-^ : "' K 

JO JO CgC I ai;bl=0 o m ,6 m =0ci,di=0 c„,d„=0 



X n (2a l )!( J R t -a 4 )!(26 J )!(i? 4 -fe 4 )! 

pr (fli + Ci-l)!(fli + di-l)! 

X 11 (2c« + 1)! {Ri - a)\ (2d t + 1)! {Rt - di)\ 

y (3/2)„ i+ni {3/2) Vi+bi pr (3/2)„> +Ci (3/2)„, +di 

X fJ (2)„ (+0( (2)„, +6i 11 {l) vl+ct {l) vl+dt 
m n 

TT TT t, i+ 1 ';+ a i+ c j+ 1 "i+"j+a»+dj+l fi+'f^+bi+Cj+l ^i+Wj+bi+dj+l 
1111 ' % z ij "'<./ 

i=l j=l 

xa o 1 (C) 6 / 3l (C) . . . a ^(C) h r^(C) cT {C) d ^C) . . . c 7n(C)^„(C) j e(cr) . . . ^ 



Clearly, the inner multiple sum separates in terms of simple sums on oi, b\, . . . , a m , b m 
and ci, di, . . . , c n , d n . Moreover, all the simple sums arising this way have one of the forms 
(6.7) or (6.8). Moving the 4mn-fold integral sign inside the outer multiple sum above one 
obtains (6.6). □ 
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The sums (6.3) and (6.4) can conveniently be expressed in terms of hypergeometric 
functions 6 . Indeed, for integer R the upper summation limits in (6.3) and (6.4) can be 
replaced by oo without affecting the definitions, due to the factors (—R) a and (— R) c m 
the summands. Using that (x) v+a = (x) v (x + v) a , (6.3) becomes 



{*,v,x)- R {2)v ^ {1)a{1/2)a{v + 2)a UJ 



1 (3/2)™(-R) a (R) a (v + 3/2) a 



R (2)„ 



a=0 



(l)-(l/2) 



a (v + 2) a UJ 



x {/ n a(a - 1) • • • (a - n + 1) + f n -\a{a - 1) • • • (a - n + 2) H h /ia + Jo}, 

(6.9) 

where the coefficients /j are defined so that the last factor in the summand equals a n (in 
particular, /„ = 1 and /o is the Kronecker symbol S n o). Since a(a— 1) • • • (a — fc + l)/(l) a = 
l/(l)a-fe 5 for a > fc, we have 



(-■R), (fl)„ (v + 3/2) „ fa^« 
^ (l)„(l/2) (t; + 2) a 



E 



(^) a a(a-l)...(a-fc + l) 
(-iZ)fc (i2)fc (v + 3/2) fc /z\ fe y, + fc ) a _ fc (g + fc ) a _ fc (t, + fc + 3/2) a _ fc /£\»- fc 



(l/2) fc (« + 2) fc 



(-iZ)fc (i2)fc (« + 3/2) fc /rex* 



(l/2) fc (« + 2) fc 



(9 



i—k 



3F2 



-R + k, R + k, %+v + k x 



± + fc,2 + v + fc 



Substituting this into (6.9) we obtain the first part of the following result. 
Lemma 6.2. We have 



T^ n \R,v;x) = 

1 (3/2) " (-R) k (R) k (v + 3/2) h 



R (2) v 



fc=0 



(l/2) fc (w + 2)fc 



3^2 



n/(") 



(R,v;x) 



1 (3/2), " (-i?) fe (i?) fe (« + 3/2) fc /x 



i? (1). 



E/< 



(3/2) fc (w + l)fc 



/ x\ 

(4) » Fa 



-i? + fc, i? + fc, §+-u + fc_x 
± + fc,2 + u + fc '4 



-i? + fc, i? + fc, § + v + fc _ 
§+fc, l + v + k 



(6.10) 



(6.11) 



6 The hypergeometric function of parameters ai , . . . , a p and bi,...,b q is defined by 

01 , . . . , Op 



61, ■ • • ,b t 



where (a)o := 1 and (a) fe := a(a + 1) • ■ ■ (a + fc — 1) for k > 1. 
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where the fk 's are as in (6.9) (in particular f n = 1). 

Proof. Starting from (6.4), we obtain (6.11) by nearly the same calculation that proved 
(6.10). □ 

Since by Proposition 6.1 the boundary-influenced correlation w& is expressed in terms 
of the moments M aii/ g lj ... iami| g m;7li 5 li ... ;7nj 5 n given by (6.6), which in turn depend on 

T^(R, v; x) and T'^ n \R, v;x), to determine the asymptotics of u>b we need to under- 
stand the asymptotics of the T^'s and T'^'s. In the next section we deduce these latter 
two asymptotics from a result — stated in Proposition 7.2 — whose technical proof we defer 
to Section 9. In Section 8 we show that the asymptotics of the T^'s and T'^'s found in 
Section 7 can be used to obtain the asymptotics of M aii| g li ... iami( 3 m;7li 5 1; ... i7n; 5 n . 

7. The asymptotics of the TW's and T' (n), s 

Given that in Theorem 2.1 the coordinates of the quadromers approach infinity as 
indicated by (2.3), we need more specifically to find the asymptotics of T^(R, qR + c; x) 
and T l<<n \R, qR + c;x), as i? — > oo. These are given by the following result. 

Proposition 7.1. Let q > be fixed rational number, and let n > and c be fixed 
integers. Then for any real number x G (0,1], we have 



T^ n \R,qR + c-x 



X COS 



T' (n) {R, qR + c; x) 



i?arccos ^1 — — ^ — - 
1 



R 3/2 



R 



arctan — 

q 



X COS 



R arccos 



('-f) + 5 



arctan — 

q 



R 



nir 
~2~ 



< MR n ~ 5/2 



(7.1) 



x (n — 1)tt 



< M'R n - 5 / 2 , 
Vx 



(7.2) 



for R > Rq, where i?o, M and M' are independent of x G (0, 1]. 

In our proof of the above statements we make use of the following result, whose proof 
is presented in Section 9. 

Proposition 7.2. Let p(t) and Q(t) be functions depending on the parameter x G [0, 1], 
defined on (0, 1) by 



Pit) 

Q(t) 



-qlnt 

t l 



i i xt 

i arccos 1 



(4 - 2xt) a 



(l-t)V2 (4-xt) b 
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(7.3) 
(7.4) 



where 0<geQ,0<agZ, -1/2 < b £ \l and I £ \L are all fixed. Then 

r-l 



-Rp(t) 



Q(t)dt- F(R,x) 



< MR~ 3/2 , 



for R > Ro, with Rq and M independent of x £ [0, 1], and F(R,x) given by 



F (R,x) = ^ (4 ~ 2x)a/(4 = X) efo""*"* 1 -*)-* 



arctan -| \f-%z 



(7.5) 



(7.6) 



COROLLARY 7.3. Under the assumptions of Proposition 7.2, we have 



-Rp(t) 



Q(t)dt 



< MRT 1 ! 2 , 



(7.7) 



where M is independent of x £ [0, 1]. 

Proof. The absolute value of the part of the right hand side of (7.6) not containing R 
can clearly be majorized, for x £ [0, 1], by a constant independent of x. The statement of 
the Corollary follows by combining this observation with (7.5). □ 

Proof of Proposition 7.1. By [24, (10), p. 58], a p+ \F q+1 hypergeometric function can be 
written as an integral of a p F q as 



p+l-Fg+l 



(3, a P 
[3 + a, p q 



r(/3 + a ) 

r(/3)r( 



4 / t?-\i-ty-\F q 

o) Jo 



Up ;zt 
Pi 



dt, (7.8) 



provided p < q + 1, Re/3 > 0, Reer > 0, and \z\ < 1 if p = q + 1 (here a p stands for a 
p-tuple, yO g for a q-tuple of parameters). 

Taking (3 = v + k + 3/2 and a = 1/2, (7.8) yields 



3^2 



-R + k, R + k, | + v + k_ x 
\ + k,2 + v + k '4 



l> + fc + 2) 



r(« + fc + 3/2)r(i/2) 



x / t v+k+1 ' 2 {l-t)- 1 ' 2 2 F 1 
Jo 



-R+k, R + k xt 



dt. (7.9) 



On the other hand, from the definition of a 2 F\ hypergeometric function it readily follows 
that 



d F 



Repeated application of this shows that 



ai, a 2 


a\a 2 


b ' ,Z _ 


— , 2-Fl 

b 



iF\ 



-R + k, R + k xt 
\ + k ; 4 



(1/2)* 



ai + 1, a 2 + 1 
6+1 



{-R) k (R) k (x/A) k dt k 
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2*1 



-R, R xt 
1 > T 



(7.10) 



However, the latter 2*1 evaluates exactly (see for instance [14, p. 1055, #1]): 



2*1 



— R. R 

1 \z 
2 



cos [i?arccos(l — 2z)] . 



(7.11) 



Replacing (7.9) and (7.10) into (6.10) and using (7.11) we obtain 



(-fl) fc (R) k (3/2) v+k ( x \ k T(v + k + 2) 
T (n )m „. x _ 1 V , (l/2)fc(2)„ +fc UJ I> + fc + 3/2)r(l/2) 

UJ (l/2) fc 



fc=0 



r r+fe+i/2(i - fri/2 ^ c ° s 

1 r +fe+l/2 (1 _ f) -l/2^ cos 



i? arccos 




H)] 


i? arccos 









eft 



dt 
(7.12) 



(for the second equality we also used T(l/2) = y/n and the recurrence T(a; + 1) = xT(x)). 

By Lemma 7.4, the successive derivatives with respect to t of cos(i?arccos(l — xt/2)) 
can be written as 



dt k 



cos (R arccos (1 — xt/2)) = cos 



R arccos 1 



xt\ nir 



R 



0{R k - 1 ) 



At - xt 2 



(7.13) 



where the terms of the omitted linear combination are of the form R^x s cos(i?arccos(l — 
xt/2)) or R^x s sin(7?arccos(l — xt/2)) times a function of type (7.4), with < j < k — 1 
and s > 0. (In fact, this is how the family of functions (7.4) was chosen, to contain all 
functions arising this way from successive derivatives of cos(i?arccos(l — xt/2)), and the 
analogous functions resulting when doing the same analysis to the T'^'s.) 

Since clearly \x\ < 1 throughout the range x £ (0,1], one sees by Corollary 7.3 that 
for v = qR + c, q > and for any fixed k € {0, . . . , n}, the total contribution of the 
lower order terms in (7.13) to the sum on the right hand side of (7.12) is bounded in 
absolute value by Lfci? fe_3 / 2 for R > Rk, where Rk and Lk are independent of x e [0, 1], 
k = 0, . . . , n. Therefore, the combined contribution to T^> (R, qR + c; x) of all lower order 
terms in (7.13), for k — 0, . . . ,n, is bounded in absolute value by LR k ~ 5 / 2 for R > po, 
where po and L are constants independent of x e [0,1]. 

On the other hand, by using Corollary 7.3 as in the previous paragraph, one sees 
that the combined contribution of the leading terms of (7.13) for k = 0, . . . , n — 1 to 
T( n \R,qR + c;x) is again bounded in absolute value by KR n ~ 5 / 2 for all R > pi and 
x e (0, 1], for some constants K and pi independent of x e (0, 1]. 
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Taking into account the leading term of (7.13) for k — n, we obtain by (7.12) and the 
previous two paragraphs that 



T^\R,qR + c;x)-^J^t qR 



j-n+c+1/2 
(1 -i) 1/2 



R 



At - xt 2 



X cos 



xt 



\ nir 


dt 


) + T 







< M R 



n-5/2 



(7.14) 



for R > p2, with p2 and Mo independent of x € [0,1]. However, it is readily seen by 
Proposition 7.2 that 



I 



-1 t n+c+l/2 
t qR 



o (l-i)V2 VV 4t-xt 2 



cos 



i? arccos ^1 — 



4 - x 



-1/4 



4 - x 



X COS 



i? arccos ^1 — 
< Mii?"" 5 / 2 , 



1 i 

arctan — , / 

2 gV 4 -^ 



+ ■ 



717T 



n7r 



(7.15) 



for R> p 3 , with p 3 and Mi independent of x <G [0, 1]. Relations (7.14) and (7.15) imply 
(7.1). 

We now turn to proving (7.2). Because of the requirement Recr > we cannot apply 
(7.8) directly to the 3F 2 's of (6.11), so we write first 



3^2 



-R + k, R + k, | 
| + 1 + v + k 



v + k x 
5 4 



E 

a>0 



(-R + k) a (R + k) a (1/2 + v + k) a (l + l/2 * v + k 



3 r 2 



a\ (3/2 + k) a (l + v + k) a 
1 t v t ft; 

' 4 

(-R + fc)i (i2 + fc)i (1/2 + w + fc)i /x 



(!)' 



-R + k, R + k, ± + v + k x 
^ + k, 1 + v + k 



+ 



(!) 



1/2 + w + fc (3/2 + fc)i (1 + v + k)i 

{-R + k + !)„_! (i? + fc + 1) _! (3/2 + w + fc) _! 



E 

a>l 



3-^2 



(a- 1)! (5/2 + fc) a _ 1 (2 + W + fc) a _ 1 
+ 



(!)' 



-R + k, R + k, \ + v + k x 



x (-R + k){R + k) 
4(3/2 + k)(l+v + k) 



3F2 



' 4 



-i? + fc + l, ii + fc + l, f+w + fcx 



+ fc, 2 + v + k 



(7.16) 
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The two 3-F2 's on the right hand side of (7.16) have the same form (the second is obtained 
from the first by replacing k by k + 1). Equality (7.8) is applicable to them and yields 



3^2 



-R+k, R + k, \+v + k x 
| + l + v + k '4 



1) 



T(v + k + l/2)r(l/2) 



f 1 t v+k-l/2 {1 _ t yl/2 2Fi 

Jo 



-R + k, R + k xt 



§ + * 



dt. 



Repeated application of the relation just before (7.10) shows that 



2F1 



-R + k, R + k xt 
% + k ; T 



(3/2)* 



(-R) k (R) k (x/4) k dt k 



iF 1 



-R, R xt 

3 ; t 



(7.17) 



(7.18) 



To continue our analysis we need a closed formula for the 2F1 on the right hand side of 
(7.18). We obtain it as follows. By [14, 8.962, #1] one has 



P^\o 



(-!)"!>+ 1 + /3) 
n!T(l + /3) : 



-n, n + a + f3 + 1 1+x 
1 + /3 



where P^ ,,i \x) is the nth Jacobi polynomial of parameters a and (3. Taking a 
fi = 1/2, this gives 

n!r(3/2) D (-3i, 



2F1 



-n, n 1 + x 
3 I 7. 
2 z 



(-!)"!> + 3/2) 



\x). 



-3/2, 
(7.19) 



On the other hand, from [14, 8.962, #4] we obtain 

P { n hh \x) 



r(2)r (n + 3/2) ci 



T(n + 2)r(3/2) 
T(2) T(n + 3/2) sin[(n+ l)arccosa;] 
~ r(3/2) T(n + 2) sin(arccosa;) ' ( ' ' 

where is the nth ultraspherical polynomial of parameter A, and at the last equality in 
(7.20) we used [14, 8.937, #1]. 

However, by [14, 8.961, #8] the Jacobi polynomials satisfy the recurrence 

(2n + a + f3)Pi a - 1 ^(x) = (n + a + 0)P^\x) - (n + /?)p£f \x). 

By this the explicit formula (7.20) for P n 2 ' 2 {x) yields one for P n 2 ' 2 (x), which in turn, 

by another application of the above recurrence, yields an explicit formula for P n 2 ' 2 (x). 
Substituting this into (7.19) one obtains after simplifications that 



2*1 



— n, n t 
3 ; 7 



2n 



4n 2 - 1 



t 



■ sin 
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/; an ci is ( 1 — ^ 



4n 2 



• cos 



n arccos 1 — 



(7.21) 



Expressing the 3 F 2 's in (6.11) with the use of (7.17), (7.18) and (7.21), (6.11) becomes 
T' (n) (R,v;x) = 

1 ™ 2U , ; [-) !.'(,■ + /,- l."2il.'(i/2! 



k=0 
~1 



(-R) k (R) k (3/2)„ +fc (x\ k T(v+k+l) 

v+k \i) r(t)+ fc+i/2)r(i/2) 
tt; ; ~ (f) 



(-R) k (R) k (x\ k 



2R 



4-xi 



4R 2 - 1 



AR 2 - 1 V xt 
_/t arccos 1 



■ sin 



xt 

i? arccos 1 — — 



dt 



n 



(-R) k (R) k (3/2) v + k (x \k I>+fc+2) 

(3/2) fc (i)„ +fc U/ r( v +fe+3/2)r(i/2) a; (-R + k)(R + k) 



k=0 

"1 



(-fl)fc+l (fl)fc+l (x^ 1 

(3/2) fc+1 U/ 



/ ^+fe+l/2 (1 _ t )-l/2 
JO 



2R 



d k+i 



dt^ 1 ) 4R 2 - 1 



1 



4.R 2 - 1 



// arccos 1 



4(3/2 + fc)(l + i; + fc) 



4 - xt f xt 

sin R arccos 1 — — 

xt \ 2 

xt , , 

dt, 



which when substituting v = qR + c becomes after simplifications 

a « fl .fc+c-l/2 



T' (n) (i?, gi? + c;x) 



dt k 



4 - xt 



xt 



■ sin 



i? arccos 1 — 



1 

cos 

2R 



R arccos 1 — 



xt 



+ nAR 2_ 1 2^hJ o 1 (1 _ t) l/2 



d fc+1 



4 -art 
xt 



■ sin 



i? arccos 1 



xt 



2i? 



cos 



i? arccos 1 



fit 



| dt. 

(7.22) 



Let.F(i) = ^(4 - xt)/(xt) sm(i?arccos(l-xf/2))-l/(2i?) cos(i?arccos(l-xt/2)). Lemma 
7.4 implies that 



jp(k)( t ) = 



4 -art 



xt 



■ cos 



xt\ (n — 1)tt 



/? 



4t - xt 2 



+ 0(i? fc - 1 ), 
(7.23) 



where each term of the omitted linear combination is of the form R J x s cos(i?arccos(l — 
xt/2)) or Rix s sin (R arccos (1 — xt/2)) times a function of type (7.4), with j < k — 1 and 
s > -1/2. 
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Since there are only a finite number of omitted terms in (7.23) for each k = 0, . . . , n, 
we obtain by (7.22) and Corollary 7.3 that their contribution to the first sum in (7.22) 
is bounded in absolute value by L k R k ~ 3 ^ 2 / \fx for all R > R' k and x G (0,1], where 
L' k and R' k , k = 0, . . . , n, are some constants independent of x G (0, 1]. Therefore, the 
overall contribution of the non-leading terms in (7.23) to the first of the two terms on the 
right hand side of (7.22) is bounded in absolute value by L' 1 R n ~ 5 / 2 /y/x for all R > p' 
and x € (0,1], where L\ and p' are constants independent of x G (0,1]. Similarly, the 
overall contribution of the non-leading terms in (7.23) to the second term on the right 
hand side of (7.22) is seen to be bounded in absolute value by L' 2 R n ~ b / 2 / 'y/x for all 
R > p" and x G (0, 1], with constants L 2 and p" independent of x € (0, 1]. Therefore the 
total contribution of the non- leading terms in (7.23) to T' {n) (R,qR + c; x) is bounded in 
absolute value by L'R n ~ 5 / 2 / ^/x for all R > p' and x G (0, 1], for some constants V and 
p' independent of x G (0, 1]. 

On the other hand, by using Corollary 7.3 as in the previous paragraph, one sees 
that the combined contribution of the leading terms of (7.23) for k = 0, . . . , n — 1 to 
T'( n \R, qR + c; x) is again bounded in absolute value by K'R n ^ 5/2 /^/x for all R > p\ and 
x G (0, 1], for some constants K' and p' x independent of x G (0, 1]. 

By the previous two paragraphs and (7.23) we obtain from the k = n terms of (7.22) 
that 



T' (n) (R,qR + c;x)-{ / t« R 



t - 1 / 2 U-xt 



wR J (1 - i)V2 V xt 

xt 



R 



xt 



X cos 



R arccos 1 — 



+ 



4-xt 
(n — 1)tt 



dt 



+ 



nR 2 l 



t qR 



(l-t)V2V Xt 



t - 1 / 2 U-xt 



R 



xt 



n+l 



A-xt 



X cos 





nix 




M) 






+ T 





< —=M' R n - 5 / 2 1 



(7.24) 



for all R> p' 2 and x G (0, 1], where the constants M' and p' 2 are independent of x G (0, 1]. 

Denote by I\ and li the first and second integrals on the right hand side of (7.24), 
respectively. They clearly satisfy the hypothesis of Proposition 7.2. By (7.5) and (7.6) we 
obtain 



q 



,2 



^ VnR 
x exp I i 



-1/4 



A - x 
R arccos ^1 
< MiR^^x^- 1 



R 



A- x 



x 
2 

)/2 



1 1 

— arctan — 

2 q 



(n — l)ir 



dt 



(7.25) 
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and 



irR 2 



h - Re 



/ L 



\ V 4 - 



-1/4 



J? 



4 - x 



x exp z 



R arccos ( 1 



x\ 1 



1 / x 



■ arctan ■ . 
2 g\/ 2 



< M! 2 R n ~ b l 2 x n/2 , 



(7.26) 



for R> p' 3 , where P3, M[ and are constants independent of x € (0, 1]. 
Using the fact that 



q + l 



A-x 



4 - x 



exp 1 arctan - 



1 



q V 4 — x 



the two second terms on the left hand sides of (7.25) and (7.26) are readily seen to add 
up precisely to the second term on the left hand side of (7.2). Since in (7.25) and (7.26) 
the integer n is nonnegative and since |x| < 1, inequalities (7.25), (7.26) and (7.24) imply 
(7.2). □ 

Lemma 7.4. Let the functions F\(t), F 2 (i) and F 3 (t) be given by 





R arccos 




= cos 






R arccos 




= sin 


H)] 









A-xt 



xt 



■ sin 



R arccos I 1 — ^ 



where x € (0, 1] is a parameter. Then for any n > one has 
F[ n \t) = (r 



Ft\t) = [ n 

Ft\t) = 



A-xt 



At - xt 2 
x 

At - xt 2 
R 



cos 



sin 



xt \ mr 
it arccos | 1 — — | + 



xt\ mr 
R arccos [1 — — | + 



+ y: cvr^x^qv 

(7.27) 



(7.28) 



cos 



xt \ ~ V At - xt 2 

+ cWr^x^QP 
uev 3 



R arccos 1 



xt \ (n — l)ir 
' + " — ^ — — 



(7.29) 



where l^\l^ > 0, ll?' > —1/2, for all v, and for k — 1, 2, 3 and all v G 14 one /ias 
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,(3) 



(i) Vk is some finite set 
(fc) .(fe) ,(fe) ^ 

(m) < n — 1 

is a function of type (7.4) multiplied by either cos(i?arccos(l — xt/2)) or 
sin(i?arccos(l — xt/2)). 

Proof. It readily follows by induction on n that for n > 1 one has 

^e h W=e h WJ2 E ajh^\t)---h^(t), (7.30) 

fc=i »i — »fe>i 

il + - + i fc =n 

with coefficients a/ = ari^...,^) £ Z and an...,!) = 1. 
A similar argument shows that for n > 1 one has 

^[5(*)]" 1/2 = I.9Wr 1/2 E E /3/bW]- 4 °.9 (n) W---3 fe) W, (7-31) 

fc— 1 io,ii,...,ifc>l 

where only finitely many of the coefficients f3j = P(i Qt i 1 i k ) € Q are nonzero. 

Choosing h(t) = ii?arccos(l — xt/2), the derivatives on the left hand side of (7.27) can 
be expressed as 

F<»>(t)=Re|^«. (7.32) 

For our choice of h we obtain 



h'(t) = iR^j^—p = iR^lgit)]- 1 / 2 , (7.33) 

with g(t) = At — xt 2 . Since g'(t) = 4 — 2xt, g"(t) — — 2x, and the higher derivatives of g 
are zero, (7.33) and (7.31) imply that for j > 1 we have 

where the sum is finite, V £ Q and a v ,b v ,s v > are integers. From the expression of 
h'(t) it is apparent that (7.34) holds also for j = 0. 
Substituting (7.34) in (7.30) we obtain for n > 1 that 

in n ~ 1 

±- e Ht) ^ e Ht) [h > {t)] n + e h(t)j2 J2 a T h^\t)---h^\t) 

fe=i «i,...,i fc >i 

i 1 + ... + i fe = n 
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where the sum is finite, n — 1 > n v G Z, < a„, s„ € Z, < 6^ € ^Z, and 7^ G Q. Relation 

(7.35) clearly holds also for n = 0. By (7.32), taking the real and imaginary parts in (7.35) 
one obtains (7.27) and (7.28), respectively. 

To prove (7.29), apply Leibniz's formula for the derivatives of the product /i(i)/2(i) 
defining F 3 (t), with fi(t) = ^/(4 — xt)/(xt), f 2 (t) = sin(i?(arccos(l — xt/2)). It is clear 
by (7.28) that the highest power of R in the resulting terms is n, and it occurs in only 
one of these terms, namely in /i(i) times the highest order term in /^(t). By (7.28) it 
also follows that the leading term in R of F^ n \t) is the first term on the right hand side 
of (7.29). 

Furthermore, (7.28) implies that the successive derivatives of the second factor f 2 {t) 
have the form of the summand on the right hand side of (7.29). Therefore, to finish the 
proof it is enough to show that the derivatives of f\(t) are also of this form. 

We have 

d n U-xt 1/2 d n [4 

- — \ = x~ 1/2 — \ - -x. 7.36 

dt n V xt dt n \ t y ' 

One easily obtains by induction on n that for n > 1 

^b(i)] 1/2 = b(i)] 1/2 E E 6il9(tT i °9 iil \t)---9 iik Xt), (7.37) 

k=l i ,i 1 ,...,i k >l 

where the sum is finite and the coefficients Si = /3(i .i lj ....i fc ) € Q. 

Choose g{t) = \-x. Since g^{t) = t^'" 1 , j > 1, it follows by (7.37) and 

(7.36) that for n > 1 



/A \~ s " 

EE 5 - (7-^ ^ ^ 7 - 38 ) 



k=l v 

where < s„, l v £ Z. Since the terms in (7.38) have the form of the summand on the right 
hand side of (7.29), and since this is also true of the left hand side of (7.38) for n = 0, the 
proof of (7.29) is complete. □ 

8. Replacement of the T^'s and T'^'s by their asymptotics 
does not affect the asymptotics of the moments M 

Denote the approximants of the T^'s and T'^'s in Proposition 7.1 by 
F k (R,q;x) = (r 

{l 2 + 4^) 



X cos 



/ x\ 1 1 / x kir 

R arccos 1 arctan - * / 1 

V 2/2 q\l 4-x 2 
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(8.1) 



and 



F' k {R,q-x) = 



1/4 



R 3/2 



R 



4 - x 



x cos 



R arccos 



(, x \ 1 1 f 

1 — — + - arctan - . / ■ 
V 2/ 2 q\ 



i r~r~ t (k- i)tt 

g V 4 - a; 



+ 



(8.2) 



Proposition 8.1. Consider the moments M aij/ g lj ... iami| g m;7li 5 li ... ;7n; 5 n defined by (6.6). 
^4s i/ie variables Ri, Vi, Rj, v'j, i = 1, . . . , m, j = 1, . . . , n approach infinity as specified by 
(2.3), we /iai>e 



M, 



/■l »i m n 

i,/Ji,...,a ra ,/J rai 7i,«i,...,7»,«»= / •••/ II XX^im^ij^iiY'' 1 ''' 1 '^'^ '''' J 

Jo Jo i=lj=1 

n n 

x F ai (i?i, gi; JJ xijyij) F [h {Ri,qi; J| zytoij) ■ ■ ■ 

3 = 1 3 = 1 

n n 

m m 

x ^ (izi , gi ; n z *i ) (^i » «i ; n ) • • • 

i=l i=l 

■ • • F 7n (R' n ,q' n ; Y[ XinZin) F' 5n (R' n ,q' n ; J} yinW in ) dx n ■ ■ ■ dw mn 

+ ^E™l( Q «+ft)+i:*Li(7 I +'5 I )-3m-3r l -l^ ^ 

w/iere the integration variables are Xij, yij, Zij, Wij, i = 1, . . . , m, j = 1, . . . , n. 
Proof. By Proposition 7.1, for any fixed k > we have 

\T {k) {R,qR + c;x) -F k (R,q;x)\ < MR k ~ 5 ' 2 
\T' {k) (R, qR + c;x)- Fj,(R, q;x)\< x~ l ' 2 M' R k ~ 5 / 2 



(8.3) 



(8.4) 
(8.5) 



for R > Ro, where the constants i?o, M and M' are independent of x € (0, 1]. 

From (8.1) and (8.2) it is clearly seen that there exist constants Mi and M{ so that 



\F k (R,q;x)\ < M^" 3 / 2 
\F} t {R,q;x)\<x- 1 '*M' 1 R!'-W, 



(8.6) 
(8.7) 



for all k > 0, R > and x £ (0, 1]. 
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By (8.4)-(8.7) it follows that for any fixed k > 

\T^ (R, qR + c;x)\< M 2 R k - 3 ' 2 (8.8) 
\T' {k) (R,qR + c;x)\ < i-WM^R*- 3 ' 2 , (8.9) 

for R> R' , for some constants R' , Mi and M' 2 that are independent of x e (0, 1]. 
By (6.6) we have 

~\ ~l m n 

M auPl ,... tami p m71uSu ..., ln , Sn - ■ ■■ [J H(x ij y ij z ij w ij r' R ^ R '^ +c J +1 

Jo Jo j=lj . =1 

n n 

x F ai (i?i,gi; JJ xijyij) F/3 1 {R 1 ,q 1 ; zyWy) • • • 

n n 
' -^a m (-^m; 9m: ^mjUrnj) ^/3 m (-^m: 9m; J^J ^mj^rnj) 



x (^i , ?i ; n ^ ) > ; n ^1^1 ) • • • 

m m 



,1 m n 



/"■■■ / nn(ww) 8W ^ +i 



J=l 3=1 

n n 



x > , «i ; J] IH z a ) T' (Sl } , „J ; JJ ) • • • 

2 = 1 i=l 

m m 

■ ■ ■ T' iln) (R' n ,v' n ; ft x m z m ) T'^\R' n ,v' n ; [] y m w m ) 



F ai (Ri ,qi ,Y[ xijVij ) F /3i ( R i . 9i ; II ■ ■ ■ 

3=1 3=1 

n n 
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x F-n (R'i . q'i ; II XilZil ^> F ^ ( R 'i ' £ ; n ) • • • 

i=l i=l 

m m 

■ • • F j n {R'm In, II ^i" 2 ™) F <L (K> 9n5 II y™"'™)} dx H ' ' ' dw ™n 
i=l i=l 

(8.10) 

Clearly, for any 21 quantities fi, g i} i = 1, . . . , I one has 

fi • • • fl ~ 9i ■ ■ ■ 9l = fi ■ ■ ■ fl-i(fl - 9i) + fi ■ ■ ■ fl-2{fi-i - 9i-i)9i 

+ fi ■ ■■fi-s{fi-2 - 91-2)91-191 H 1- (fi - 91)92 ■■■gi- 

(8.11) 

Apply (8.11) to the expression E in the curly braces on the right hand side of (8.10): we 
have I = 2m + 2n, the /j's become T( fe )'s or T'^'s, and the g^s become F k ''s or F' k 's. 

This results in expressing E as a sum of 2m + 2n products of the form h\ ■ ■ ■ h 2m +2n, 
where exactly one of the hi's, say /i, , is a difference — Fk or T 1 ^ — F' k , all others 
being of the form T^ k \ T' (k) , F k or F' k . 

Furthermore, exactly 2n of the /ij's are of the form T' (fe) , F k , or T l{k) - F k , and the set 
of their x-arguments is always 

{mm m m \ 

(8.12) 
i=l i=l i=l i=l J 

Since the 2n subsets of variables that are multiplied together in the elements of (8.12) form 
a partition of the set of our 4mn variables Xij, yij, Zij, Wij, i = 1, . . . , m, j = 1, . . . , n, 
inequalities (8.4)~(8.9), and the fact that |/ij | can be bounded using the sharper (8.4) or 
(8.5), imply 

m n 

\h! ■ ■ ■ h 2m+2n \ < M 2m +^R^M^)+J:: =1 h i+ S,)-3 m -3n-l -Q J J ,,,,,,„;„„,, , •/*. 

i=l 3=1 

(8.13) 

for all R > Rq, where the constants Rq and M are independent of x £ (0,1]. Taking into 
account also the contribution of the double product in (8.10) to the integrand of (8.10), 
by (8.13) the absolute value of the latter is majorized by 



(2m + 2n)M 2m+2n R^=^ + ^+^=^ +s ^- 



3m — 3n— 1 



m n 

xnn^f^^) ?ifl * + ^ +ci+cJ+i/2 - 

i=1.7=l 



Since by (2.3) the exponent of XijyijZijWij is non-negative, this shows that (8.3) holds 
with implicit constant (2m + 2n)M. □ 
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9. Proof of Proposition 7.2 



For fixed x G [0, 1], the asymptotics of the integral in (7.5) can be readily obtained by 
Laplace's method for contour integrals, as it is presented for instance in [28, Ch. 4, §6.1]. 
However, what (7.5) states is the existence of a uniform error bound for the Laplace 
approximation, for x G [0, 1]. We obtain this by extending the arguments in the proof of 
[28, Theorem 6.1, p. 125] to the case when the integrand depends on a parameter. 

For t e (0, 1), write the factor (l-t) 1/2 in the denominator of Q(t) in (7.4) as (l-t) 1/2 = 
—i(t — l) 1 / 2 , where the square root on the right hand side has its principal determination. 
The integral in (7.5) becomes 

,o 

I(R) = -iJ e- Rp ®Q(t)dt =: -iI(R), (9.1) 
where p(t) is given by (7.3) and 

Q<'>= (t 44-^ - (e(o ' i) ' (9 - 2) 

where as in (7.4) < a G Z, —1/2 < b G ^Z, I G and the square root on the right 
hand side of (9.2) has the principal determination. 

Regard the integral I(R) defined by (9.1) as a contour integral over the path V = [1,0] 
(a line segment) in the complex plane. Choose the principal determinations for all the 
multiple-valued maps in the expressions (7.3) and (9.2) defining p(t) and Q(P). Note that 
for any fixed x G [0, 1] the integral I(R) satisfies the following properties: 

(i) p(t) and Q(t) are independent of i?, single-valued and holomorphic in the pointed 
open disk D{1, 1) = D(l, 1) \ {1}. 

(ii) V is independent of R, and ^(1,0) (i- e -> the path V less its endpoints) is contained 
in i>(l,l). 

(in) For t e D(l, 1), the functions p(t) and Q(t) can be expanded in convergent series 

as 

00 

p(t)= P (i) + Y / Ps(t-iy + \ 

where p(l) = —i arccos(l — x/2) and po = — q — i^Jx/(A — x), and 

00 

Q(t) = Y / qs(t-iy- 1 /\ 

s=0 

where qa = (4 — 2x) a / (4 — x) b and (t — l) 1 / 2 has its principal determination. 

(iv) I(R) converges at absolutely and uniformly with respect to R > —l/q and x G 
[0,1]. 

(v) Re{p(t) — p(l)} is positive when t G (0, 1), and is bounded away from uniformly 
with respect to x G [0, 1] as t — > along V. 
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To avoid interruption in proving Proposition 7.2, we phrase three facts we need in the 
proof as Lemmas 9.1-9.3, and include them at the end of this section. 
Consider the map 1 1— > v(t), \t — 1| < 1, given by 

v(t) = p(t) — p(l) = — q\nt — i arccos ^1 — + i arccos ^1 — ^ . (9.3) 

Let U x be a neighborhood of t = 1 and D a disk centered at v = satisfying the 
statement of Lemma 9.1 (in particular, D is independent of x e [0, 1]). By Lemma 9.1(a), 
(9.3) maps U x conformally onto D. Thus the inverse function v i— » t(v) is also holomorphic, 
and hence t — 1 can be expanded in a convergent series 

oo 

t-l = J2 c sV s , veD, (9.4) 

s=l 

where the coefficients c s are expressible in terms of the p s ; in particular, c\ = l/po- 

By Lemma 9.2, one can choose k e [0, l)nC/, with U as in Lemma 9.1(b), k independent 
of R and independent of x e [0, 1], such that the disk \v\ < \p(k) — p(l)\ is contained in 
D for all x € [0, 1]. Then the portion [l,k] of V may be deformed to make its w-map a 
straight line [0,/C], for all x e [0, 1], without changing the value of the integral I(R) (by 
(i) and (ii)). Making the change of variable v — p(t) — p(l) on this deformed portion [1, k] 
of V we obtain 

J e- Rp(t) Q(t)dt^e- R P^ e- Rv f(v)dv, (9.5) 

JC=p(fc)-p(l), /( v ) = Q(t)^ = |||, (9.6) 

and the path of integration on the right hand side of (9.5) is a straight line. 
By (iii) and (9.4), for v e D, f(v) has a convergent expansion of the form 

f( v ) = gjg^iff^L + 2 -- ; (9 . 7) 

po + 2pi (J27Li c sv s ) + 3p2 (J27=i c * vS ) 

the branch of the square root being the principal one. By the binomial theorem one obtains 
from (9.7) that f(v) can be expressed as a convergent series 

/(«) = a v- 1/2 + a lV 1/2 + a 2 v 3/2 + • • • , 

where v 1 / 2 has its principal value and the coefficients a s can be expressed in terms of the 
p s 's and q s 's as in [28, Ch. 3, §8.1]; in particular, 



where 



(9.8) 



the square root having its principal determination. 
Define fi(v) by the relations /i(0) = a\ and 

f(v)=a Q v- 1 / 2 +v 1 / 2 f 1 (v) (vjtO). (9.9) 

By Lemma 9.3, we may assume that 

|/i(«)|<M 2 , |v| <K, xe [0,1], (9.10) 

where the constant M 2 is independent of x e [0,1]. (Indeed, by Lemma 9.2, k can be 
chosen close enough to 1 so that in addition \JC\ = \p(k) —p(l)\ <p,i£ [0, 1], for the p of 
Lemma 9.3.) 

Using (9.9), rearrange the integral on the right hand side of (9.5) as 

£ e-*>f<v)do = ^1°° e-yy-^dy - ^ jT e^y^'dy + £ e~ R ^h{v)dv 
= a R- 1 ' 2 T(l/2) - ei (R) + e 2 {R), (9.11) 

where 

ei(i?) = a R- 1/2 r(l/2,ICR) (9.12) 
e 2 {R)= I e- Rv v 1 / 2 h{v)dv. (9.13) 

Note that in (9.11) y 1 ^ 2 has its principal value (since y = Rv, and u 1 / 2 does so) and ICR is 
not on the negative half- axis (since Rc/C > by (v)), so the incomplete Gamma function 
in (9.12) also takes its principal value. 

The absolute value of ei(i?) can be bounded as follows. By [28, Ch.4, (2.02), (2.04)], for 
real a the incomplete Gamma function T(a,z) satisfies 

r(a,z) = e- z z a - 1 {l + e(z)} (9.14) 

where 6* = phz (i.e., z — re 10 for some r > 0), (i £ (—n,ir) is arbitrary, 

a(p)= sup ^-^ln|l + i| (9.16) 

p ht=-0 \t\ 

and z is restricted by \6 — /3\ < n/2, \z\ cos(6» - (3) > a{(3). 

Apply (9.14) and (9.15) for our case, a = 1/2, z = KR. Choose /3 = 0. We have 
9 = ph/Ci? = ph/C, so the first condition on z, \B — (3\ < tt/2, is met by property (v). The 
second condition on z is also met, because by (9.16) 

3 

<r(0) = sup- — ln(l + t) < 0. 
t>o 2t 
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By (9.14) and (9.15) we obtain 

T{1/2,JCR) = e- KR {JCR)-^ 2 {l + e}, 

with 

w< 1/2 - 1 



i?|/C|cos(ph/C) 2R{-q\nkY 

for all R > 0. Since ReJC = — q\nk > is independent of x, by the last two relations it 
follows that 

|r(l/2,/CR)| < N'hR- 1 , R> n, a; e [0,1], (9.17) 
for some constants Mi and n independent of x. On the other hand, by (9.8) and (iii) 

(4-2x) a / r^=— \-V2 



(4 - a;) 6 V V 4 - x 

so |a | can clearly be bounded above uniformly for x G [0,1]. By (9.17) and (9.12) it 
follows that 

\ei{R)\< M[R- 3/2 , R>n, (9.18) 

for all x G [0, 1], where the constants M[ and n are independent of x. 

Next, we bound the absolute value of C2{R)- The substitution v = JCt in (9.13) implies 

e 2 (R) = JC 3/2 f 1 e- RKT T X l 2 f x (KT)dT. (9.19) 
Jo 

We have 

Re{-RJCr} = -Rt Re{p(k) - p(l)} = -Rrrjk, 
where rjk — — q\nk > is independent of x e [0, 1]. Using this and (9.10) we deduce 

/>oo 

|e 2 (i?)| < M 2 |/C| 3/2 / e- fl " fcT r 1 / 2 dr = M 2 |/C| 3/2 (i?r/ fe )- 3 / 2 r(3/2). (9.20) 
Jo 

Since 

|/C| = \p(k)-p(l)\ < -qlnk + | arccos(l - xk/2) - arccos(l - x/2)\ 

and the two functions inside the absolute value sign are continuous in x € [0, 1], K, can be 
bounded uniformly and (9.20) implies 

\e 2 (R)\ < M^R- 3/2 , R>0, (9.21) 

for all x € [0, 1], where M' 2 is independent of x. 
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Substituting (9.18) and (9.21) in (9.11) and using | e -^(i) | = l for all R e R and 
x € [0, 1], we obtain by (9.5) 

rk 



J e- R ^Q{t)dt - V^a R- 1/2 e- R ^ 



<M 3 iT 3/2 , R>r u 



(9.22) 



for all x e [0, 1], where M 3 and ri are constants independent of x. 

We now turn to bounding the absolute value of the tail of I(R) omitted by the integral 
on the left hand sides of (9.5) and (9.22). We have 

Re{p(t) -p(l)} = -qlnt > 77 > 0, t e [fc,0), 

for all x G [0, 1], where r/ is independent of x. Therefore for R> ro one has 

Re{i?p(i) - Rp(l)} = {(R - r ) + r } Rc{p(t) - p(l)} 

> (i2 - r )?7 + Rc{r p(t) - r p(l)}. 



We obtain 



-iip(t) 



< e" 



1)1 / ' 

Jk 
ik ' 



- Re{ilp(t)-Jip(l)} 



|Q(t)|dt 



-r p(t) 



\\Q(t)\dt 



(we also used |e- flp W| = 1 and |e^ 0?>( : i > | = 1). Choosing r = -l/q, since | e -^oP(*) | = t <jr 0) 
the last integral is uniformly bounded for 2 € [0, 1] (see property (iv)). We deduce 

fO 



Ik 



-Rp(t) 



Q{t)dt 



<M 4 iT 3/2 , R>r 2 



(9.23) 



for all x e [0, 1], where the constants M 4 and r 2 do not depend on x. 
By (9.22) and (9.23) it follows that 

I(R)-V^a R- 1/2 e- R P (1 '> <MR-^ 2 , R> i? , 

uniformly for x € [0, 1]. This can be rewritten by (9.1) as 

I(R) - {-i)^a R- l ' 2 e- Rp ^ < MR~ 3 ^ 2 , R>R Q , 

for all x € [0, 1], where the constants M and R do not depend on x. 
However, by (9.8) and (iii) 

-1/2 



(9.24) 



a 



(4 - 2x) a 
{A-x) b 



where the phase of the quantity in the large parentheses has its principal value — 7r + 
arctan 1/q^/x/ (4 — x). This implies 



a 



(4 - 2x) a 
(4 - x) b 



-1/4 



4 - x 

This shows that the approximant of I{R) in (9.24) has precisely the expression (7.6), and 
the proof of Proposition 7.2 is complete. 
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Lemma 9.1. Let t v x (t) = v(t) be the map defined by (9.3), all multivalued maps in 
(9.3) taking their principal values. Let x G [0, 1] be arbitrary. 

(a) . There exists a disk D centered at in the v-plane, independent of x G [0, 1], and a 
neighborhood U x containing 1 in the t-plane, so that t i— > v x (t) maps U x univalently onto 
D(0,p), and v' x {t) ^ on U x . 

(b) . D and U x can be chosen in part (a) so that there exists a disk U centered at 1 in 
the t-plane, independent of x, with U C U x , for all x G [0, 1]. 

Proof. We have 

j t v x (t) = -| - T — ^ 0, xe [0, 1], |t - 1| < 1, (9.25) 

so the last condition in the statement of part (a) is met. To complete the proof of part 
(a) it suffices to find S, p > independent of x G [0, 1] so that 7 

v x {t) — vo = has a unique solution in the disk \t — 1| < S, for any v G D = D(0, p). 

(9.26) 

Suppose 8 > is independent of x G [0,1] and satisfies 

v x (t) = has exactly one root in \t — 1| < 5, for all x G [0, 1] (namely t = 1). (9.27) 

Since the map (x,t) \— > |wa;(t)| is continuous and non-zero on the compact set [0, 1] x {\t — 
1| =5}, there exists m > independent of x G [0, 1] so that 

\v x (t)\ > m > 0, |t - 1| = 6, x G [0, 1]. 

Choose p so that < p < m. Then for any v G D(0, p) we have |«o| < P < m < l w a;(t)|, 
for \t — 1| = 5, and Rouche's theorem (see e.g. [1]) implies that v x (t) and v x (t) — vq have 
the same number of roots inside \t — 1| = 5. By (9.27) we obtain that the 5 of (9.27) and 
our choice of p satisfy (9.26). 

To finish the proof of part (a) we need to prove the existence of some 6 satisfying (9.27). 

Let xo G [0, 1] be fixed. Since v Xo (l) — 0, and by (9.25) v' (1) ^ 0, there exists 6 Xo > 
so that v XQ (t) has a unique root in \t — 1| < S Xo . Set b Xo = min| t _ 1 | =5xo |^x (*)l > 0- Write 

v x (t) = v xa (t) + (v x (t) - v X0 (t)). (9.28) 

We have 

Vx{t)— v x„(t) = — i[arccos(l— xt/2)— arccos(l— ccot/2)]+i[arccos(l— x/2)— arccos(l— xq/2)]. 

(9.29) 

Let f{u) — arccos(l — u/2). 

7 The fact that in (9.26) 8 is independent of x is not necessary here, but will be needed in the proof of 
Lemma 9.3. 
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Consider first the case x^ ^ 0. Since /'(it) = (4m — w 2 ) -1 / 2 , its absolute value is 
bounded as long as u is bounded away from 0. By the mean value theorem, it follows from 
(9.29) that there exists an open interval U Xo containing x so that \v x (t) — v Xo (t)\ < b XQ 
for x E U XQ and \t — 1| = S Xa . Rouche's theorem applied to (9.28) shows then that v x (t) 
has a unique root in \t — 1| < S Xo (namely, t — 1) for all x <E U Xo . 

For x = we have 

v x (t) — v Xo (t) — — i[arccos(l — xt/2) — arccos(l — x/2)], 

and by the continuity of f(u) at u = one obtains again that Rouche's theorem is appli- 
cable, and there exists a neighborhood Uq of x — in [0,1] so that v x (t) has a unique 
root in \t — 1| < 5 Xo = 5q for all x e U~o, namely the root t = 1. 

Since [0, 1] is a compact set, it is covered by a finite sub-collection, say U Xl , . . . , U Xn , of 
{U x )xe[o,i\- Then 5 = min™ =1 8 Xi satisfies (9.27). 

To prove part (b), consider the inverse function v h-» t x (v), v G D. Let <5 and p be as in 
part (a). Since the original domain of (9.3) is the disk \t — 1| < 1, we clearly have 5 < 1. 
By (9.25), 

^K(t)) = 1 -^=, xe [0,1], |t-l| <5. 



The denominator above is bounded for x and t as indicated. Furthermore, by part (a), 
for any fixed x € [0, 1], the range of v x (t), \t — 1| < S, contains the disk \v\ < p. It follows 
that there exists A > independent of x so that |^f-(f)| > A, for all \v\ < p and x e [0, 1]. 

For any < r\ < p, consider the circle in the w-plane centered at and having 
radius r\. Since t x (v) is a conformal map, it follows that t x (C v ) is a simple closed curve 
in the £-plane containing t = 1 in its interior. By the above lower bound on |^jf (w)|, all 
points of t x (C n ) are at least at distance A77 from t = 1. Therefore, t x (C v ) must contain 
the disk \t — 1| < Xi] in its interior. Thus, U x = t x (D) = t x (Uo< v < P C v ) contains the disk 
\t- 1| < Xp, for all x e [0,1]. □ 

Lemma 9.2. Let the function p(t) be given by (7.3). Then for any p > one can choose 
k e [0, 1), k not depending on x e [0, 1], so that the disk \v\ < \p(k) — p(l)\ is contained in 
D(0,p), for all x e [0,1]. 

Proof. We have 

— p(l) = — q\nt — i[arccos(l — xk/2) — arccos(l — x/2)], 

and since the first term on the right hand side is independent of x, it suffices to show that 
there exists some k € [0, 1) so that 

I arccos(l - xk/2) - arccos(l - x/2)\ < p/2, x e [0,1]. (9.30) 

For f(u) = arccos(l — u/2) one has f'(u) = (4u — u 2 ) -1 / 2 , u £ (0,1], and thus f(u) is 
increasing on [0, 1]. Therefore the absolute value in (9.30) equals in fact 



g(x) := arccos(l — x/2) — arccos(l — xk/2). 
57 



One readily sees that g'(x) = {Ax - x 2 )- 1 ' 2 - (Ax/k - x 2 )- 1 / 2 > 0, k e (0, 1), x G [0, 1], 
so (9.30) is equivalent to g(l) = arccos(l — 1/2) — arccos(l — k/2) < p/2. Since one can 
clearly choose k to satisfy the latter condition, the proof of the Lemma is complete. □ 



Lemma 9.3. Let fi(v) be defined by (9.9), where f(v) is given by (9.6), (9.2) and (7.3), 
and all involved multiple valued maps have their principal determination. Then there exists 
p > and a constant M, both independent of x e [0, 1], so that 

\fi(v)\<M, all\v\<p,xe[0,l]. 



Proof. By (9.26) and the proof of Lemma 9.1, there exist < S < 1 and p > 0, both 
independent of x, so that for any fixed x € [0, 1], all z's in the w-plane with \z\ < p can be 
written as z = v(t), with \t — 1| < S. Therefore, to prove the current Lemma it is enough 
to show that |/i(u)| = \fi(v(t))\ < M, for \t - 1| < 5 and x e [0, 1], for some constant M 
independent of x. 

By (9.9) and (9.8) we have 

^>=»- v2 {iH»- i/2 } 

= »-" 2 |^('-l)- 1/2 -^"- 1/2 }. 

where 
Write 



k p'(i) ; p'(iy/ 2 ' 
H(t) = {t-l) 1 ' 2 Q(t). (9.32) 



v = m t _f-\ t-\)=Tfm-l), (9.33) 
where £ is on the line segment [l,t]. Then (9.31) becomes 

AW _^_ 

= 1 r _ £0) \ , (go) _ go) 
(*- iy(o 1/2 lUw p'W / U'w p'(i) i /2p' ( i ) i/2 



+ 



p'(t) 1 /2 1 /(l)V2 p'^)l/2y( 1 )l/2 

' [H(t)-H{l)] + 4^L[U{t)-U{l)] 



)} 



p'(l)V2 



+1^5 TO -W]}. (9-34) 



where 



u ®=7m- (9 - 35) 
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We have 

H(t)-H(l) = (t-1)H'(£ 1 ) 
U(t)-U{l) = (t-l)U'fa) 
U(t)-U(0 = (t-OU'(r 1 ), 

where £1, £ 2 and 77 are on the line segment [l,t] (the latter because £ is on this segment). 
Therefore (9.34) implies 



\t-w(o i/2 \ i b'(*)i 1 vs; " 1 

\H(1)\ 



+ 



|p'(l)l/2| 



< L_ J 1 \„,(c y , 1^(1)1 ■ y | 1^(1)1 irr/^l 

(where at the last step we used that £ is on the line segment [l,t]). Therefore, to finish 

the proof we need upper bounds that are uniform in x E [0, 1] for and for |p'(i)| _1 , 
\H'(t)\, and \U'{t)\, for \t - 1| < 5. 
By (9.32) and (7.4), we have 

H(t) = t l {A-2xt) a /{A- xt)\ (9.36) 

where < a e Z, -1/2 < b e |Z and Z e ±Z. Since the right hand side of (9.36) is 
continuous in (x,t) on the compact set [0,1] x {\t — 1| < 6}, \H(t)\, and in particular 
|-ff(l)| can be majorized uniformly with respect to x € [0, 1]. 

Since by (9.35) U'(t) = -l/2\p'(t)]-^ 2 p"{t), all remaining uniform upper bounds will 
follow from such bounds for \H'(t)\, l/\p'(t)\ and \p"(t)\. 

By (9.36) we have 



l _ 1 (4-2xt) a l -2ax(4-2xt) a - 1 (4- xt) b + bx{4-2xt) a (4- xt)*- 1 

H ^ = lt (4-^ +t i^tr 



It follows from (7.3) that 



•"n—i-'Jxhp (9 - 38) 
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Since the right hand sides in (9.37) and (9.39) are continuous in (x, t) on the compact set 
[0, 1] x {\t — 1| < o*}, it follows that \H'(t)\ and \p" (t)| can be majorized as desired. Since 
\t - 1| < (5, and by (9.38) \p'(t)\ > q/\t\, one has l/\p'(t)\ < (1 + <5)/g, and the proof is 
complete. □ 

10. The asymptotics of a multidimensional Laplace integral 

By Proposition 8.1, the asymptotics of the moments M aii| g lv .. iami| g m;7l; 5 1; ... i7n; 5 n is given 
by the multiple integral in (8.3). In determining the asymptotics of the latter we will use 
the following result. 

PROPOSITION 10.1. Let h,a : D — > C, (0,1]™ C D C R n , be two functions. Assume that 
J[o i]n N < 00 an d that there exists a neighborhood V of (1, . . . , 1) so that h G C^V) and 
a e C 2 (V). Then for fixed q\,...,q n > we have 

I(R)= [ ■■■ f 1 x? qi ■■■x^h(x 1 ,...,x n y Ra(x »-^dx 1 ---dx n 

(10.1) 



1 h(l l) e iRa(l,...,l) 

iJ n [«i + iox 1 (l,...,l)]-"[gn + *Ox n (l,..-,l)] 



where a Xi is the partial derivative of a with respect to the variable Xi . 

We deduce the above result from the following. 

Lemma 10.2. Let k e (0,1) and let V C R n be an open set containing [k,l] n . Let 
h, a : V — > C be functions so that h e C 1 (F) and a <G C 2 (V). Then for fixed q\, . . . ,q n > 
we have 

h(R) = f - f ^ ■ ■■^ q "h(x u . . .,X n y Ra{x ^-^dx 1 ■ --dXn 

= i f- /^•■;^(M 2 ..., ;B ) e „ (1 „ 2 , )fc2 „, &n 

RJk Jk qi+ta Xl (l,x 2 ,...,x n ) 

+ 0(R- n - 1 ). (10.2) 

Proof. Write I k (R) as 

h(R)= f - f h{x ll ... 1 x n )e Rb ^-^dx 1 ---dx n , (10.3) 
Jk Jk 

where 

6(xi, . . . ,x n ) = qi hixi H h g„ lnx„ + ia(a:i, . . . , x„). (10-4) 

Apply integration by parts 8 with respect to the variable x\ in (10.3). We have 

d 



_^_ e Rb(xi,...,x„) = R (JL b { Xll . . . , In )j e «*(«.-.-»). 



(10.5) 



8 The idea of using integration by parts was suggested to the author by Andras Vasy. 
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By (10.4), 

-— b(x u ...,x n ) = — +i-—a(x 1 ,...,x n ), (10.6) 
oxi xi ox i 

so in particular b Xl {x\, . . . , x n ) ^ for all (xi, . . . , x n ) e [k, 1]™. Therefore (10.5) can be 
rewritten (omitting the arguments for brevity) as 

e m = 1 dm (10?) 
Rb Xl d Xl v ' 

Writing the second factor of the integrand in (10.3) as in (10.7) and applying integration 
by parts with respect to the variable x\ we obtain 

f 1 f 1 ( f 1 1 de Rb ) 



f 

Jk 



k Jk I Ri>xi 



k 



e Rb_^_ j | dxi dx 2 ■ ■ ■ dx n 



dxi \Rb xl 



= i r r m i ,-, >•,.. ,,.,, : .„ v , ll , 1 ... lh . i 

RJk Jk b Xl {l,x 2 ,...,x n ) 

_I f 1 ... f 1 k R ^x Rf » ■ ■ ■ X2 ' ' ' ' ' Xn \ e iRa ^-^) dx2 ... dXn 

RJk Jk b Xl {k,x 2 ,...,x n ) 

_ i r 1 r 1 h ^- f^ e H bdxi ... dXn (10 . 8) 

RJk Jk ip Xl f 

(where at the last equality we used (10.4) for the second term on the right hand side). 
By (10.6), 

\b Xl {x u ...,x n )\ > — > q u (x u ...,x n )€[k,l] n . (10.9) 
xi 

Therefore, the absolute value of the second term on the right hand side of (10.8) is ma- 
jorized by 



— k Rqi f 1 ■■■ f 1 
RQi Jk Jk 



\h(k, x 2 , ■ ■ ■ , x n )\dx 2 ■ ■ ■ dx ri 



which in turn, by the presence of the exponential and since h is continuous on [fc, 1]™, is 
less than Mii? - ™ -1 for R > n, for some suitable constants Mi and r\. 

On the other hand, the absolute value of the third term on the right hand side of (10.8) 
can be majorized, according to (10.9), by 

R ( qi y J "J \ h *i b xi ~ hb XlXl e Rb \dx 1 ■ ■ ■ dx n 

777 v? / x i ' x n 9 ™ l^i ^xi hb xixi \ dxi ■ ■ ■ dx n . 

R{qi) Jk Jk 
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Since h Xl b Xl —hb XlXl is continuous on [k, 1]" and Jj fe 1 , n irf 91 - ■ ■ x Rqn < J, Q 1 j„ x Rqi - ■ ■ x Rqn = 
n™=i(^« + there exist constants M2 and r<i so that the right hand side of the above 
equation is majorized by MiR~ n ~ x , for all R>r2- 

By (10.8), (10.6) and the last two paragraphs we obtain (10.3). □ 

Proof of Proposition 10.1. Choose k E (0, 1) so that [k, 1]™ C V. We note first that in 
order to prove the Proposition it suffices to show that Ik (R) has the asymptotics given by 
the right hand side of (10.1), where I k (R) is defined by (10.2). Indeed, I(R) = I k {R) + 
J(R), where 

J(R)= I x Rqi ■■■x^h{x u ...,x n y Ra{xi '-^ ) dx 1 ---dx n . 

•/[0,l]™\[fc,l] n 

Since throughout the integration range above one has Xi < k for at least one index i, 
denoting q = min(gi, . . . , q n ) > we obtain 

\J(R)\< [ k Rqo \h(x 1 ,...,x n )\dx 1 ---dx n 

V[0,l]"\[fe,l]" 

< k Rqo / \h(xi,...,x n )\dx\ ■ ■ ■ dx n - 
J{o,i]« 

By hypothesis the integral above is finite, so there exist constants M and r so that 
I J(R)\ < M i?(7" _1 for all R> r . Therefore addition of J(R) to h{R) does not change 
its asymptotics. 

By Lemma 10.2, the asymptotics of Ik(R) is the same as 1/R times the asymptotics 
of the integral on the right hand side of (10.2). In turn, the latter integral meets the 
hypotheses of Lemma 10.2. Applying Lemma 10.2 to it we obtain from (10.2) that 

f 1 x Rqs ■■■x R ^h(l,l,x 3 ...,x n ) 

J k [qi +ia Xl (l,l,x 3 , . . .,x n )][q 2 +ia X2 (l, l,x 3 , . . .,x n )] 
x e iRo,{i,i,x 3 ,...,x n ) dXi ... dXn + 0{R- n - x ). 

Lemma 10.2 can be applied again to the integral on the right hand side above, yielding 
an (n — 3)-fold integral that meets its hypotheses. All the successive applications in this 
manner of Lemma 10.2 lead to multiple integrals to which it is again applicable. After n 
applications one obtains (10.1). □ 

By (8.1) and (8.2), the multiple integral (8.3) leads to integrals having the form of I(R), 
but with the exponential function replaced by a cosine. As indicated below, this type of 
integrals can easily be handled by Proposition 10.2. 
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Corollary 10.3. Let h,a,c: D R, (O^fcDcl", be three functions. Assume that 
I[o i] n \h\ < 00 an d that there exists a neighborhood V of (1, . . . , 1) so that h,c G C (V) 
and a £ C 2 (V). Then for fixed q\, . . . , q n > we have 

K{R)= / •••/ %i qi - ■ -x^ n h{xu ■ ■ ■ ,x n ) cos[Ra(x 1 , . . . ,x n ) + c(x 1 , . . . ,x n )]dxx ■ ■ ■ dx n 
Jo Jo 

1 h(l, ...,1) 



R n ^{ qi f + (a xi (1, . . . , W ■ ■ ■ ^(q n y + (a x Jl, I)) 2 
x cos 



Ra(l, . . . , 1) + c(l, . . . , 1) — arctan Qxi — ' ' ' ' ^ — ... — arctan & ' 



91 9n 

+ o(ir™- 1 ), (10.10) 

w/iere <z Xi zs i/ie partial derivative of a with respect to the variable Xi . 

Proof. Since h, a and c are real- valued, cos(Ra + c) = Ree iiia+lc , and thus 

X"(fl)=Re /"■■■/ if'-xf'"^ i^'^e^* 1 '">(iii-(ii„. 

Jo Jo 

Therefore, by Proposition 10.1 

1 Ml l) e i[Ra(l,...,l)+c(l,...S)] 

K(R) = — Re T 1 = +0(R- n - 1 ). (10.11) 

V ! R n [ qi +ia xl (l,...,l)]---[q n +ia Xn (!,...,!)] V ; V ' 



Let aj = a x . ( 1 , . . . , 1 ) and write qj +iaj = \J Qj + a "j ( cos 0j + i snl Qj ) i with 6j = arctan , 
for j = 1, . . . , n. The contribution of the exponential and the denominator to the fraction 
in (10.11) becomes 

e i[fla(l,...,l)+c(l,...,l)-0i e n ] 



y/qj+aj ■ ■ ■ ^/qj + a? n 
Substituting this in (10.11) gives (10.10). □ 

11. The asymptotics of Wj,. Proof of Theorem 2.2 

In this section we use the results of Section 10 to deduce the asymptotics of u>b from its 
expression (6.5). To this end, we need to determine first the asymptotics of the moments 
M au u ..., am ,0 m . t7uSl ,..., 7n ,6 n defined by (6.6). 

Recall that the arguments of u>b — and therefore those of the moments M — approach 
infinity as specified by (2.3). In particular 9 , Vk = quRk + Ck, k = l,...,m and v'j — 
q'jR'j + Cj, j = 1, . . . ,n. 



9 In this section doubly indexed variables appear alongside the complex number i. For this reason, we 
use here k and j as opposed to the more familiar i and j as indices of these doubly indexed variables. 
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Substituting (8.1) and (8.2) into (8.3) we obtain that 



*)2m 



Mr, 



1 



io Jo JJ 



i 



2 , n"=i ^fc 3 yfcj 



yfe t 4-n" =1 Zfejjufcj 



><n 



n 4/ /2 , Yl-k =1 Xk 3 z~ 4/ ,2 , nr=i^fcj"'fcT" 



nr=i ^fcj^fcj / n fcLi vkjw kj 

4 ^lirli x kjZ kj V 4-H™ =1 y kj w kj 



*;=i V 

m ^ 

n 

*;=i 1 

X COS 



rife=i x kjZkj 



1i 



-Rah 



4 -n"=i^fej^- 

nm 



4 - nr=i ^ ■ / v J V 4 ~ nr=i ^ 



cos 



Rk arccos 1 



,•=1 arfcjZ/fcj \ 1 



1 / U^iXkjVkj a k TT 

I — - arctan — t / 1 — — 

y 2 g fc y 4-}{ j=1 x kj y kj 2 



><n 



cos 



Rh arccos 



Rj arccos 



, iYj=i z kjw k j\ i i / n"=i ^fe^fej /?fc7T 

1 — ^ - x arctan — \ / ^= 1 — — 



Ik V 4 - n?=i ^fej Wfei 2 



! _ m^MlM.) + I arctan 1 J + fa ~ 1} " 



X cos 



m n 



J2 - arccos 1 - llfc - x — i + - arctan —W - 1X *-^ — ^— + y 

' V 2 / 2 «5 V 4 - n fe =i ykjwkj 2 



x n ii 1 -'-./^^./"-^ 



(ii.i) 



One readily proves by induction on s that 



cos 0i • • • cos# s = — ^ cos(ei#i + • • • + e s s 

ei,...,€,=±l 



Since by (2.3) we have R k = A k R, fc = 1, . . . , m and i?^ = j = 1, . . . , n, the product 
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of the 2m + 2n cosines in the integrand of (11.1) becomes by the above formula 



22m+2n / / 

€l,...,e2T7i+2n = ±l 



COS • 



R 



( A (, Y\ n j=i x koVk 3 \ ( Y\ n 0=1 z k jw k j 

} / [ e 2 k-\A k arccos I 1 — I + € 2 kA k arccos I 1 J - — 



k=l 



+ £ (e^+v-iBj arccos (l - IIk|MfM^ + £2m+2 jBj arccos _ UZiV^k^ 



El £2fe-l , 1 / V^=\ x ki,Vk 3 e 2k 1 / n"=l z *:j w fcj 
— - — arctan — * / h — — arctan — i / 

fe=1 V 2 * V 4 ~ Uf=i ^fej'J/fcj 2 % V 4 - n j= i ^fejWfej 



El £2m+27-l 1 / Ili— 1 x kjZ k j , £2m+2j 1 / II*— 1 Vkj w ki 
^— arctan— W . ffm + — arctan — W 1 ^ — J — 



+ ^2(e2k-ia k + e 2k (3 k )- + ^(e 2m +2j-i7j + e2m+2j5j)- - ^(e 2TO +2j-i + e 2m +2j)-^ \ ■ 

k=l j=i j=i 



(11.2) 



The 4mn-fold integral in (11.1), with the product of cosines in the integrand written 
as in (11.2), becomes the sum of 2m + 2n integrals of the form of the multiple integral in 
Corollary 10.3. Our current functions h satisfy J, Q 1 j 4m „ \h\ < oo. Indeed, since the product 

of the 2n products (8.12) is Ylf^iYl^i XkjVkjZkjWkj , one readily sees that \h\ can be 
majorized by KY\^ =1 Y\^ =1 {x k jy k jZ k jW k j) Ck+c i +1 ^ 2 , where if is a constant independent 
of x k j,y k j, z k j,Wkj- Since by (2.3) the exponent in the previous double product is non- 
negative, it follows that Jj ^ 4mn \h\ < oo. 

The other conditions in the hypothesis of Corollary 10.3 are clearly satisfied. The a- 
functions of Corollary 10.3 are now the coefficients a(xu, . . . , w mn ) of R in the argument 
of the cosines in (11.2). 

Note that there are exactly two terms of each of our a(xn, . . . , uw)'s containing any 
given x k j, y k j, z k j or w k j, and that all variables are set to 1 on the right hand side of 
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(10.10). Since d/dx arccos(l — xu/2) — (Ax/u — x 2 ) 1 / 2 , we obtain 



da 


(1,- 


..,1) = 




dx k j 




da 


(1,- 


..,1) = 


t2k-\Ak + (2m+2jBj 


dykj 




da 


(1,- 


..,1) = 


£2kAk + £2m+2j-lBj 


dz k j 




da 
dwkj 


(1,- 


..,!) = 


(■2k Ak + f.2m+2]B J 

V3 



(11.3) 



By (11.1)-(11.3), Corollary 10.3 implies 



J«ai,/3i,->am,/3m;7i>«i.->7n,«n - ^4mn 11 ^ ^3 J 11 \ / 



E 



■ cos 



£)€i,. ..,e 2m +2» I 3 

ei,...,e2m+2n = ±l 



fe=l j=l 



^ Ak(e2k-1 + t2k) + Bj(£2m+2j-l + £2m+2j) 
fc=l j=l 

m n 1 

E£2fe-l+e2fe J- . e 2m+2j-l + £2m+2j J- 
arctan = + > arctan = 
fe=1 2 »V5 ^ 2 

m ^_ n ^ n ^ 

+ ^(e2fe-iafe + £2fc/3fc)-^ + E^ 2m + 2 J- 17 -?' + e 2m+2j5j)- - y^(e 2TO+ 2j-i + e 2m +2j)- 

fe=l j=l i=l 

I £2k-lAk + C2m+2j-lBj C2k-lAk + €2m+2jBj 

- > > arctan ^ + arctan 

t^ij^A (ikAk + q'jB^VS (q k A k + q'B^VS 

£2kAk + €2m+2j-lBj ^2kA k + £2m+2jBj \ I 

+ arctan — + arctan —pr > 

(q k A k + q'jB^VZ (qkA k + ^B^y/Z J J 

+ 0(i?~ 4m ™+S™=i( Qfc + /3fc )+S" = i(7j+5:i)-3m-3ri-l^ (H-4) 



where 



E 3- 1 rWgg 
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and 



k=ij= 



(q k A k + q'jBj) 2 + 


{tl k -\A k + €2 m +2j-lBj) 


3 


1 {q k A k + q'jBj) 2 + 


{t2k-iA k + e 2m +2jBj) 2 


3 


l( qk A k + q rB 3 ) 2 + 


(e2 k A k + C2m+2j-lBj) 2 


3 



x v(g^+^-) 2 + (c2fci4 * +c r +!wBj)a 



(11.6) 



Define 



rfei,...,€2m+2n 

ai,/3i,...,a m ,/3 m ;7i,5i,...,7 n ,(5 n ' 

TT — = ) TT ( ~P J x (summand of the (2m + 2n)-fold sum in (11.4)). 

fe=i VV3/ fJi\V3j (117) 

We call (ei, . . . , e 2 m+2n) balanced if £2j-i + £2j = 0, j = 1, . . . , 2m + In. 

Lemma 11.1. Fix ej € {—1, 1}, j = 1, • • • , 2m + 2n. ^4s m Section 6, /eiC 6e £/ie collection 
of terms obtained by expanding out the left hand side of (6.2), and for any C e C let e(C), 
a k (C), (3 k {C), 7j(C) anrf Sj[C) have the same significance as in (6.2). 
T/ien unless (ei, . . . , £2m+2n) *s balanced, one has 

e ( C )^aV(C)A^),...,a m (C),/3 m (C); 7 i(C),5 1 (C),...^(C'),5„(C) = °" ( 1L8 ) 

cec 



Proof. Suppose there exists j € [l,m] with £2j-i = e 2j- Without loss of generality we 
may assume ei = £2- Then it follows from (11.7) that 

C e lv'! e 2m+2n c^l , ■ ■ ■ ,^2m + 2n ^ -| 

ai,/3i,a2 ) /32v) a m)/3m;7i)5i,...,7 n ,5 n ai — l,pi + l,a2,02,---,& m ,0 Tn ;yi,5i,...,'y n ,d n ' 1 — 

(11.9) 



Partition the terms in the expansion C of the left hand side of (6.2) into four classes C\, 
C2, C3, C4, according to which of the terms of the factor (a\ — bi) 2 = (a 2 — a\b\ — b\a\ +b\) 
(in order, from left to right) is chosen when expanding. 

By (11.9), the restriction of the sum (11.8) to C\ is canceled by its restriction to C2, 
and the restriction to C 3 is canceled by the restriction to C4. This proves (11.8). 

The case when e 2 j-i = e2j, j € [m + l,m + n], is treated in a perfectly analogous 
manner. □ 
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Define Ma 1 ,p 1 ,a2,P2,...,a m ,p m ;n,s 1 ,...,-y n ,s n to be § iven b y tlie expression on the right hand 
side of (11.4), when the summation range is restricted to balanced (ei, . . . , £2m+2n)'s. 
Clearly, for (ei, . . . , £2 TO +2n) balanced we have 

{t2k-lAk + £2m+2j-lBj, 62k-lAk + €2m+2jBj, 62fc^4fc + £2m+2j-lBj , €2kAk + £2m+2jBj} 

= {Ak + Bj,Ak — Bj, —Au + Bj, —Au — Bj}, 

for all k = 1, . . . , m, j = 1, . . . , n. It follows that for (e 1; . . . , e 2m +2n) balanced one has 

(i) the double sum of arctangents in (11.4) is 0, and 

(ii) the double product fl fl '-> f2 »+ 2 » in (11.6) equals 



fe=i j=i 



( (Zfe A fc + g ;-B J ) 2 + -( J 4 fe -B,) J 



1 



(g fc A fc +^B,) 2 + -(A fe +B, 



an expression independent of (ei, . . . , £2m+2n). 

By (i) and (ii) above we obtain from our definition of the Mq's that 



(11.10) 



M° 



E 



ai,/9l,a2,/32,.") a m)/3m;71i5l,...,7n,5n 

X ^ COS 

e2i-i=±l, (=l,...,ro+n 



n 



Ri 



=i 



afc+/9fc ™ 



n 



/ \ 7j+<5j 



2- 5 n + 2^ 



.7=1 x v ' 

e2m+2j-l(7j 



i=i 



(11.11) 



By (6.2), when our parameters depend on R as in (2.3), the order in R of e(C) plus 
SfeLi( Q; fe(C) + Pk{C)) + Sj=i(7j(C) + ^j(C)) i s always less or equal than the number of 
factors of the product £ of (6.2), namely 2m+2n+4(™)+4(") = 2m 2 +2n 2 . Thus, the omit- 
ted part of the approximation of e(C , )M Ql(c . ) ^ l(c)i ... jQm(c) ^ m(c);7l(c)i(5l(c) ^.. j7ri(c)i( 5 ri(c) 
resulting from (11.4) has at most order 2m 2 + 2n 2 — Aran — 3m — 'in — 1 in R. Therefore, 
by Lemma 11.1 and our definition of the Mq's, (6.5) implies 



Ri 



Rm . R'l R' n \ _ 



X 



i(C),Pi(C),...,a m (C),P„ 



(C); 7 i(C)A(C),..., 7 „(C),5„(C) 



cec 
+ 0(R 



2rn 2 +2n 2 -Amn-2m-l\ 



where 



X = X 2m,2n [J R k ]J ~ 1/2) (^ + 1/2), 

k=l j=l 

with x given by (4.3). 



(11.12) 
(11.13) 



Lemma 11.2. We have 



M° 

I "ai,/3i,...,a m ,/3 m ;7i,<5i,...,7„,<5„ 



a k m 

n 



-iRi 



0k ™ 



iR 



, , \ vo / , i V v3 / i V vo / , \ v3 

fc=i x ' fe=i x ' j=i x ' j=i 

if ai = ft (mod 2), . . . ,7„ = <5„(mod2) 
otherwise, 



where 



Y 



)m-non 



n; =1 + 1 



(11.14) 



nr=i ny=i[(9*^ + + - B 3 m qk A k + ^ + + b y\ (11 15) 

Proof. If ai = ft (mod 2), . . . ,7„ = 5„(mod2), all the constants multiplying 7r in the 
argument of the cosine in (11.11) are integers. Since cos(— nit + x) = cos(n7r + x) for all 
n € 1 and all ir, it follows that all 2 m+ ™ terms in the sum (11.11) are equal to the term 
corresponding to (ei, e 3 . . . , e 2m +2n-i) = (1, . . . , 1), which is 



cos 



V 



ak - ft 



^ 2 

. fe=i 



7j 



-7T = (-i)S. 



= JJi a *(-i) ft JJ*^(-t)^. 
fe=l j = l 

Using this and substituting in (11.11) the expressions for E and D given by (11.5) and 
(11.10) we obtain (11.14). 

Assume next that there is an index k so that a k and (3k have opposite parity. Since 
cos(— nir + x) = — cos(nn + x) for all n G \ +Z and all x, it follows that the €2k-i — 1 part 
of the multiple sum in (11.11) cancels the C2k-i = — 1 part, so the multiple sum is 0. The 
same argument applies if there is an index j so that -fj and Sj have opposite parity. □ 



By (11.12) and Lemma 11.2 we obtain 
(Rx Rm R' } ...K\ =XY 



L e{c) \R^) Rv^-) 33 33 

+ o(R 2m2+2n2 - 4mn - 2m - 1 ), (11.16) 
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where for any monomial /i = IlfeLi x k k y k k IljLi z ] 3 w j 3 > the angular brackets denote (/i) = 
/i in case a k = (3 k (mod2), 7^ = 8 j (mod 2), for all k and j, and (/i) = otherwise. 

Note that if in (11.16) we had (/z) = it for aZZ monomials /x, the sum in (11.16) would 
evaluate to the specialization of the expression £ of (6.2) when 

ak= vf' 6fc = ~7f' Cj= vT^vT (1L17) 

fc = 1, . . . ,m, j = 1, . . . ,n. 

This observation, together with the fortunate situation presented in the result below, 
allows us to evaluate the sum (11.16). 

Lemma 11.3. Under the substitutions (11.17), all contributions to the expansion of £ com- 
ing from monomials [ic = e(C)a 1 1 ^6f 1 ^ ■ ■ ■ a^^b^^cl 1 ^^ 1 ^ ■ ■ ■ cZ"^ dfr^ C \ 
C GC, for which not all pairs (a k (C), f3k(C)), (7j(C), 5j(C)) have components of the same 
parity, cancel out. 

Proof. We define first the following map C C' on the subcollection Co consisting of 
those C G C for which not all pairs (afc(C), f3k{C)), (7j(C), 5j(C)) have components of the 
same parity. 

Totally order the disjoint union of the set of indices in the pairs (otk(C), f3k(C)) with 
the set of indices in the pairs (7j(C), Sj(C)). Let C e Co, and consider the smallest pair 
index for which the components have opposite parity. Assume that this smallest index k 
occurs in a pair of the first type, (a/ So (C),/3fe (C)). Our arguments apply the same way to 
the case when the smallest such index occurs in a pair of type (-fj(C), Sj(C j). 

By (6.2), C is obtained by selecting a signed term from each of the 2m+2n+4(™ ) +4(") 
factors of £. Define C as being obtained by making the following selections: 

(i) From the m pairs of factors (ak — bk)(ak — bk), k = 1, . . . , m, make the selection as 
follows: 



Hows: 

(a) if / ^ ko, select in C the same signed terms as in C; 

(b) if k = ko, select in C the signed term of each factor of 



' (a ko - b k „)(ak - b ko ) that 

was not selected in C . 

(ii) From the n pairs of factors (cj — dj)(cj — dj), j — 1, . . . , n, select in C the same signed 
terms as in C. 

(Hi) From the (™) groups of four factors 

((vj - v k ) + a,j - a k )((vj - v k ) + aj - b k )((vj - v k ) + bj - a k )((vj - v k ) + bj - b k ), 

1 < k < j < m, make the selection as follows: 

(a) if ko {k,j}, make in C the same selection as in C; 

(b) if k = ko, and the selected terms in C from the four factors 



(("j - v ko ) + aj - a ko )((vj -v ko ) + aj - b ko )(( Vj - v ko ) + bj - a ko )((vj - v ka ) + bj - b ka ) 
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are the fcith, k2th, k^th and feith, respectively, then select in C the feth, fcith, A^th and 
k^th terms of the above factors, respectively; 

(c) if j — fc , and the selected terms in C from the four factors 

(Ko - V k ) + Ofeo - afcXKo ^ w fe) + fl feo - frfcXKo - w fe) + b ko - d k )((v ko - Wfe) + frfco - &fc) 

are the kith, fc 2 th, fc 3 th and fctth, respectively, then select in C the ^th, A^th, fcith and 
feth terms of the above factors, respectively. 

(iv) From the final (™) groups of four factors make the same selections in C as in C. 
Let us now compare the monomials 

m m 

Mc = e (c)nc (c) ^ (c) ri c ? (c) 4' (c) 
fe=i j=i 

and 

^ = e(C") J] a- (C,) 6^ (C,) J] cf^df^ 

k=l j=l 

generated by the selections C and C . 

It is clear from our construction that the portions of C and C covered by step (Hi) 
produce contributions to /ic and fie whose only difference is that the roles of a ko and b ko 
are interchanged. 

The same is true for step (i). Indeed, for k = ko, to the four possible selections 
0000, ( a k„)(-b ko ), {-h )(a ko ), (-b ko )(-bk„) in C there correspond the four selec- 
tions {-b ko )(-b ko ), (-b ko )(a ko ), (a ko )(-b ko ), (a ko )(a ko ) in C , respectively So the four 
possible contributions to fi c and nc are a\\P ka , -a ko b ko , -a 2 k b ka , a° ko b 2 ko and »° 6f , 
-ak b ko , -a 2 ka b ko , a 2 ko b° ko , respectively 

Since steps (ii) and (iv) do not involve the index fc , and the portions of C and C 
covered by these steps coincide, it follows from the previous two paragraphs that 

? ®ko 5 b k Q , . . ■ , dm ? b m , c\ , d\ , . . . , c n , d n ) 

= Hc(ai,bi, . . . ,b ko ,a ko , . . . , a m , b m , c\,d\, . . . , c n ,d n ). 

Since by (11.17) a ko and b ko become the negatives of each other, and since a ko (C) and 
(3 ko (C) have opposite parity, it follows that under the specialization (11.17) the monomials 
lie and [ic cancel out. 

Our map C ^ C clearly maps Co to itself, and sends C back to C. By the previous 
paragraph, when making the specialization (11.17), all the monomials corresponding to 
term selections in C cancel out in pairs. This proves the Lemma. □ 
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Corollary 11.4. We have 



L0 b 



Ri R m i?j 



^™ 1 - XY 



V3 



n 



-iRk 
V3 



Ac(C) n /,-D/N7i(C) n 

S(7t) 11 



/ \ <5 3 (C) 



2m 2 +2n 2 -4mn-2m-l 



(11.18) 



where X and Y are given by (11.13) and (11.15). 

Proof. This follows directly from (11.16) and Lemma 11.3. □ 

Our results allow us now to obtain the asymptotics of thus proving Theorem 2.2. 

Proof of Theorem 2.2. By Corollary 11.4, it suffices to evaluate the sum in (11.18). 
However, by (6.2), this sum is just the specialization (11.17) of the product £ in (6.2). 
The product of the first 2m + 2n factors of £ specializes to 



mrm 



I A\m+n m n 



(11.19) 



fc=i j=i 



The product of the next 4(™) factors of £ specializes to 



n 

l<k<j<m 



x (vj -v k )+[- 



- n 

l<k<j<m 





iRk 


V3 


V3 


iRj 


iRk 


V3 


V3 



("j - «*) + 7= + 



iRj iRk 



(vj - v k ) 2 + -(Rj - Rtf 



{v 3 -v k Y + -(R j + R k ) 



V3 V3 

2 



(11.20) 



while the product of the last 4(") factors of £ specializes to 



n 

l<fe<j<n 



VV3 V3 
V3 



n 

l<k<j<n 



1 



K-o 2 + x(i?;- + i?D 2 



(11.21) 
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In the expression Y given by (11.15) write, using (2.3), 
T 1 



3 R k 



1 



1 



(QkRk) 2 + ^Rl = Wk \l(v k -c k ) 2 + -Rl 



1 1 



1 



and 



(q k A k + q' j B j ) 2 + -(A k ±B j y 



1 

i? 2 

1 

i? 2 



1?' 



{q k A k R + q'.BjR) 2 + -(A k R±B 1 R) 2 



(v k +v' j -c k -c' j ) 2 + -{R k ±R' j f 



Substituting the resulting expression for Y and formula (11.13) for X into (11.18), we 
obtain by Corollary 11.4 and (11.19)-(11.21) that 



i?i R m R[ R' 
u>b • ' ' ) i • • • 



im-noTi 



x 2m , 2n n ^ n #;(^ - 1/2) (i?; + 1/2) 

fe=l j=l 

n^ lv /K-4) 2 + i(^) 2 



nr=i n;=jK + ^ - <* - ^ + i(R k - ^hk + ^ - <* - c ;.) 2 + + ^) 2 i 



(-4) 



m+n 



m n 



>< n 

l<k<j<m 



fe=i i=i 



(«j -Wfe)^^^-^) 2 



x n 

l<fe<j<n 

->2m 2 +2n 2 -4mn-2m-l 



1 



K-^) 2 +3(^-fife) 2 



( t , J .- Wfc )2 + _(i2 J .+iJ fc )S 



K-<) 2 + 3(^+fife) 2 



+ 0(R 2 



(11.22) 



By (2.3), the parameters R\, . . . , R m , v\, . . . , v m and R[, . . . , R' n , v[, . . . , v' n approach in- 
finity as R — > oo, while c\, . . . , c m and c[, . . . ,c' n arc constant. It follows from this that 
the difference between the product on the right hand side of (11.22) and what it becomes 
when one omits the additive constants ±1/2 on its first line and the constants c k and c'j 

is o(i? 2m2+2 " 2 - 4 " m - 2m - 1 ). The latter is readily brought to the form given by (2.9) and 
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Figure 12.1. E 2 [l, 2, 3, 4; 1,2,3,4,5,6]. Figure 12.2. Region denoted 

#[2,3,4,5,6,7] [2,3,4,5] ( 6 ) in [5]. 



(2.10). By the assumption on the distinctness of the pairs (Ak,qk), k = l,...,m, and 
(Bj,q'j), j = 1, . . . ,n in (2.3), the expression on the right hand side of (2.9) has degree 
2m 2 + 2n 2 — Amn — 2m in R. Therefore, (2.9) does indeed give the asymptotics of uib- □ 

12. Another simple product formula for correlations along the boundary 

By Proposition 3.2, the asymptotics of the joint correlation at the center oj will follow 
provided we also work out the asymptotics of the boundary-influenced correlation u>b 
defined by (3.3). We present this in Section 13. But first we need an analog of Proposition 
4.1 corresponding to the regions En used to define Qf,. 

Let E be the region determined by the common outside boundary of the regions 

En ( ^ 1 • • • ^ m \^} ■ ■ ■ ) defined in Section 3, for fixed m, n and N. Then E is 

the half-hexagonal lattice region with four straight sides — the southern side of length 
N + 2n + 1, southeastern of length 2iV + Am — 1, northeastern of length 2N + An, and 
northern of length N + 2m — followed by N + 2n descending zig-zags to the lattice point 
O, two extra unit steps southeast of O followed by one step west, and N + 2m — 1 more 
descending zig-zags to close up the boundary (an example can be seen in Figure 12.1). In 
addition, the N + 2m— 1 dimcr positions weighted 1/2 in the regions En are also weighted 
so in E. 
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As in the case of the regions W of Section 3, the vertical, jagged boundary of E can 
be viewed as consisting of bumps — in the present case, pairs of adjacent lattice segments 
forming an angle that opens to the east: N + 2m — 1 bumps below O, and N + 2n above 
O. Label the former by 0, 1, . . . , N + 2m - 2 and the latter by 0, 1, . . . , N + 2n - 1, both 
labelings starting with the bumps closest to O and then moving successively outwards. 

In the description of E the parameters m and n always appear with even coefficients. 
We re-denote, as in Section 3, 2m by m and 2n by n, for notational simplicity. Therefore 
we consider the four straight sides of E to have lengths N + n, 2N + 2m — 1, 2N + 2n 
and N + to, while the number of bumps below and above O is N + to — 1 and N + n, 
respectively. 

We allow removal of any bump exactly like in Section 3: above O, place an up-pointing 
quadromer across it and discard the three monomers of E it covers; below O, use down- 
pointing quadromers. 

Define En\k\, . . . , k m ; h,. . . , l n ] to be the region obtained from E by removing the 
bumps below O with labels 0<fci<fc 2 <---<fc m <A r + TO — 2, and the bumps above O 
withlabelsO< h < h < ■ ■ ■ < l n < N+n-l. Figure 12.1 shows E 2 [l, 2, 3, 4; 1, 2, 3, 4, 5, 6]. 

The correlations of the removed bumps on the boundary of the regions En turn out 
to be exactly the same as the corresponding correlations we found in Section 3 for the 
regions Wm- 

Proposition 12.1. For m, n > and fixed integers < fci < fc 2 < • • • < k m and 
< h < h < ■ ■ ■ < In we have 

M(E N [k 1 ,...,k m ;h,...,l n }) 

hm — : 7- = 

N^ooM(E N [0, . ..,to- 1;0, . ..,n- 1]) 

n n Vs-ii) 

n (3/2)fc i T-r (3/2)^ l<Kj<m l<i<j<n 

' {)ki ^ {)h f[f[ {ki +l j + 2 ) (12-1) 



where Xm,n is given by (4.3). 



Proof. We proceed in analogy to the proof of Proposition 4.1. The results of [5] provide 
an explicit formula for M(£jv[fci, ■ ■ ■ ,k m ;h, . . . ,l n ]) as well. Indeed, in the notation of [5] 
we have 

E]y[ki, . . . , k m ; h, ■ ■ ■ ,l n ] = 

^[l,...,JV+n]\[/ 1 + l,...,; n + l],[l,...,Ar-l+m]\[/si+l,...,fe m + l](^ + m ) (12-2) 

(this is illustrated by Figure 12.2, which is just the 180° rotation of Figure 12.1). 

Proposition 2.1 of [5] and formulas [5, (1.1), (1.3), (1.5)] provide an explicit formula for 
M(i? Pj q(a;)), for any pair of lists p = [pi, . . . ,p 8 ], 1 < pi < ■ ■ ■ < p s and q = [q lt q t ], 
1 < qi < ■ ■ ■ < qt, and any nonnegative integer x > qt — p s — t + s — 1. 

75 



Written as a constant times a monic polynomial in x, this formula takes the form 



M(i2 p , q (a;))=Cp, q Gp, q (a:) ) (12.3) 

where by [5, (1.3)] we have 

r -oCVh'TJ 1 TT 1 ni<i<j<sfa-PWi<i<i<t(*i-*) n9d) 

p,q l\(^y-l\(^-^ nunU(Pi+^) " ( } 

Furthermore, it is straightforward to check, using (1.1) and (1.5), that the polynomials 
Gp, q (x) satisfy 



G p io iq (ar) 



(x-pi+p s )(x+pi+p s - s + t + 2), for 1 < i < s (12.5) 



i 

(x + Qi +p s + i)(x - qi +p s - s + t + 1), for 1 < i < t 



(as in Section 4, pi') denotes the list obtained from p by increasing its i-th element by 1, 
and is defined only if — li > 2). 
Consider first the limit 

lim M (E N [ki, kj-uh + 1, kj+u . . . ,fc m ;/i, . . .,?„]) 7 ^ 
w^oo M (E N [ki, . . . , ki-i, hi, k i+ i, . . . , k m ;h, . . . ,/„]) 

for ki + 1 < fcj+i. 

Use (12.2) to view the regions involved in this fraction as i? Pjq (x)'s. Comparing formulas 
(12.4) and (4.6), and using ( s ~*) - t= ( t ^ s ) - s, one sees that 

c p,q = c q,p- 

For the lists on the right hand side of (12.2) this implies 

C[l,...,iV+n]\[i 1 +l,...,; n +l],[l,...,JV-l+m]\[fei+l,...,fe m +l] 

= C[i,...,iV-l+m]\[fci+l,...,fc ro +l],[l,...,JV+n]\[Ji+l, ...,«„+!]• ( 12 - 8 ) 

Therefore, the contribution to the fraction in (12.7) coming (via (12.2)) from the c Piq -parts 
of (12.3) follows from our work in Section 4. The only difference from that case is that 
in the c Piq of (12.8) the largest element of p is now (N — 1) + m, as opposed to N + m. 
However, substitution of N by N — 1 in (4.10) makes no difference in the limit N — > oo, 
so the contribution of the c Pjq -parts to the limit (12.7) is, asymptotically for N — > oo, 
precisely the same as (4.10). 
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On the other hand, one readily sees by (12.2) and (12.6) that the contribution of the 
G Pj q(x)-parts to the fraction in (12.7) is 



1 

(2N + m + n + ki + 2) (2N + m + n - k % - 1) ' 

a fraction whose asymptotics as N — > oo is clearly the same as that of (4.11). 

It follows from the above two paragraphs that the limit (12.7) is equal to the product 
on the right hand side of (4.12), so it has exactly the same expression as in the case of the 
regions Wjy treated in Section 4. 

A similar analysis shows that 

M (E N [ki, . . . , k m ; h, ■ ■ ■ , h-i,lj + 1, h+i, . . . ,l n }) 
n^oo M (£?jv[fci, . . . , k m ; h,..., k-i,k, h+i, . . . , l n ]) 

M (Wjy[fci, . . . ,fc m ;li, ■ ■ .,h-i,h + l,h+i, ■ ■ ■ ,l n \) 
jv^oo M (Wjv[fci, . . . , k m ; h,..., k-i,h, k+i, ■ ■ ■ , l n ]) 

so that also decrementing a single element from the second index list has exactly the 
same effect for the regions E N as for the regions Wn of Section 4. Formula (12.1) follows 
therefore by Proposition 4.1. □ 

13. The asymptotics of Cb b . Proof of Theorem 2.1 

Our reasoning from Section 5 applies with no change to the regions Em defined in 
Section 3. Each missing quadromcr can be placed back at the expense of performing a 
Laplace expansion in the corresponding Gessel-Viennot determinant. The resulting 2x2 
cofactors have precisely the forms (5.5) and (5.13). The analog of (5.19) we get is 

m ( En ( vi ■■■ R vZ ; = 2m+ " n r ^ - i/2m + i/2) n ^ 



... (-i)E™i( a '+^)+sr=i( c »+*) 

<di<i?; 

- ai)(Ri + <n - l)\(Ri + bi- 1)! 



0<ai<6i<i?i 0<a m <b m <R m 0<c 1 <d- L <R' 1 0<c n <d n <R 



X J] (2at + i)! _ a .)i (26- + l)! (jj. - 6,)! 
tt & - Ci)(J^ + Cj - 1)! + - 1)! 

x sgn(v[ + a,v' 1 +di,...,v' n + c n ,v' n + d„)sgn(«i + ai,vi + bi,...,v m + a m ,v m + b m ) 
x M(E N [{ Vl + ai ,...,v m + 6 m }<; {v[ + Cl ,...,v' n + d„}<])| , (13.1) 

where the summation extends only over those summation indices for which both lists of 
arguments in the region on the right hand side have distinct elements, and sgn denotes 
the sign of such a list when regarded as a permutation. 
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After removing the forced dimcrs, the normalizing region at the denominator of (3.3), 



,13 2m - 1 1 3 In - 1 , , > 

En \o ()••• o -0 ()••• o )- (L, -> 

is readily seen to be precisely Ejy[0, 1, . . . , 2m — 1; 0, 1, . . . , 2n — 1]. 

Moreover, by the sameness of the right hand sides of (12.1) and (4.2), (5.21) implies 

,. M (EnHvx + ai, . . .,v m + 6 m }<; {v[ +c 1 ,...,v' n + d„}<]) _ 
nZo M(£ JV [0,...,2m-l;0,...,2n-l]) ~ X2m ' 2 " 

x sgn(t>i + oi, ui +6i, . . .,« m + a ro ,u m + b m )sgn(v[ + d,v[ + d 1 , . . .,v' n + c n ,v' n + d n ) 

fr (3/2)„ i+ai (3/2)„ i+6i A (3/2)„, +Ci (3/2)„, +di 
X fi (2) B4+04 (2)„ i+i4 

x [ | (vj -Vi + aj - ai)(vj -Vi + aj - bi)(vj -Vi + bj - di)(vj -Vi + bj - h) 

l<i<j<m 

X 1 [ ( V 'j ~ V 'i + °i ~ C »)K' ~ V 'i + °i ~ di ^ V 'j ~ V 'i + d 3 ~ C »)K' ~ V 'i + d - d i) 



l<i<j<n 

UZl( a i - b i)U7=l( c i ~ d i) 

U4L1 n?=i( u y + a i + c j)( u ij + a i + dj)(v,ij +h + Cj)(uij +h + d 3 ) ' 



(13.3) 



where x is given by (4.3) and Uij = Vi + v' 3 , + 2, i = 1, . . . , m, j — 1, . . . , n. 

Divide (13.1) by the matching generation function of the region (13.2) and let N — > 00. 
By (3.3) and (13.3) we obtain 



L0 b 



Rl Rm . R[ R' n 

v 1 v m ' v[ v' n 



2 m+n Xrn,a f[ R t (R t ~ l/2)(i2< + 1/2) [] R[ 



j2 ■■■ J2 J2 ••• (_i)E£ 1 («i+bi)+E?=i(c i +d i ) 

0<ai<hi<fli 0<a m <b m <R m 0<c 1 <d 1 <R[ 0<c n <d n <R' n 

X 11 (2a, + 1)! (ifc - a t )\ (26, + 1)! (R t - h)\ 

-A- (fl t + c t - 1)1(^ + ^-1)! 
X 11 (2ci)! (ifc - Cj)! (2*)! (i?i - di)! 

x -pV (3/2) Ui+a< (3/2)„ i+6i -pr (3/2)„> +c< (3/2) <+di 
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x ] I (vj —Vi + dj - di)(vj —Vi + aj - bi)(vj —Vi + bj - ai)(vj —v% + bj - hi) 

l<i<j'<m 

X Y[ ( V 'j ~ V 'i + C - C i)( V 'j ~ V i + c 3 - d i)(v'j - v i + d j - c i)( v 'j ~ v i + d j - di) 
l<i<j<n 

XlZM-^) 2 K=M-dtf 



II™ i + a i + c j)( u ij +ai + dj)(uij +bi + Cj)(uij + hi + dj) 



(13.4) 





Ri 


R-m 








Vl 


v m 


Ri 




R m 


E 




E 


01,61= 





a m ,6 m = 






m 






*II 






1=1 

n 






*n 









where the summation range is restricted to those summation variables for which v\ + 
ai,vi + bi, . . . ,v m + a m ,v m + b mi as well as v[ + ci,v[ + di, . . . , v' n + c n ,v' n + d„, are 
distinct. 

By the argument in the last paragraph of the proof of Lemma 5.1, the restrictions 
a,i < bi, i = 1, ...,m, and c, < di, i = l,...,n can be dropped at the expense of a 
multiplicative factor of i/2 rn+n . We obtain the following result. 

Lemma 13.1. For fixed R\ , . . . , R m ,R[,..., R' n > 1 and v\ , . . . , v m , v[ , . . . , v' n > we 
have 



■ ■ V ) = X2m.2n ]J Ri(Ri ~ 1/2)^ + 1/2) ]J R[ 
n ' i=l i=l 

1[ R' n 

P ... ^ (-l)^™i( a '+ b ')+S?=i( c '+ d ') 
ii=0 c n ,d n — 

(Ri + a t - 1)! {Rj + h- 1)! (3/2)„ i+ni (3/2)„ i+6< 
\ (2a t + 1)! (ifc - a t )\ (26, + 1)! (R t - h)\ (2) Vi+ai (2) Vi+bi 

^ (R l + c t -l)\{R l + d l -l)\ (Z/2) <+Ci (Z/2) <+di 
J (2c,)! (ifc - Ci)! (2d,)! (ifc - *)! (1)„ J+C1 (l)„j +di 

x |[ (vj -Vi + aj - a l )(v j - Vi + a 3 - bi)(vj - Vi + bj - a{)(vj - v l + bj - b { ) 

l<i<j<m 

x f] (v'j -v'i + Cj - Ci)(v'j - vl + Cj - di)(v'j - v'i + dj - a)(v'j - v\ + dj - di) 

l<z<j<n 



UT=i( a i- b i) 2 Uti( c i-di) 2 



Y\T= i ]Tj=i( u ij + a i + c j)( u ij +ai + dj)(uij +bi + Cj)(uij + bi + dj) 



(13.5) 



where 



Uij = Vi + v'j + 2 



for i = 1, . . . , m, j = 1, . . . , n, and \ is given by (4.3). 

Note that expressions (13.5) and (5.1) are nearly identical. There are only two differ- 
ences: first, the roles of Ri and i?- are interchanged in the products on the first lines of their 

79 



right hand sides; and second, in the denominators on the third and fourth lines of (13.5) 
one has the expressions (2a^ + 1)!, (26j + 1)!, (2cj)! and (2di)\, while in the corresponding 
positions in (5.1) one has (2dj)!, (26j)!, (2cj + 1)! and (2di + 1)!, respectively. 

The fact that these are the only differences allows us to obtain the asymptotics of uib 
from the asymptotics of Wb worked out in Sections 6-12. 

The arguments of Section 6 apply equally well to yield an expression for u>b analogous 
to the expression (6.5) for u>b- By the the first two formulas in the proof of Proposition 6.1, 
besides interchanging the roles of Ri and R\ in the products on the first line of (5.1), the 
only effect of the difference between (13.5) and (5.1) is that the fractions 1/2 and 3/2 at 
the denominators of (6.7) and (6.8) swap places for the sums T(") and T' [n) — the analogs 
of T(") and T' (n) . More precisely, define 

Then by the arguments that proved Proposition 6.1 we obtain the following result. 
PROPOSITION 13.2. The boundary-influenced correlation H>b can be written as 

/ R 7? R> R' \ -^r " 

Qb u! • • • v ; v' • • • j = x ^ n r ^ V2)(^ + 1/2) n ^ 

Z— 1 2 — 1 



^ e(C , )M Ql(c . )j/ 3 l(c)i ... !Qm(c)!/ 3 m(c)m(c) ^ l(c) _. j7n(c) ^ n(c) 

CGC 



(13.8) 



where \ is given by (4.3), the collection C and e(C), a>i(C), Pi(C), Jj(C), Sj(C) are as in 
(6.2), and the "moment" M ai ,p 1 ,...,a m ./3 m ;-y 1 .S 1 ....,-y n ,8 n equals the 4mn-fold integral 

Ma 1 ,/? 1 ,...,a m ,/3 m ;7iA,-..,7nA, = / • / UY\(XijyijZi j W ij ) Vi+V i + 1 

Jo Jo 

n n 

x T^iRuw, J] x ljVlj ) T&XRuV!-, J] z ljWlj ) ■ ■ ■ 

3=1 3 = 1 

n n 

i=i i=i 

m m 

x T' (7l) „J ; [] a^) T' (5l) ^; [] i/a^i) • • • 

i=l i=l 

m m 

■ ■■T' hn \Rl l ,v' n ;Y[x m z m )f' iSn \R' n7 v' n: Y[y l nW ln )dx 11 ■ ■ ■ dw mn , 

»=i »=i (13.9) 
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withT (n \R,v;x) andT l{n) {R,v-x) defined by (13.6) and (13.7). 
Furthermore, the calculations that proved Lemma 6.2 show that 
T (n) (i?,v;x) = 



1 (3/2), 
R (2), 



k=0 



{-R) k (R) k (v + 3/2) k ( x\k 



(3/2) k (v + 2) k 



(!) 



3-^2 



-R- 



k , R -\- k , ^ 
| + fc,2 + <H 



T' v "'(i2,t;;a;) = 



1 (3/2), 



R (l)v 



E/< 

fe=0 



(-i?)fe (fl) fc (v + 3/2) fc /x 



(l/2)fc(u + l) fc 



(!) 



3-^2 



-i? + fc, i? + k, § + 1> 
\ + k,l + v + k 



k x 
5 4 



fc x 
5 4 



(13.10) 



(13.11) 

where the f k 's are as in (6.9) (in particular /„ = 1). These expressions are nearly the same 
as (6.10) and (6.11): the only difference is that, compared to and T" , the fractions 
1/2 and 3/2 at the denominators (including the denominator parameters of the 3F 2 's) 

swap places in the "barred" versions and f^ n \ 

Using (7.8) to write the 3-F 2 's of (13.10) as integrals of 2-F1 's, the derivation rule for 2-F1 's 
displayed after (7.9), and the 2-F1 evaluation (7.21), we obtain after simplifications 

f^(R,qR + c;x) = ^-J2^ / ^+ c+fe+1 / 2 (1 - i)^ 2 



2R 



d k 



dt k 4i? 2 - 1 



4 - xt 
xt 



■ sin 



i?arccos 1 



xt 
~2 



1 



AR 2 - 1 



cos 



i?arccos 1 



xt 
~2 



dt. 



(13.12) 



- (n) 

To obtain a similar expression for T' we need first an analog of (7.16). One readily sees 
that the calculation that proved (7.16) also shows that 

3 



3^2 



-R + k, R + k, § + v + k x 



± + k,l + v + k 
x {-R + k)(R + k) 



4 (1/2 + &)(!+« + A) 



4 

3F2 



3F2 
-R- 



R+k, R + k, ±+v + k x 
5 4 



ij + fc, l + v + k 



k + 1, R + k + 
§ +k, 2 + v 



1 2 

2 



; 4_ 

Express the 3-F 2 's of (13.11) using the above identity. Use (7.8) to write the resulting 
3^2 's as integrals of 2-Fi's. In turn, express the latter, using the derivation rule for 2-Fi's 
displayed after (7.9), in terms of the 2-F1 evaluation (7.11). We obtain after simplifications 

^fe+c-l/2 



T' {n) (R,qR + c; x) 



r) '» rl fk+c— 1/ 

n „i 



fc=0 

2 



/ 2 dt k 



cos 



i?arccos [ 1 — — 
2 



dt 



-J ^fc+c+l/2 



(1 - i)V2 



COS 
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i?arccos 1 — 



xt 



dt. 



(13.13) 



The asymptotics of the T^'s and T /(n) ' s can be obtained using the same approach we 
employed in Section 7 for the fW's and T'^'s. Indeed, by Lemma 7.4, when (7.23) is 
substituted in (13.12) the omitted terms in (7.23) give rise to integrals like in Proposi- 
tion 7.2. By the argument in the proof of Proposition 7.1, we deduce that the asymptotics 
of T( n \R, qR + c; x) is generated by the main term on the right hand side of (7.23) for 
k = n. More precisely, we obtain 

1 /■! f n+c+l/2 



X cos 



(l-t)V2 

i? arccos ( 1 - 



+ 



(n — 1)tt 



dt 



< M R n - 7/2 , 



(13.14) 

for R > Ro, with and Mo independent of x e [0,1]. Approximating the integral in 
(13.14) by Proposition 7.1 we deduce the first part of the following result. 

PROPOSITION 13.3. Let q > be a fixed rational number, and let n > and c be fixed 
integers. Then for any real number x e (0, 1], we have 



T (n \R,qR + c; x) 



1 



q 2 + 



4—x \l 4—x 



R 



X cos 



R arccos ( 1 — 



V 2/ 2 



arctan — 



X (n — 1)"7T 



< -—MR n ~ 7/2 



fi (n) (R, qR + c; x)-^{\q 2 



X cos 



R arccos 



(-1)4 



x R 1 / 2 



R 



x 



1 

arctan - 

q 



x 



A-x + ~ 



< M'iT~ 3/2 , 



(13.15) 



(13.16) 



for R > Ro, where Ro, M and M' are independent of x e (0, 1]. 

Proof. To prove (13.16), substitute (7.13) into (13.13). By Lemma 7.4, all omitted terms 
in (7.13) give rise to integrals of the type addressed by Proposition 7.2. By the arguments 
in the proof of Proposition 7.1, the asymptotics of T'^ n \R, qR + c; x) comes from the term 
on the right hand side of (7.13), when k = n. More precisely, we obtain 

n,(n) 



T' y ' L '(R,qR + c; x) 
\<lq 



vRJo 



xt 



xt 



cos 



R arccos 1 — 



xt 



+ 



mr 



-1 £0-1/2 



n+l 



xt 



A-xt 



cos 



U arccos ( 1 — ^ ) - 



< KR n - 3 / 2 , 
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for all R > p and x G (0, 1], where the constants K and p are independent of x € (0, 1]. 

The two integrals on the right hand side of (13.17) can be approximated by means 
of Proposition 7.1. Carrying this out and using the equation displayed after (7.26) one 
obtains (13.16). □ 

By Proposition 13.2, the analysis of the asymptotics of Ld b is a repeat of the analysis of 
the asymptotics of u b , with the TW's and T' (n) 's being replaced by the fW's and T' (n) 's, 
respectively. 

Comparing Proposition 13.3 with Proposition 7.1 shows that the only differences be- 
tween the formulas approximating the T^'s and T'^'s on the one hand and the ones 
approximating their unbarred versions on the other are: 

(i) The numerators 2 and 1 of the first fractions in the approximations are swapped in 

(13.15) (13.16) as compared to (7.1)-(7.2); 

(ii) The factor ^Jj^ m the denominator of the approximant of T'^ is moved to the 

denominator of the approximant of the "unprimed" ; 

(Hi) The powers of R in the approximants (13.15)— (13.16) are i?~ 5 / 2 and i? -1 / 2 , as 
opposed to both being R~ 3 / 2 in (7.1)-(7.2); 

(iv) Decrementation of the argument of the cosine by ir/2 occurs for in (13.15)- 

(13.16) , as opposed to occuring for T' (n) in (7.1) (7.2). 

Find the asymptotics of the moments M Qlji g 1) ... jamj| g m;7lj 5 1; ... j7n; 5 n by the same reasoning 
we used for the moments M in Sections 8 and 11. Relations (11.1) and (11.2) change 
slightly for the present case, to reflect differences (i)-(iv) above. Relation (11.3) stays 
unchanged. The resulting analog of (11.4) is 

E m ( R \ ai+l3i n ( R' ^ 7j+Sj 



Afai,ft,...,a m ,/J mi 7i,«i,...,7».«» ~ ^i^T II (~^J II 




Ai(e 2t -i + e 2i ) + ^2 Bj(e 2m +2j-i + e2m+2j) 
i=i j=i 

e2m+2j - 1 + £2m+2j arctan- 1 



z}V5 

+ ^2(e2z-ia t + e 2i Pi)^ + ^(e 2 m+2j-i7j + e 2m +2jS j )'^ - ^(e 2 j-i + e 2j )^ 
i=i j=i j=i 



m n 



/ £2i-iAj + e 2m+2 j-\Bj e 2i -iAi + e 2m+2 jBj 
I arctan r-— -(- arctan r-— 



- y^y^ arctan- , .... 

arctan ^^+6 2m +2 arctan 1 
+ 0(i? _4mn +S™ i("i+ft)+E" = i(7j+<5j)-3m-3n-l^ (13.18) 
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where 



E = 



2 + i 
^ 3 



* nr =1 i 



(13.19) 



and D«i.-.e2m+2n i s gi ven by (11.6). 

All consequences of differences (i) and (ii) are reflected in the change of exponents of 
2 and 3 in (13.19) versus (11.5). All consequences of difference (Hi) are reflected in the 
changed exponents of Ri and R\ in (13.19) versus (11.5). The combined effect of differences 
(i)-(iii) on (11.4) is thus a multiplicative factor independent of a\, . . . , j n , S n . 

The only effect of difference (iv) is that the last sum on the fourth line of (13.18) is 
E"=i( e 2j-i+ e 2j)f , as opposed to ^™ = i(e2m+2j-i+e2m+2j)f in (H-4). Furthermore, the 
statement of Lemma 11.1 is clearly valid also for the multiple sum on the right hand side of 
(13.18). Therefore, in the multiple sum of (13.18) summation can be restricted to balanced 
(ei, . . . , £2m+2n)'s without changing its value. However, this restriction eliminates the only 
difference between the multiple sums in (13.18) and (11.4). By the arguments that proved 
(11.22) wc obtain 



Ri R m R 1 R' n 



V m V\ 



XlmJln 



\{R i (R i -l/2)(R i + l/2)J{R' j 



i=l 



3 = 1 



i 

lit i n "=i + v'j -oi- Jj) 2 + \{Ri- R'j) 2 ] [(vi + v> -a- cj) 2 + k(Ri + R'if 



i=i j=i 



>< n 

l<2<j<m 



(v j -v i f + -{R j -R i y 



1 



>< n 

l<z<j<n 

->2m 2 +2n 2 -4mn-2m-l 



K-^) 2 + 3^-^) 2 



(v j -v i ) 2 + -(R j +R i ) 2 



[v' i -v' i ) 2 + -{R' j + R' i ) 2 



+ 0(R 2 



(13.20) 



This leads to the following result. 

Proposition 13.4. Let the parameters Ri, . . . , R m , vi,...,v m and R[, . . . , R' n , v[,...,v' n 



depend on R as in (2.3). Then as R — ► oo the asymptotics of Wb is given by 



x II [(^--Ri) 2 + 3(^-Wi) 2 ][(^+-Ri) 2 +3(^-^) 2 

l<i<j<m 



>i-Vi) 2 ][{Rj 

" ut i n"=i[(^ - + 3^- + v^} [{X. + Rtf + 3(4 + ^) 2 ] 

+ 0( J R 2m2+2 " 2 - 4m "- 2m - 1 ), (13.21) 

w/iere 

2 , 3 -(| C -l)V4 + (3 fc -l)/4 *-l (2) ^ ( J+2)fc 

- — ^ — n mm* }i 7372Y' (13 - 22) 

Proof. Since by (2.3) the parameters . . . , i? m , «i, . . . , v m and R[, . . . , R' n , v[, . . . ,v' n 
approach infinity as R — > oo, while c±, . . . , c m and c' 1; . . . , c' n are constant, it follows that 
the difference between the product on the right hand side of (13.20) and what it becomes 
when one omits the additive constants ±1/2 on its first line and the constants Cj and c'j 

is o(R 2m2+2n2 - imn - 2m - 1 ). This proves (13.21). Since by assumption the pairs (A l ,q l ), 
i = 1, . . . , m of (2.3) are distinct, as well as the pairs (Bj, q'j), j = 1, . . . , n of (2.3), the 
expression on the right hand side of (13.21) has degree 2m 2 + 2n 2 — 4mn — 2m in R. 
Therefore, (13.21) does indeed give the asymptotics of u>b- □ 
The proof of Theorem 2.1 now follows readily. 

Proof of Theorem 2.1. The statement follows directly from Proposition 3.2, using The- 
orem 2.2 and Proposition 13.4. □ 

14. A conjectured general two dimensional Superposition Principle 

The choice of the side-lengths of the regions (2.1) might seem unmotivated at first, 
but is in fact quite natural. Indeed, one readily sees that for any lattice hexagon on the 
triangular lattice, the difference between the lengths of all pairs of opposite sides is the 
same. Furthermore, this common difference is equal to the difference between the number 
of up-pointing and down-pointing unit triangles enclosed by the hexagon. Therefore, if 
we want to choose a lattice hexagon H to contain all our plurimers and possess dimcr 
coverings after their removal, the difference between opposite sides of H has to match the 
difference 4n+ 1 — 4m between the number of up-pointing and down-pointing unit triangles 
in the union of the plurimers. It follows that the sides of H must have the form indicated 
in Section 2. 

By a result of Cohn, Larsen and Propp [8] , a dimer covering of a large regular hexagon 
sampled according to the uniform distribution on all dimer coverings has maximal entropy 

85 



statistics at its center (and only at its center). This suggests that it is natural to place the 
hexagon H enclosing the plurimers so that they stay at its center when one scales and lets 
the size of H grow to infinity. By this observation and the previous paragraph, it follows 
that from this point of view, up to a translation by an absolute constant (independent 
of the size of H), (2.1) arc the most natural regions to use in the definition (2.2) of 
the plurimer correlation. (We note that translations by a constant vector that keep the 
symmetry about I can be handled by exactly the same approach we presented in Sections 
2-14.) 

In the results of this paper we assumed that all plurimers are triangular. Further- 
more, we assumed that all plurimers, except for a single monomer u, have side-length 2, 
are distributed symmetrically about the say vertical symmetry axis I of u, and that the 
two plurimer orientations are separated by any horizontal h through u, i.e., all plurimers 
above h point upward, and all below h point downward. However, we conjecture that the 
Superposition Principle (2.6) is valid for arbitrary plurimers. 

To make this conjecture precise, fix a vertex O of the triangular lattice. Consider a 
rectangular system of coordinates centered at O, with the unit on the horizontal axis equal 
to the side of a unit triangle, and the unit on the vertical axis equal to twice the height of 
a unit triangle (these units are chosen so as to maintain consistency with Section 2). 

Let Pi, . . . , P n be arbitrary plurimers on the triangular lattice (i.e., finite unions of unit 
triangles), and let pi be a distinguished lattice point (a base point) of Pi, i = l,...,n. 
Write Pi(Ri,Vi) for the translation of the plurimer Pi that takes the base point to the 
lattice point of coordinates (Ri,Vi). 

As in Section 2, let the charge ch(P) of the plurimer P be the number of the up-pointing 
unit triangles of P minus the number of its down-pointing ones. Let Hn be the lattice 
hexagon centered at O and having side- lengths alternating between TV and N — k, where 
k = ch(Pi) + • • • + ch(P„). Define the correlation u(Pi(Ri, vi), . . . , P n {Rn,v n )) by 

t r> I r, \ r> I r> \\ v M(Hn \ Pi (i?l , Vl ) U ■ ■ ■ U P n (R 

.(Pi(Ri,vi), P n (R n , Vn} ) = ^hm M{HNXPi{aiM)u ... uPn{aM y (^D 

where ai and hi, i= 1, . . . , n, are some fixed integers specifying a reference position of the 
plurimers. 

We conjecture that the following generalization of Theorem 2.1 holds. 

Conjecture 14.1. Suppose the coordinates Ri and Vi, i = 1, . . . , n are expressed in terms 
of the integer parameter R as 

Ri = AiR 

Vi = q.R, (14.2) 



where < A. t E Q, < q t E Q, i = 1, . . . , n are all fixed. 
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Then the asymptotics of the plurimer correlation is given by 



w(P 1 (R 1 ,v 1 ), P n (R n , «„)) = c H [(Rj - R,f + 3( Vj - ^) 2 ] ch(Pi) ch «>/ 4 

l<i<j<n 

+ 0(JjEi< i<J <„ch(P i )ch(P ;j )/2-l^ 

= C n d(mR l ,v t ),p 3 (R 3 ,v 3 )) ch(Pi)ch{P3)/2 

l<i<j<n 

+ 0( jR Ei< i<3 <nCh(Pi)ch(P j )/2-l- ); 

'(14.3) 

where d is t/ie Euclidean distance, and c depends just on the shapes of the plurimers 
Pi, . . . , P n , and not on their coordinates (Ri, v{). 

Writing, for the sake of notational brevity, Pi = Pi(Ri,Vi), it follows from (14.3) that 
the correlation u satisfies 

u(P u P 2 ) = c d(P u P 2 ) ch ^ ch ^/ 2 + OiR^^^^ 2 - 1 ) (14.4) 

U (P 1 ,...,P n )=cf H wiPuPj) +0(RZ i<i<i<» cMP^chfP,)^-!^ 

l<i<j<n (14.5) 

where c' is a constant independent of R. When taking the logarithm in (14.4), the contri- 
bution of the constant is negligible as R — > oo and we obtain 

In^fi.^- ^^ d^,^), P-oo. (14.6) 

Similarly, provided ^ Kj< . <n ch(Pj) ch(Pj) ^ 0, when taking the logarithm in (14.5) one 
can neglect the contribution of the constant and we obtain 

Inw(Pi,...,P„)~ Yl ^( p i,Pj), R^^- (14.7) 

l<i<j<n 

Equations (14.6) and (14.7) show that In a; satisfies the characterizing properties of the 
two-dimensional electrostatic potential — Coulomb's law and the Superposition Principle. 
Since all classical electrostatics can be deduced from these two properties (see e.g. [11, Ch. 
4]), Conjecture 14.1 implies that our random tiling model described in Section 2 indeed 
models classical two-dimensional electrostatics. 

We believe that in fact (14.3) holds in still larger generality. There are two ingredients 
to this extension. First, allow the family of regions used to define the plurimer correlation 
to be any family with the property that in the scaling limit the plurimers are situated in 
the region where dimcr coverings have maximal entropy statistics. We believe that the 
correlation defined by means of any such family of regions satisfies (14.3). One instance 
of this situation is presented in [6], where it is shown that (14.3) holds for two charges of 
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magnitudes 2 and —2. An extension found by the author of the result in [6] to an arbitrary 
number of even-side plurimers will be presented in a sequel of the present paper. 

Second, based on our result (14.9) below and the conjectured rotational invariance 
of monomer-monomer correlations (sec [12]), we conjecture that the natural analog of 
Conjecture 14.1 on the square lattice 1? also holds. 

More precisely, let P\,...,P n be arbitrary plurimers on the square lattice (i.e., finite 
unions of unit squares). Fix a chessboard coloring of the square lattice, and define the 
charge ch(P) of plurimcr P to be the difference between the number of its white and black 
unit squares. As we did for the triangular lattice, consider a base lattice point pi in each 
Pi, and denote by Pi(xi, ?/,) the translation of Pi taking to (xi,yi) (in order for Pi(xi, yi) 
to "look like" Pi , we assume all such translations to be color preserving) . 

Let ARn be the "Aztec rectangle region" of sides N and N + k centered at the origin 
(i.e., the lattice region dual to the corresponding Aztec rectangle graph defined in [3]), 
where k — ch(Pi) + • • • + ch(P„). Define the correlation of the n plurimers by 

, p , \ P / n i- M(AR N \ Pi(xi,yi) U ■ ■ ■ U P n (x n ,y n )) 

u{Pi{xi,yi),...,P n (x n ,y n )) = hm — — — — , 

n^oo M(AR N \ Pi(oi,6i) U • • • U P n (a n ,b n )) 

where a, and 6j, i — 1, . . . , n, are some fixed integers specifying a reference position of the 
plurimers. 

Conjecture 14.2. Suppose the coordinates x t and y i; i — 1, . . . , n are expressed in terms 
of the integer parameter R as 

x% — A^R 

y% = QiXi, 

where < £ Q, < Pj £ Q, i = 1, . . . , n are all fixed. 

Then the asymptotics of the plurimer correlation is given by 

W (Pl (*!,»!), . . .,P n (x n ,y n )) - C [] ~ X *f + (Vi - y,) 2 ] Ch(Pl)Ch(Pj)/4 

+ 0(JjEi< 4<j <„ch(F 4 )ch(F J )/2-l^ 

= c [] ^Pi{x i ,y i ),Pi{x j ,y j ))^ p ^ p ^ 2 

+ 0(ij£i<*<i<» ch ( F ') ch ( f i)/ 2 - 1 ), 

'(14.8) 

where c is a constant depending just on the types of the plurimers Pi , . . . , P n . 

More generally, we believe that any bipartite planar, periodic (i.e., invariant under 
translations by two non-collinear vectors) graph has an embedding so that joint corre- 
lations of plurimers are given by Coulomb's law (using the Euclidean distance) and the 
Superposition Principle. 



For the square lattice, the special case of (14.8) when n — 2 and the two plurimcrs 
are in fact monomers of opposite color was suggested to hold by Fisher and Stephenson 
in [12]. The further specialization when the second monomer is adjacent to a lattice 
diagonal through the first was proved by Hartwig [17]. A natural extension of Hartwig's 
result would be to study the correlation of an arbitrary collection of monomers along two 
consecutive lattice diagonals. We state below (see (14.9)) a result we found which is a 
close analog of this. The proof will appear elsewhere. 

For this purpose, rather than phrasing everything in terms of lattice regions and their 
dimer (domino) coverings, it will be convenient to use the dual set-up of graphs (duals of 
regions) and their perfect matchings. 

Consider therefore the square lattice 1? and regard it as a graph. It will be convenient 
to draw it so that the lattice lines form angles of 45° with the horizontal. 

We say that a vertex v of the square lattice has been split if v is replaced by two 
new vertices v' and v" slightly to the left and right of v, respectively, and the four edges 
incident to v arc replaced by two edges joining v 1 to the former two neighbors of v on its 
left, and two joining v" to the former two neighbors of v on its right. Figure 14.1 contains 
an example of a split vertex. 

Color the vertices of the square lattice 1? black and white in a chessboard fashion. As 
above, regard removal of a vertex of 1? as creation of a unit charge, of sign determined 
by its color. In addition, it is natural to consider the operation of splitting a vertex as 
creating a unit charge of opposite sign to the one that would be created by removing that 
vertex (this is readily seen to be justified if one views perfect matchings as being encoded 
by families of non-intersecting lattice paths). 

Let to, n > be integers, {vi, . . . , v m } and {wi, . . . , io„] two disjoint sets of nonnegative 
integers, and N an even nonnegative integer. We define ARn(vi, . . . , v m ;wi, . . . , w n ) to 
be the subgraph of Z 2 described as follows. 

Consider the Aztec rectangle (see [3]) of width N and height N + m — n, i.e., a subgraph 
of 1? consisting of an N by N + to — n array of 4-cycles touching only at vertices. Let I 
be its vertical symmetry axis, and let O be its vertex on i that is N lattice segments away 
from its base. Label the vertices on I starting with for O and continuing with consecutive 
integers as we proceed upward. We define AR^(vi, . . . , v m ;w\, . . . , w n ) to be the graph 
obtained from our Aztec rectangle by removing the vertices labeled Vi, i = 1, . . . , to, and 
by splitting the vertices Wi, i = 1, . . . , n (an example is shown in Figure 14.1). 

As mentioned above, we regard both kinds of altered vertices as unit charges. We take 
the reference position of these charges to be when they are packed next to O, and thus 
define the joint correlation of the charges as 



u(vt, . . .,v m ;wi, . . .,w n ) := 

lim MjARNjv!, ...,ti m ;wi,.. .,w n )) 

Af^oo M(AR N (0, 1, . . . , to - 1; to, to + 1, . . . , to + n - 1)) ' 



Using results from [3] and doing the asymptotic analysis of the expressions they lead to, 
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as the distances between charges grow to infinity keeping constant mutual ratios, where 
c TO; „ is some explicit constant depending only on m and n. The above formula shows that 
the Superposition Principle holds on the square lattice for any distribution of unit charges 
along a lattice diagonal. 

As a particular case of (14.9), we obtain, after completing the fairly laborious task of 
working out the constant, that 

"M-j^TeA d^oo (14.10) 



where A = 1.282427... is the Gleisher-Kinkclin constant [14], 



1 1 2 2 ■ ■ ■ n n ~ ^4 n ™ 2 /2+"/2+l/12 e -n 2 /4_ 



Formula (14.10) is a counterpart of the result of Hartwig [17], which addresses the case of 
vertices of opposite colors removed from adjacent diagonals. 

Remark 14.2. After scaling, the above set-up places the momomers whose correla- 
tion is measured exactly in the center of the scaled Aztec rectangle, which in the limit 
approaches an Aztec diamond. According to a result of Cohn, Elkies and Propp [7], a 
uniformly sampled tiling of an Aztec diamond has maximal entropy statistics at its center, 
and only at its center. The above set-up seems therefore natural. 
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15. Three dimensions and concluding remarks 



As far as the author knows, the question of working out the joint correlations of an arbi- 
trary number of holes for the dimer model on a bipartite lattice — and thereby establishing 
the emergence of electrostatics in this way — was not considered before in the literature. 
Perhaps the case of removing several regions from a lattice and study their joint correla- 
tions was not considered before due to the complexity the problem presents already in the 
case of two removed monomers (see [12]). This might have been also the reason for not 
considering the case of two removed regions of arbitrary charges, a situation that already 
points to the phenomena of electrostatics. Another possible reason is that the analysis 
of [12] was done on the square lattice, while the hexagonal lattice is the one on which 
natural regions to remove — triangular regions — can have arbitrary charges. (Indeed, on 
the square lattice to get a region of charge s one needs for instance to line up s diagonally 
adjacent monomers. Removing a whole dimer is natural on the square lattice, and this 
was studied extensively in [12]. However, since dimers have charge zero, the study of their 
correlations does not detect the Superposition Principle!) 

We mention that there are results in the physics literature outlining some connections 
of the Ising model, a well studied statistical physical model, and electrostatics (e.g. Hurst 
and Green note in [18] that the relations satisfied by certain matrices they employ to solve 
the Ising model "are similar to the relations satisfied by fermion emission and absorption 
matrices," and in [30], which presents exact calculations of 2-point spin-spin correlations 
in the Ising model, Wu, McCoy, Tracy and Barouch are led to certain integral equations 
that arose before in work of Myers [26] in the context of electromagnetic scattering from 
a strip; in [25] McCoy, Perk and Wu give recurrences for n-point spin correlations, but 
without analyzing their asymptotics). Since the Ising model on a lattice can be equivalcntly 
phrased as a dimer model on a modified, suitably weighted lattice (see e.g. [19]), this yields 
a connection between electrostatics and the dimer model. However, this turns out to lead 
to a totally different object of study in the dimer model than the joint correlation of holes 
we considered in this paper: two spins being correlated in an Ising state turns out to 
correspond to requiring two faces in the corresponding dimer lattice to be on the same 
side of the union of certain polygons corresponding to the Ising state. Moreover, the 
fundamental difference (4) mentioned in the Introduction between our model and the ones 
surveyed in [27] holds also between our model and the Ising model. 

One key feature of the work described above is that even though the quantities that 
need to be studied asymptotically are not "round" (i.e., do not factor as products of small 
factors) in general, they can be expressed as multiple sums of round terms, courtesy of the 
exact enumeration formulas from [5] and [3]. These multiple sums can then, in certain 
situations, be reduced to single or double sums that in turn lead to special functions 
whose asymptotics can be obtained using specific techniques such as Laplace's method. 
This view supplies additional motivation for the already well-represented study of regions 
whose tilings are enumerated by simple product formulas. 

The question of studying joint correlations of missing vertices in lattice graphs can 
naturally be phrased also in dimensions other than two. The really compelling case is that 
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of three dimensions. Does there exist a three dimensional analog of our model that would 
model classical electrostatics in the physical, three dimensional space? 

We present such a possibility in Question 15.1 below. In the case of an affirmative 
answer to this question, we show in Remark 15.2 how a natural new parameter could be 
introduced in our model so that the parallel between the three dimensional analogs of 
(2.13) and (2.14) holds for any temperature T. 

Based on the two dimensional case, we are guided in the phrasing of (15.2) by the 
assumption that in three dimensions as well an analogous random tiling model would 
manifest the phenomena of electrostatics. 

The square lattice version of our model is easiest to parallel in three dimensions. Con- 
sider R 3 divided into unit cubes by the lattice Z 3 . Color the unit cubes black and white so 
that adjacent cubes have different colors. Regard the unit cubes as monomers, and their 
finite unions as plurimers. 

Consider two monomers a = (0, 0, 0) and b r = (0, 0, 2r — 1) of opposite color and include 
them in a large cube Cn, where N is even. Based on our guiding physical intuition, we 
wish to define their correlation co(a, b r ) to be proportional to M(CV \ {a, b r }), for large N, 
in such a way that — lnw(a,6 r ) behaves asymptotically like the potential energy of two 
unit charges of opposite sign at distance 2r, i.e., like a positive constant times — l/(2r). 
In particular, we should have ui{a, b r ) — > 1, as r — ► oo. 

This suggests that in the definition of u(a,b r ) we should normalize by letting the 
monomers be far apart. Based on this we define plurimer correlation as follows. 

Fix a lattice point O, and consider a rectangular system of coordinates centered at O. 
Let Pi be an arbitrary plurimer, and let pi e Z 3 be a base point of Pi, for i = 1, . . . , n. 
For any integer R>1, define RPi to be the translation of Pi that takes p» to Rpi. 

Enclose the n plurimers in a large cube Cn of side N, centered at O. For the sake of 
definition simplicity, assume that ch(Pi) + • • ■ + ch(P„) = (this condition is not essential 
for defining the correlation below and phrasing Question 15.1 — it can be circumvented by 
replacing the enclosing cubes Cat by a suitable family of regions Dn with the property 
that ch(Djv) = ch(Pi) + • • • + ch(P„); for 1? , this was accomplished in Section 14 by 
choosing Aztec rectangles AR N ). 

Assuming the following limit exists, define the plurimer correlation oj(Pi, . . . , P n ) by 

ip r r M(Cat \ P 1 U ■ ■ ■ U P n ) 
uj(Pi, . . . , P n = hm hm — — — ; -, 15.1 

V ' ' ' R^ooN^ooM(C N \RP 1 U---URP n ) ' 

where M{D) denotes the number of dimer coverings of the lattice region (on Z 3 ) D. (The 
inside limit should exist for any R — it is just the ratio of two joint correlations of plurimers 
in a sea of dimers; we are assuming in addition that its value approaches a limit as R — > oo.) 

The possibility we referred to above for modeling three dimensional electrostatics is 
phrased below. As in two dimensions, the charge ch(P) of a plurimer P is defined to be 
the difference between the numbers of its white and black unit cubes. 

Question 15.1. Is it true that 

co(RP u ...,RP n ) = l-k 3 ]T C ^ Ch ^ +0(R- 2 ) (15.2) 
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as R — > oo, where k 3 > is an absolute constant? 



If the answer to this question is affirmative, taking the logarithm in (15.2) we obtain 

lnu(RP 1 ,...,RP n )~-k 3 ]T C , h /^ )Ch i P W R ^°°- ( 15 - 3 ) 

^ d(Rpi,Rpj) 



l<i<j<n 



This would show then that indeed In a; behaves like the electrostatic potential in three 
dimensions. 

Remark 15.2. Assuming the answer to Question 15.1 is affirmative, we can introduce a 
new parameter in our model that allows it to parallel electrostatics for any temperature. 

The new parameter, denoted by x, ranges over the odd positive integers. For any such 
x, refine each unit cube c of Z 3 into x 3 equal smaller cubes, and properly color the latter 
black and white so that the smaller cubes fitting in the corners of c have the same color 
as c (this can clearly be done, since x is odd). Moreover, if we regard the subdivided cube 
c as a plurimer on the lattice (jZ) 3 , its charge is readily seen to agree with the charge of 
c, regarded as a monomer on Z 3 (again, this is due to x being odd). 

Let Pi and pj, i = 1, . . . , n, be as above. Define the correlation u) x (P\, . . . , P n ) by 

m P x ,. ,. M x (C N \P 1 U---UP n ) 

^ (Pl ' ' • ' ' ^ = i^c a!™ M^^NflPxU-UflPn) ' (15 ' 4) 

where for a lattice region H in Z 3 , M X (H) denotes the number of its dimer coverings when 
regarded as a lattice region in ( jZ) 3 . (Since as noted in the previous paragraph our lattice 
refinement preserves charge, H has the same number of black and white fundamental 
regions when considered in (^Z) 3 , provided it does so when considered in Z 3 .) Clearly, 
existence of the limit (15.1) implies existence of the above limit. 

The asymptotics of u> x (Pi, . . . , P n ) can be deduced from (15.1) as follows. Scale down 
the lattice Z 3 by a linear factor of x, keeping our lattice point O fixed, and view the 
plurimers and their base points as being in the scaled down lattice: the plurimers become 
Ff,i = l,...,n, and their base points pf = xpi, i = 1, . . . , n. It follows from the definition 
that 

w x (Pi,...,P n )=w(i?,...,0 

Clearly, we also have d(i?pf , i?pj) = xd(Rpi, Rpj), for all 1 < i < j < n. Furthermore, 
the plurimer (RPi) x is the same as the plurimer RPf. Therefore, we obtain by (15.2) that 

Taking the logarithm in both sides we deduce from (15.5) that 



l<i<j<n 
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Repeating the analysis at the end of Section 2, (15.6) implies that in the limit N — > oo, 
when one samples uniformly at random from all dimer coverings of the x-subdivided 
Cn \ Pi U • • • U P n , the relative probability of having the plurimers at preassigned distances 
dij versus is, by the fundamental theorem of statistical mechanics (sec e.g. [10, §40)]), 
exactly the same as the relative probability of having n corresponding electrical charges at 
those distances, at temperature 

T = xql/(4TTeokk 3 Lp), (15.7) 

where q e is the charge of the electron, k is Boltzmann's constant, Lp is the Planck length 
(<~ 1CP 35 , smallest non-zero length that "makes sense" physically; in our analysis above, 
we express all physical distances using L P as a unit 10 ) and e is the permittivity of empty 
space. 

The free energy per unit volume in this model equals the free energy per Z 3 -site in the 
dimer model on the lattice (j^) 3 . This energy is F — l^x 3 , where l 3 <~ 0.446. . . is the 
three-dimensional dimer constant (see e.g. [16] [4]). 

As mentioned in the footnote at the end of Section 2, there is an extra "calibration" 
parameter one can consider in our model — a fixed positive integer a so that a physical 
elementary charge corresponds to a plurimcr of charge a in our model. The effect of this 
extra parameter in (15.7) is to multiply its right hand side by 1/a 2 . Expressing x from 
the resulting relation and substituting into F = Z 3 x 3 we obtain 

f = ^kh Lpa y hT3 (158) 

This is similar in form to the formula for the free energy per unit volume, for high tem- 
perature T, obtained in quantum field theory. Note that the constant of proportionality 
in (15.8) is <~ (10 _30 a 2 ) 3 , so the value of F depends crucially on a; for small a the value 
is negligible even for very large temperatures T, while for a ~ 10 15 — a conceivable value, 
given that Lp ~ 10 -35 — it becomes quite significant. The free energy per unit volume F 
in our model is reminiscent of the cosmological constant, which also encodes in some sense 
the overall energy stored in the vacuum of space. The counterpart in our model of the 
fact that the cosmological constant is believed to be positive but too small to be detected 
by current experiments would then be that the calibration parameter a satisfies a <C 10 15 . 

It is important to remark that while the Superposition Principle seems to hold indepen- 
dently of the background bipartite lattice, the effect of lattice refinement on both (2.13) 
and its three dimensional analog that (15.2) would imply depends crucially on the lattice. 
For instance, (14.8) would imply that (2.13) holds also on the square lattice. However, in 
this case x-fold refinement does not provide an extra parameter: for even x the plurimers 
become neutral on the refined lattice, while for odd x the plurimer charges are invariant 

10 Strictly speaking, one needs to express physical distances in terms of Lp also in (2.13), but since 
(2.13) is invariant under changing the unit for distance, this was not necessary in the two dimensional 
case. 
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under refinement, a fact which together with the invariance of (2.13) under distance scaling 
shows that (2.13) does not change in this case under refinement. Similarly, for bipartite 
lattices in three dimensions the charge is not invariant in general under refinement, so 
different lattices lead to different analogs of (15.7) and (15.8). 

An interesting question is the following: if bipartite lattices lead to electrostatic effects, 
what do non-bipartite lattices lead to? Since the latter do not have a black and white 
coloring, it is natural to expect uniform behavior — either universal attraction, or universal 
repelling. Is one perhaps led to the effects of some other fundamental physical force? 
The most natural non-bipartite plane lattice, the triangular lattice, is considered in [9], 
and an analysis of the monomer- monomer correlations on it, paralleling that of Fisher 
and Stephenson [12] on the square lattice, is carried out. The first 15 correlations of two 
monomers in adjacent lattice lines are computed, and their values [9, Table III] provide 
convincing evidence for an exponential decay to a limiting positive constant. Changing 
the normalization in the definition of correlation in [9, p. 8] so that one normalizes by the 
reference position of two infinitely separated monomers (as in (15.4)), this is equivalent 
to the modified correlation exponentially approaching zero. Then the "potential" would 
be In uj ~ — kr, with k > a constant and r the separation between monomers, and the 
"field" hiw — or more directly, according to "discrete calculus" (see e.g. [15]), u(r + 
l)/w(r) — 1 — would approach — k. This is similar to the case of inter-quark force, which 
approaches a constant as the separation between the quarks approaches infinity. It would 
be very interesting to study the joint correlation of several plurimers on the triangular 
lattice in more detail, in particular to determine whether it is a function depending just 
on the pair correlations (as it is the case in the presence of a Superposition Principle), and 
to find out how it depends on the plurimer sizes. 

We conclude by mentioning that since the set of all tilings of a region is equivalent to 
the set of all families of non-intersecting lattice paths with certain starting and ending 
points — which is in turn cquinumerous, by the Gessel-Viennot theorem, with the set of 
all (appropriately signed) families of lattice paths with these starting and ending points — 
averaging over all such tilings is reminiscent of the "sum-over-paths" interpretation of 
particles in quantum mechanics due to Feynman. The results in this paper show that 
if the different tilings of a region with holes are associated with different ways for pairs 
of virtual particles and their antimatter companions to annihilate, then the microscopic 
frenzy of quantum-mechanical fluctuations in the vacuum of empty space that Feynman 
once jokingly described as "Created and annihilated, created and annihilated — what a 
waste of time," is seen to actually generate the phenomena of electrostatics 11 . 



11 A consequence of quantum mechanics is that the quantum fluctuations of the vacuum drive the 
intrinsic strength of the electric field of charged particles to get larger when examined on short distance 
scales. This offers the opportunity to perform a test for the parallel between our model and electrostatics. 
Explicit numerical calculations based on results of this paper and [6] confirm that the quantity u)(r + 
l)/u>(r) — 1 (where ui(r) is the correlation between a fixed hole and another hole of fixed shape at distance 
r), the analog in our model of the electric field, when divided by the two dimensional electric field intensity 
1/r, does indeed yield a quantity that grows as r — > 0. Details will appear in a sequel to this paper. 
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